Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,966,023 A→T 6.3% S57T (TCG→ACG)  flhB ← flagellar biosynthesis protein FlhB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,966,0230AT 6.3% 290.9 / 7.7 96S57T (TCG→ACG) flhBflagellar biosynthesis protein FlhB
Reads supporting (aligned to +/- strand):  ref base A (44/46);  new base T (3/3);  total (47/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.80e-01

GAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTG  >  NC_000913/1965875‑1966171
                                                                                                                                                    |                                                                                                                                                    
gaAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGAc                                                                                                                                                     <  1:787151/150‑1 (MQ=255)
    tGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCCTGCCCGACAATc                                                                                                                                                 >  2:1212038/1‑150 (MQ=255)
          ccaccaGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACggg                                                                                                                                           <  1:1329994/150‑1 (MQ=255)
            accaGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATCAAGCCCAGCGGAGAGCATGCCCGACAATCTACGGGcc                                                                                                                                         >  1:1378566/1‑150 (MQ=255)
              caGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGTCGATTGTCGGGcatg                                                                                                                                       <  2:529402/150‑4 (MQ=255)
               aGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGTCGATTGTCGGGcatgc                                                                                                                                      <  1:351339/150‑5 (MQ=255)
               aGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGc                                                                                                                                      >  1:172710/1‑150 (MQ=255)
                cgctgctgccgcTGATTAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCg                                                                                                                                     >  2:281101/8‑150 (MQ=255)
                   ccacGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGaca                                                                                                                                  <  1:73546/150‑1 (MQ=255)
                    cacGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGacac                                                                                                                                 <  1:678814/150‑1 (MQ=255)
                    cacGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGacac                                                                                                                                 >  1:1428266/1‑150 (MQ=255)
                      cgccgcTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGacacac                                                                                                                               >  1:1107528/1‑150 (MQ=255)
                       gccgcTGATTAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGAGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACAcc                                                                                                                              >  1:1073092/1‑150 (MQ=255)
                           cTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCa                                                                                                                          >  1:499100/1‑150 (MQ=255)
                              aTCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAAc                                                                                                                       <  2:1179443/150‑1 (MQ=255)
                              aTCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAAc                                                                                                                       >  1:1523299/1‑150 (MQ=255)
                                 aGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAg                                                                                                                    <  2:475134/150‑1 (MQ=255)
                                 aGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCAAGCGACACACCGCCAAACCAg                                                                                                                    >  2:549709/1‑150 (MQ=255)
                                  gCGGCAGCAGCGCCAGCATGGTTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGa                                                                                                                   <  1:768379/150‑1 (MQ=255)
                                      cagcagCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAac                                                                                                               <  2:1420564/150‑1 (MQ=255)
                                       agcagcGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAaca                                                                                                              >  2:1210124/1‑150 (MQ=255)
                                        gcagcGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAacac                                                                                                             >  1:1407854/1‑150 (MQ=255)
                                        gcagcGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAacac                                                                                                             <  1:1033020/150‑1 (MQ=255)
                                        gcagcGCAAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAacac                                                                                                             <  1:188543/150‑1 (MQ=255)
                                            cgcCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAAc                                                                                                         >  2:44790/1‑150 (MQ=255)
                                             gcCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGTAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACg                                                                                                        <  1:5566/150‑1 (MQ=255)
                                             gcCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACg                                                                                                        >  1:216527/1‑150 (MQ=255)
                                             gcCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACg                                                                                                        >  1:728258/1‑150 (MQ=255)
                                               cAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGcc                                                                                                      <  1:726740/150‑1 (MQ=255)
                                                aGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAACTGCAGCCCAGCGGAGAGCAAGCCCGACAATCGCCGGGCCAGCGACACCCCGCCAAACCAGATAACACTAACGccc                                                                                                     >  1:1512352/1‑150 (MQ=255)
                                                  cATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCccac                                                                                                   <  2:499100/150‑1 (MQ=255)
                                                     ggCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAg                                                                                                <  1:263440/150‑1 (MQ=255)
                                                      gCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGTCGATTGTCGGGcat                                                                                                                                        <  2:25336/109‑3 (MQ=255)
                                                      gCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGTCGATTGTCGGGcat                                                                                                                                        >  1:25336/1‑107 (MQ=255)
                                                       cTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCa                                                                                              >  1:362318/1‑150 (MQ=255)
                                                          ctctGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCAAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCaaaa                                                                                           <  1:549709/150‑1 (MQ=255)
                                                            ctGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATc                                                                                         <  1:1303143/150‑1 (MQ=255)
                                                              gATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATcac                                                                                       >  2:1152995/1‑149 (MQ=255)
                                                              gATCACCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATcag                                                                                       <  1:756953/150‑1 (MQ=255)
                                                                tCAGCAGAATAATCTGCCCGAGGATCAGATTCGGTTCATTGATAAAACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGCGCCAGAGACACACCGCCAAACCAGATAAAAAAAACCCCCACCAGCAAAATcagca                                                                                     >  1:1512596/1‑150 (MQ=255)
                                                                tCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATcagca                                                                                     <  2:1523299/150‑1 (MQ=255)
                                                                 cagcagAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGACAACACTAACGCCCACCAGCAAAATcagcag                                                                                    <  2:1396672/150‑1 (MQ=255)
                                                                   gcagAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACTGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTg                                                                                  <  2:1482916/150‑1 (MQ=255)
                                                                    cagAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGa                                                                                 >  2:1037326/1‑150 (MQ=255)
                                                                     agAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAg                                                                                >  1:322701/1‑150 (MQ=255)
                                                                           agcTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGt                                                                          <  2:949804/148‑1 (MQ=255)
                                                                             cTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTc                                                                        >  2:1097340/1‑150 (MQ=255)
                                                                              tGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCa                                                                       <  1:815983/150‑1 (MQ=255)
                                                                               gCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCAc                                                                      <  2:125000/150‑1 (MQ=255)
                                                                               gCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCAc                                                                      >  1:837867/1‑150 (MQ=255)
                                                                               gCCCGAGGATCAGATTCGGGTCATTGATAATACCGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCCGCGCCACACCGCCAACCCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCAc                                                                      >  1:1425471/1‑150 (MQ=255)
                                                                                 ccGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACgg                                                                    <  1:281101/150‑1 (MQ=255)
                                                                                  cGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACggg                                                                   <  2:395169/150‑1 (MQ=255)
                                                                                  cGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACggg                                                                   <  2:810560/150‑1 (MQ=255)
                                                                                  cGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACggg                                                                   <  2:925110/150‑1 (MQ=255)
                                                                                   gAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGa                                                                  >  1:341705/1‑150 (MQ=255)
                                                                                       aTCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAAcgc                                                              >  2:1155633/1‑150 (MQ=255)
                                                                                        tCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAAcgcg                                                             >  1:1523667/1‑150 (MQ=255)
                                                                                          aGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCggg                                                           <  2:216527/150‑1 (MQ=255)
                                                                                                gggTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGc                                                     >  2:851835/1‑150 (MQ=255)
                                                                                                 ggTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGcc                                                    >  1:71729/1‑150 (MQ=255)
                                                                                                 ggTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGcc                                                    <  2:71729/150‑1 (MQ=255)
                                                                                                 ggTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGcc                                                    >  2:375670/1‑150 (MQ=255)
                                                                                                  gTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGccc                                                   >  1:1357807/1‑150 (MQ=255)
                                                                                                       tGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCtt                                              >  1:77761/1‑150 (MQ=255)
                                                                                                          taataCTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTccc                                           <  1:851835/150‑1 (MQ=255)
                                                                                                          taataCTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAATGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTccc                                           >  2:222766/1‑150 (MQ=255)
                                                                                                            ataCTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCcgc                                         >  1:907479/1‑150 (MQ=255)
                                                                                                             taCTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCcgcg                                        <  2:172710/150‑1 (MQ=255)
                                                                                                                  tgATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGTAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCtttt                                   <  2:362318/150‑1 (MQ=255)
                                                                                                                  tgATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCtttt                                   >  1:489148/1‑150 (MQ=255)
                                                                                                                  tgATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCtttt                                   >  1:509478/1‑150 (MQ=255)
                                                                                                                   gATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCttttt                                  >  1:1480815/1‑150 (MQ=255)
                                                                                                                     tCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCt                                <  1:299780/150‑1 (MQ=255)
                                                                                                                     tCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCt                                <  2:1073092/150‑1 (MQ=255)
                                                                                                                       aaaaTGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAg                              >  2:311771/1‑150 (MQ=255)
                                                                                                                        aaaTGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGt                             <  2:1428266/150‑1 (MQ=255)
                                                                                                                        aaaTGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCAGCGCTTTTTCTAGt                             >  2:375255/1‑150 (MQ=255)
                                                                                                                                 ccAGCGGAGAGCATGCCCGACAATCGACGTGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGgt                    >  2:780984/1‑150 (MQ=255)
                                                                                                                                  cAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGgtg                   <  2:523643/150‑1 (MQ=255)
                                                                                                                                  cAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGgtg                   >  1:628758/1‑150 (MQ=255)
                                                                                                                                     cGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTgg                <  1:1155633/150‑1 (MQ=255)
                                                                                                                                     cGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTgg                <  2:509478/150‑1 (MQ=255)
                                                                                                                                        gtcgctggCCCGTCGATTGTCGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTggggg             >  1:1079851/8‑150 (MQ=255)
                                                                                                                                         gagCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAATGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGc            <  1:222766/150‑1 (MQ=255)
                                                                                                                                          cgctggCCCGTCGATTGTCGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCCCCAGCAAAAGCAGCAGGGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGGGTGGGGGCt           >  1:1315793/6‑150 (MQ=255)
                                                                                                                                           gCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCtt          >  1:18011/1‑150 (MQ=255)
                                                                                                                                             aTGCCCGACAATCGACGGGCCAGCGACACACCTCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGTGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCt        <  2:579682/150‑1 (MQ=255)
                                                                                                                                               gCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGt      >  1:660589/1‑150 (MQ=255)
                                                                                                                                                 ccGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGttt    <  1:1097340/150‑1 (MQ=255)
                                                                                                                                                 ccGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGttt    <  2:1407854/150‑1 (MQ=255)
                                                                                                                                                 ccGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGttt    >  2:228423/1‑150 (MQ=255)
                                                                                                                                                  cGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGtttt   <  2:837867/150‑1 (MQ=255)
                                                                                                                                                   gACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTg  <  1:44790/150‑1 (MQ=255)
                                                                                                                                                   gACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTg  <  2:341705/150‑1 (MQ=255)
                                                                                                                                                   gACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGGCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTg  <  2:592936/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTG  >  NC_000913/1965875‑1966171

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: