Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,136,706 G→A 5.4% G145S (GGC→AGC)  flgI → putative flagellar basal body protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,136,7060GA 5.4% 235.9 / 2.0 74G145S (GGC→AGC) flgIputative flagellar basal body protein
Reads supporting (aligned to +/- strand):  ref base G (32/38);  new base A (2/2);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

TTTGGACGTCAGGGGCAAACCATCGATGTGGTGGTTTCTTCCATGGGAAATGCCAAAAGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGA  >  NC_000913/1136559‑1136851
                                                                                                                                                   |                                                                                                                                                 
ttgggACGTCAGGGGCAAACCATCGATGTGGTGGTTTCTTCCATGGGAAATGCCAATAGCTTGCGTGGAGGTACGTTGTTGATGACATCGCTTAAGGGCGTTGACAGTCAGGTGGATGCGCTGGCGCAGGGCAATATTCTGGTTggcgg                                                                                                                                                  <  2:88801/146‑1 (MQ=255)
   ggACGTCAGGGGCAAACCATCGATGTGGTGGTTTCTTCCATGGGAAATGCCAAAAGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGgcggcg                                                                                                                                                >  6:184865/1‑148 (MQ=255)
    gACGTCAGGGGCAAACCATCGATGTGGTGGTTTCTTCCATGGGAAATGCCAAAAGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTggcg                                                                                                                                                   >  2:257547/1‑144 (MQ=255)
    gACGTCAGGGGCAAACCATCGATGTGGTGGTTTCTTCCATGGGAAATGCCAAAAGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTggcg                                                                                                                                                   <  1:257547/144‑1 (MQ=255)
    gACGTCAGGGGCAAACCATCGATGTGGTGGTTTCTTCCATGGGAAATGCCAAAAGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCa                                                                                                                                              <  2:289740/149‑1 (MQ=255)
       gTCAGGGGCAAACCATCGATGTGGTGGTTTCTTACATGGGAAATGCCAAAAGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTTGCGGCGCAg                                                                                                                                             >  1:88216/1‑147 (MQ=255)
         cAGGGGCAAACCATCGATGTGGTGGTTTCTTCCATGGGAAATGCCAAAAGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGCCAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGgcgc                                                                                                                                         >  8:77278/1‑149 (MQ=255)
         cAGGGGCAAACCATCGATGTGGTGGTTTCTTCCATGGGAAATGCCAAAAGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGc                                                                                                                                         >  2:168450/1‑149 (MQ=255)
              gCAAACCATCGATGTGGTGGTTTCTTCCATGGGAAATGCCAAAAGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCg                                                                                                                                    >  6:143865/1‑149 (MQ=255)
                aacccATCGATGTGGTGGTTTCTTCCATGGGAAATGCCAAAAGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCt                                                                                                                                  <  1:291637/146‑1 (MQ=255)
                  aCCATCGATGTGGTGGTTTCTTCCATGGGAAATGCCAAAAGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTgg                                                                                                                                >  3:171902/1‑149 (MQ=255)
                   ccATCGATGTGGTGGCTTCTTCCATGGGAAATGCCAAAAGCTTGCGTGGAGGGACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGc                                                                                                                               <  3:292981/149‑1 (MQ=255)
                    cATCGATGTGGTGGTTTCTTCCATGGGAAATGCCAAAAGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGa                                                                                                                                           <  4:217180/136‑1 (MQ=255)
                    cATCGATGTGGTGGTTTCTTCCATGGGAAATGCCAAAAGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGa                                                                                                                                           >  3:217180/1‑136 (MQ=255)
                    cATCGATGTGGTGGTTTCTTCCATGGGAAATGCCAAAAGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCt                                                                                                                                  <  6:368191/145‑1 (MQ=255)
                    cATCGATGTGGTGGTTTCTTCCATGGGAAATGCCAAAAGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCt                                                                                                                                  >  5:368191/1‑145 (MQ=255)
                    cATCGATGTGGTGGTTTCTTCCATGGGAAATGCCAAAAGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGGCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCg                                                                                                                              <  7:238236/149‑1 (MQ=255)
                            tggtggTTTCTTCCATGGGAAATGCCAAAAGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAgt                                                                                                                      >  6:263712/1‑149 (MQ=255)
                                ggtTTCTTCCATGGGAAATGCCAAAAGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGtt                                                                                                                   >  6:252450/1‑148 (MQ=255)
                                        ccATGGGAAATGCCAAAAGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGtt                                                                                                             <  7:299041/146‑1 (MQ=255)
                                         cATGGGAAATGCCAAAAGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAAc                                                                                                          <  8:70671/148‑1 (MQ=255)
                                           tGGGAAATGCCAAAAGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCaa                                                                                                       >  4:102512/1‑149 (MQ=255)
                                                aaTGCCAAAAGCTTCCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGTTTGGCGGCGCAGGAGCCTCAGCTGGCGGTAGCAGAGTTCAGGTTAACCAACTGaa                                                                                                  >  6:113349/1‑149 (MQ=255)
                                                    cgaaaaGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACgg                                                                                               <  1:97195/146‑1 (MQ=255)
                                                    ccAATAGCTTGCGTGGAGGTCCGTTGGTGATGAACCCGCTTAAGGGCGTTGACAGGCAGGTGTATGCGCTGGCGCAGGGCAATAATCTGGTTGGCGGCGCAGGAGCCTCCGCTGTCGGTAGCAGTGTTCAGGTTAACCAACTGAACgg                                                                                               <  6:266135/148‑1 (MQ=255)
                                                     cAAAAGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTg                                                                                             <  7:96799/149‑1 (MQ=255)
                                                      aaaaGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTg                                                                                             >  8:216957/1‑148 (MQ=255)
                                                      aaaaGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTg                                                                                             >  3:188824/1‑148 (MQ=255)
                                                         aGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACg                                                                                         <  5:143865/149‑1 (MQ=255)
                                                           cTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCTTCCGCTGGCg                                                                                                                              <  2:253150/110‑1 (MQ=255)
                                                           cTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGa                                                                                       >  3:312953/1‑149 (MQ=255)
                                                           cTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCATCCGCTGGCg                                                                                                                              >  1:253150/1‑110 (MQ=255)
                                                               cGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCAc                                                                                   <  4:188824/149‑1 (MQ=255)
                                                                 tGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCAcc                                                                                  >  5:342594/1‑148 (MQ=255)
                                                                   gAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAAt                                                                               <  6:229596/149‑1 (MQ=255)
                                                                   gAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAAt                                                                               <  5:184865/149‑1 (MQ=255)
                                                                    aGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATg                                                                              >  1:145123/1‑149 (MQ=255)
                                                                      gTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGt                                                                            >  8:225243/1‑149 (MQ=255)
                                                                       tACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCACTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTg                                                                           <  5:273091/149‑1 (MQ=255)
                                                                        aCGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGc                                                                          >  5:241499/1‑149 (MQ=255)
                                                                              ttgATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTa                                                                     <  7:352615/148‑1 (MQ=255)
                                                                                gatgaCACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTAtt                                                                   >  7:103562/1‑148 (MQ=255)
                                                                                gatgaCACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTAtt                                                                   >  7:37745/1‑148 (MQ=255)
                                                                                gatgaCACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTAtt                                                                   <  6:158134/148‑1 (MQ=255)
                                                                                 atgaCACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTg                                                                  >  3:310595/1‑148 (MQ=255)
                                                                                  tgaCACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGaa                                                                <  5:195200/149‑1 (MQ=255)
                                                                                   gaCACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGaa                                                                <  4:171902/148‑1 (MQ=255)
                                                                                    acacCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAAc                                                               <  3:281935/148‑1 (MQ=255)
                                                                                      acCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGt                                                             >  1:264940/1‑148 (MQ=255)
                                                                                      acCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTg                                                            <  7:216957/149‑1 (MQ=255)
                                                                                      acCGCTCAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCAGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGCTCACCAATGGTGCGGTTATTGAACGt                                                             >  2:174308/1‑148 (MQ=255)
                                                                                      acCGCTCAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCAGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGt                                                             >  2:174306/1‑148 (MQ=255)
                                                                                         gCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGaag                                                         <  5:263712/149‑2 (MQ=255)
                                                                                         gCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAAt                                                         <  4:312953/149‑1 (MQ=255)
                                                                                         gCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAAt                                                         <  1:168450/149‑1 (MQ=255)
                                                                                          cTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAAtt                                                        <  7:38579/149‑1 (MQ=255)
                                                                                             aaGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGCCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAATTGcc                                                     >  2:42508/1‑149 (MQ=255)
                                                                                              aGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAATTGccc                                                    <  7:225243/149‑1 (MQ=255)
                                                                                               gggTGTTGACAGTCTGGTGTATGCGCTGGCGCAGGGGAATTTTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTTACCGGTGGACGGATCTCCAATGGTGCGGTTTTTGAACGTGAATTGCcca                                                   <  8:37745/149‑1 (MQ=255)
                                                                                                     tGACAGTCAGGTGTATGTGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAATTGCCCAGCCAGt                                             >  4:229752/1‑149 (MQ=255)
                                                                                                       aCAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAATTGCCCAGCCAGttt                                           <  5:252450/149‑1 (MQ=255)
                                                                                                        cAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAATTGCCCAGCCAGTTTg                                          <  4:310595/149‑1 (MQ=255)
                                                                                                         aGTCAGGTGTATGCGCTGGCGCAGGGCAATTTTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGATCGTGAATTGCCCAGCCAGTTTgg                                         >  8:81168/1‑149 (MQ=255)
                                                                                                           tCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCCACTGAACGGTGGACGGATCACCAATGGTGCGGGTAGTGAACGTGAATTGCCGAGCCAGTTAGGCg                                       >  8:224167/1‑149 (MQ=255)
                                                                                                               gtgtATGCGCTGGCGCAGGGCAATATTGTGGTTGGCGGCGCAGGAGCCTGCGCTGGCGGTAGGAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGGGCGGTTATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCgg                                   <  5:113349/149‑1 (MQ=255)
                                                                                                               gtgtATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCgg                                   <  3:102512/149‑1 (MQ=255)
                                                                                                               gtgtATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCgg                                   <  2:264940/149‑1 (MQ=255)
                                                                                                               gtgtATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCgg                                   >  2:257672/1‑149 (MQ=255)
                                                                                                                 gtATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCggg                                  >  1:233375/1‑148 (MQ=255)
                                                                                                                    tGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAAt                               <  5:333220/148‑1 (MQ=255)
                                                                                                                     gcgcTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATAc                             >  8:72884/1‑149 (MQ=255)
                                                                                                                       gcTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGCCGGATCACAAATGGTGCGGTTATTGAAAGTGAATTGCCCAGCCAGTTTGGCGTAGGGAATAacc                           >  6:221703/1‑149 (MQ=255)
                                                                                                                       gcTGGCGCAGGGCAATATTCTGGTTGGCAGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATAccc                           <  1:174306/149‑1 (MQ=255)
                                                                                                                       gcTGGCGCAGGGCAATAATCTGGTTGGCAGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTGAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATAccc                           <  1:174308/149‑1 (MQ=255)
                                                                                                                                      tatTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACtt             >  3:156032/1‑148 (MQ=255)
                                                                                                                                       atTCTGGTTGGCGGCGCAGGAGCCTCCTCTGGCGGTAGCAGTGTTCAGCTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTaa           <  1:114159/149‑1 (MQ=255)
                                                                                                                                        ttCTGGTTGGCGGCGCAGGAGCCTCCGCTGTCTGTAGCAGTGTTCAGGATAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTaa           <  8:103562/148‑1 (MQ=255)
                                                                                                                                             gTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAAcgacga     <  1:42508/149‑1 (MQ=255)
                                                                                                                                                ggcggcGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAg   <  2:145123/148‑1 (MQ=255)
                                                                                                                                                 gcggcgCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGa  <  1:137129/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TTTGGACGTCAGGGGCAAACCATCGATGTGGTGGTTTCTTCCATGGGAAATGCCAAAAGCTTGCGTGGAGGTACGTTGTTGATGACACCGCTTAAGGGCGTTGACAGTCAGGTGTATGCGCTGGCGCAGGGCAATATTCTGGTTGGCGGCGCAGGAGCCTCCGCTGGCGGTAGCAGTGTTCAGGTTAACCAACTGAACGGTGGACGGATCACCAATGGTGCGGTTATTGAACGTGAATTGCCCAGCCAGTTTGGCGTCGGGAATACCCTTAATTTGCAACTTAACGACGAAGA  >  NC_000913/1136559‑1136851

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: