Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,136,026 G→T 5.9% V155L (GTA→TTA) ‡ flgH → flagellar L‑ring protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,136,0260GT 5.9% 286.9 / 5.1 85V155L (GTA→TTA) ‡flgHflagellar L‑ring protein
Reads supporting (aligned to +/- strand):  ref base G (43/37);  new base T (3/2);  total (46/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GCTTTGATACTGTGCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCG  >  NC_000913/1135877‑1136172
                                                                                                                                                     |                                                                                                                                                  
gCTTTGATACTGTGCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAgg                                                                                                                                                    <  2:1487623/150‑1 (MQ=255)
 cTTTGATACTGTGCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGt                                                                                                                                                   >  2:852776/1‑150 (MQ=255)
 cTTTGATACTGTGCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGt                                                                                                                                                   >  2:634245/1‑150 (MQ=255)
   ttGATACTGTGCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTAc                                                                                                                                                 <  2:480064/150‑1 (MQ=255)
    tGATACTGTGCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACt                                                                                                                                                >  1:965044/1‑150 (MQ=255)
      aTACTGTGCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTgg                                                                                                                                              >  1:621926/1‑150 (MQ=255)
       tACTGTGCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGt                                                                                                                                             >  2:458626/1‑150 (MQ=255)
         cTGTGCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCa                                                                                                                                           <  2:312379/150‑1 (MQ=255)
           gtgCCTCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAACGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGTACCTGGTCAAc                                                                                                                                         <  1:1512560/150‑1 (MQ=255)
            tgCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACg                                                                                                                                        <  2:1550685/150‑1 (MQ=255)
             gCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACgg                                                                                                                                       >  2:108130/1‑150 (MQ=255)
              ccGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGc                                                                                                                                      <  2:191928/150‑1 (MQ=255)
                gcgcTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCaa                                                                                                                                    >  2:788390/1‑150 (MQ=255)
                    tATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTg                                                                                                                                <  2:1363908/150‑1 (MQ=255)
                        tGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCAtg                                                                                                                            <  2:1257072/150‑1 (MQ=255)
                         gCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGTACCTGGTCAAccgtc                                                                                                                                     <  2:356839/140‑5 (MQ=255)
                         gCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGTACCTGGTCAAccgtc                                                                                                                                     >  1:356839/1‑136 (MQ=255)
                         gCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACAGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCAtgt                                                                                                                           >  1:409163/1‑150 (MQ=255)
                          cAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCAtgtg                                                                                                                          >  2:256855/1‑150 (MQ=255)
                             gggCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATgtggtg                                                                                                                       <  1:1525443/150‑1 (MQ=255)
                               gCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTggtggg                                                                                                                     <  2:777233/150‑1 (MQ=255)
                               gCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTggtggg                                                                                                                     <  1:489374/150‑1 (MQ=255)
                                cTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTggtgggt                                                                                                                    >  1:1028543/1‑150 (MQ=255)
                                     tGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGaaaa                                                                                                               <  2:1630194/150‑1 (MQ=255)
                                     tGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGaaaa                                                                                                               <  1:77872/150‑1 (MQ=255)
                                      ggTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCTACCTGCATGTGGTGGGTGaaaaa                                                                                                              >  1:722990/1‑150 (MQ=255)
                                      ggTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGaaaaa                                                                                                              >  2:718841/1‑150 (MQ=255)
                                        tAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGc                                                                                                                                      <  2:1348467/124‑1 (MQ=255)
                                        tAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGc                                                                                                                                      >  1:1348467/1‑124 (MQ=255)
                                         aaCGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAg                                                                                                           <  1:543432/150‑1 (MQ=255)
                                              tCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGc                                                                                                      <  1:185020/150‑1 (MQ=255)
                                                 tGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCAt                                                                                                   >  1:1373690/1‑150 (MQ=255)
                                                 tGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCAt                                                                                                   <  1:634245/150‑1 (MQ=255)
                                                   ccGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTa                                                                                                 >  1:1387105/1‑150 (MQ=255)
                                                    cGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTaa                                                                                                >  2:359808/1‑150 (MQ=255)
                                                      aTGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATc                                                                                              <  1:788390/150‑1 (MQ=255)
                                                          cGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAggg                                                                                          <  2:1512352/150‑1 (MQ=255)
                                                          cGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAggg                                                                                          >  2:452042/1‑150 (MQ=255)
                                                             aGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGCGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTAc                                                                                       >  1:488574/1‑150 (MQ=255)
                                                              gCCTCCGGTGGTAACACGTTCAACAGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTAcc                                                                                      <  2:409163/150‑1 (MQ=255)
                                                                    ggtggtAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAAttt                                                                                >  2:779603/1‑150 (MQ=255)
                                                                         tAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCg                                                                           >  1:753300/1‑150 (MQ=255)
                                                                         tAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCg                                                                           >  2:1335029/1‑150 (MQ=255)
                                                                          aaCACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGc                                                                          >  2:878969/1‑150 (MQ=255)
                                                                           acacGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCt                                                                         <  1:316656/150‑1 (MQ=255)
                                                                            cacGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCtt                                                                        >  1:796088/1‑150 (MQ=255)
                                                                               gTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTtctc                                                                     <  2:962907/150‑1 (MQ=255)
                                                                                  cAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTgg                                                                  >  1:980112/1‑150 (MQ=255)
                                                                                  cAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTgg                                                                  >  2:1466057/1‑150 (MQ=255)
                                                                                     cGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGt                                                               <  2:1577579/150‑1 (MQ=255)
                                                                                      ggAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTg                                                              >  1:394129/1‑150 (MQ=255)
                                                                                              cGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATcc                                                      <  1:779603/150‑1 (MQ=255)
                                                                                                  gCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGGTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGc                                                  <  2:1462464/150‑1 (MQ=255)
                                                                                                        gCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATc                                            >  1:1196733/1‑150 (MQ=255)
                                                                                                         ccAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCa                                           >  1:64031/1‑150 (MQ=255)
                                                                                                         ccAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCa                                           >  2:333137/1‑150 (MQ=255)
                                                                                                          cAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAg                                          >  2:70814/1‑150 (MQ=255)
                                                                                                           aGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGc                                         >  1:755085/1‑150 (MQ=255)
                                                                                                            gCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCg                                        <  2:552854/150‑1 (MQ=255)
                                                                                                            gCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCg                                        <  2:717330/150‑1 (MQ=255)
                                                                                                              aaTACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGgc                                      <  1:245033/150‑1 (MQ=255)
                                                                                                               aTACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGgca                                     <  1:544715/150‑1 (MQ=255)
                                                                                                                  ccTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCTATTAATTAGTGTAACGAATTTATTCGATTCTATGTCGTTTTTAATCCAAGAACTATCAGCTgcagca                                  <  2:1205532/150‑1 (MQ=255)
                                                                                                                  ccTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGgcagca                                  >  1:376890/1‑150 (MQ=255)
                                                                                                                   cTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCaa                                 >  1:131373/1‑150 (MQ=255)
                                                                                                                      tAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAAtac                              <  2:1373690/150‑1 (MQ=255)
                                                                                                                          ggCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAAtaccgta                          >  1:1476582/1‑150 (MQ=255)
                                                                                                                           gCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAAtaccgtac                         <  1:836478/150‑1 (MQ=255)
                                                                                                                           gCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAAtaccgtac                         >  1:1516309/1‑150 (MQ=255)
                                                                                                                              cGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAAtaccg                            >  2:857007/1‑144 (MQ=255)
                                                                                                                              cGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAAtacag                            <  1:857007/144‑3 (MQ=255)
                                                                                                                                ttGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATc                                            >  1:843439/1‑126 (MQ=255)
                                                                                                                                ttGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATc                                            <  2:843439/126‑1 (MQ=255)
                                                                                                                                  gacggtgacggtTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTAc                  >  2:633987/1‑150 (MQ=255)
                                                                                                                                    cggtgacggtTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTc                <  1:718841/150‑1 (MQ=255)
                                                                                                                                      gtgacggtTGACCAGGTACTGGTCAACGGCTACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAg              <  2:722990/150‑1 (MQ=255)
                                                                                                                                       tgacggtTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAgg             >  2:382390/1‑150 (MQ=255)
                                                                                                                                         acggtTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTg           <  2:131373/150‑1 (MQ=255)
                                                                                                                                          cggtTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTgg          <  2:980112/150‑1 (MQ=255)
                                                                                                                                            gtTGACCAGTACCTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCg        >  1:1450937/1‑150 (MQ=255)
                                                                                                                                            gtTGACCAGTACCTGGTCAACGGCAACCGGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGTTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCAACGCACTATCAGCGGCAGCAATACCCTACAGTCTACTCAAGTGGCg        >  2:653402/1‑150 (MQ=255)
                                                                                                                                             ttGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCgg       >  1:429554/1‑150 (MQ=255)
                                                                                                                                              tgcccAGGTACTGGTCAACGGCAACCTGCATGTGGAGGGTGAAAAACACATTGCCATTTAACAGGCAACCGAATTTATTCGCTACTCTCGCGTGCTAAATCCACGCACTAACAGCGGCAGCAATACCGTACCCTCTACTCAGGTGGCCGa      >  2:1480752/4‑150 (MQ=255)
                                                                                                                                                aCCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATg    >  2:1122486/1‑150 (MQ=255)
                                                                                                                                                  cAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATgcg  >  1:69639/1‑150 (MQ=255)
                                                                                                                                                  cAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGCGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATgcg  <  1:359808/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCTTTGATACTGTGCCGCGCTATTTGCAGGGGCTGTTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACGGCAACCTGCATGTGGTGGGTGAAAAACAGATTGCCATTAATCAGGGTACCGAATTTATTCGCTTCTCTGGCGTGGTTAATCCACGCACTATCAGCGGCAGCAATACCGTACCGTCTACTCAGGTGGCGGATGCG  >  NC_000913/1135877‑1136172

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