Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,966,030 A→T 6.3% R54R (CGT→CGA flhB ← flagellar biosynthesis protein FlhB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,966,0300AT 6.3% 281.9 / 7.7 95R54R (CGT→CGAflhBflagellar biosynthesis protein FlhB
Reads supporting (aligned to +/- strand):  ref base A (42/46);  new base T (3/3);  total (46/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 = 3.75e-01

CCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCT  >  NC_000913/1965885‑1966179
                                                                                                                                                 |                                                                                                                                                     
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                                                                                                                               >  1:1523299/1‑150 (MQ=255)
                    aTCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAAc                                                                                                                               <  2:1179443/150‑1 (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:728258/1‑150 (MQ=255)
                                   gcCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACg                                                                                                                >  1:216527/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                                                                                                                                                >  1:25336/1‑107 (MQ=255)
                                            gCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGTCGATTGTCGGGcat                                                                                                                                                <  2:25336/109‑3 (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                                                                              >  1:837867/1‑150 (MQ=255)
                                                                     gCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCAc                                                                              <  2:125000/150‑1 (MQ=255)
                                                                     gCCCGAGGATCAGATTCGGGTCATTGATAATACCGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCCGCGCCACACCGCCAACCCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCAc                                                                              >  1:1425471/1‑150 (MQ=255)
                                                                       ccGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACgg                                                                            <  1:281101/150‑1 (MQ=255)
                                                                        cGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACggg                                                                           <  2:810560/150‑1 (MQ=255)
                                                                        cGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACggg                                                                           <  2:395169/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                                                            <  2:71729/150‑1 (MQ=255)
                                                                                       ggTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGcc                                                            >  2:375670/1‑150 (MQ=255)
                                                                                       ggTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGcc                                                            >  1:71729/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                        <  2:509478/150‑1 (MQ=255)
                                                                                                                           cGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTgg                        <  1:1155633/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            >  2:228423/1‑150 (MQ=255)
                                                                                                                                       ccGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGttt            <  1:1097340/150‑1 (MQ=255)
                                                                                                                                       ccGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGttt            <  2:1407854/150‑1 (MQ=255)
                                                                                                                                        cGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGtttt           <  2:837867/150‑1 (MQ=255)
                                                                                                                                         gACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTg          <  2:341705/150‑1 (MQ=255)
                                                                                                                                         gACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTg          <  1:44790/150‑1 (MQ=255)
                                                                                                                                         gACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGGCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTg          <  2:592936/150‑1 (MQ=255)
                                                                                                                                              tCGACGTGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTgtcgtc     <  1:780984/150‑1 (MQ=255)
                                                                                                                                                 aCGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGct  >  2:1362241/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCT  >  NC_000913/1965885‑1966179

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: