Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,158,513 C→A 8.2% N215K (AAC→AAA ptsG → glucose‑specific PTS enzyme IIBC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,158,5130CA 8.2% 280.0 / 10.5 85N215K (AAC→AAAptsGglucose‑specific PTS enzyme IIBC component
Reads supporting (aligned to +/- strand):  ref base C (37/41);  new base A (4/3);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

TGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCC  >  NC_000913/1158368‑1158660
                                                                                                                                                 |                                                                                                                                                   
tgCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTAc                                                                                                                                                 <  1:894012/150‑1 (MQ=255)
 gCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTAcc                                                                                                                                                >  1:1296718/1‑150 (MQ=255)
   tGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCtt                                                                                                                                              <  1:613258/150‑1 (MQ=255)
    gTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCttt                                                                                                                                             >  1:1533461/1‑150 (MQ=255)
      ccTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTcc                                                                                                                                           <  1:905022/150‑1 (MQ=255)
      ccTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTcc                                                                                                                                           <  2:238710/150‑1 (MQ=255)
       cTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCa                                                                                                                                          >  1:163729/1‑150 (MQ=255)
       cTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCa                                                                                                                                          >  2:1092567/1‑150 (MQ=255)
           aTTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATGGATCGTTGCCTGGTACCGTTTGGTCTGCACGACATCTGGAAAGGTACGTTCCAGATg                                                                                                                                      <  2:956888/150‑1 (MQ=255)
           aTTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATg                                                                                                                                      >  2:121068/1‑150 (MQ=255)
           aTTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAAAGGTACGTTCCAGATg                                                                                                                                      <  1:732824/150‑1 (MQ=255)
                 ccgccgATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGAtt                                                                                                                                >  2:108545/1‑150 (MQ=255)
                 ccgccgATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAAAGGTACGTTCCAGATgtggtgc                                                                                                                                <  1:476624/150‑7 (MQ=255)
                   gccgATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGATGCTTACGAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTgg                                                                                                                              <  2:719379/150‑1 (MQ=255)
                   gccgATTGGTTCTGCAATCCAGACCTTCTCTCAGTCGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTTCCTGGTACCGTTTGCTCTGCACCACATCTGGAACGTCCCTATCCAGATGCAGATTgg                                                                                                                              >  2:1025394/1‑150 (MQ=255)
                    ccgATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGt                                                                                                                             >  2:761809/1‑150 (MQ=255)
                     cgATTGGTTCTGCAATCCAGACCTTCTCTCAGTTGGCTGCTTACCAGAACCCGGTATTTGCTTTTGGCTTTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTTGTCTTCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTg                                                                                                                            <  2:658583/150‑1 (MQ=255)
                       aTTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGaa                                                                                                                          <  1:108545/150‑1 (MQ=255)
                          ggTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATac                                                                                                                       >  2:462067/1‑150 (MQ=255)
                           gTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATaca                                                                                                                      >  2:1501051/1‑150 (MQ=255)
                           gTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATaca                                                                                                                      >  1:567883/1‑150 (MQ=255)
                              cTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCa                                                                                                                   >  2:991607/1‑150 (MQ=255)
                              cTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCa                                                                                                                   <  2:700582/150‑1 (MQ=255)
                               tGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCaa                                                                                                                  >  1:30087/1‑150 (MQ=255)
                               tGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCaa                                                                                                                  >  2:62071/1‑150 (MQ=255)
                               tGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCaa                                                                                                                  >  1:597414/1‑150 (MQ=255)
                                gCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAAc                                                                                                                 <  2:1627617/150‑1 (MQ=255)
                                gCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAAc                                                                                                                 <  2:1555427/150‑1 (MQ=255)
                                  aaTCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACgc                                                                                                               >  1:880218/1‑150 (MQ=255)
                                   aTCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACgca                                                                                                              >  1:1541894/1‑150 (MQ=255)
                                     ccAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTTGCATTTACTGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGTATTCACCaatctgca                                                                                                            <  1:385985/150‑7 (MQ=255)
                                     ccAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACgcagc                                                                                                            <  1:238491/150‑1 (MQ=255)
                                      cAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACgcagca                                                                                                           <  2:350274/150‑1 (MQ=255)
                                      cAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACgcagca                                                                                                           <  2:528252/150‑1 (MQ=255)
                                        gACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACTGTTTTGTTTTCAGATAATCTGTAACGTACCTTTTAATATGCAGATTTGTGAATACATAAACGCAGctgg                                                                                                         <  2:886402/150‑1 (MQ=255)
                                         aCCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGcaggt                                                                                                        <  1:1501051/150‑1 (MQ=255)
                                            ttCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGcaggtcag                                                                                                     <  1:890493/150‑1 (MQ=255)
                                              ctctcAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGTTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGcaggtcaggt                                                                                                   <  1:1025394/150‑1 (MQ=255)
                                                 tcAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGtttt                                                                                                <  2:1141530/150‑1 (MQ=255)
                                                  cAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTc                                                                                               <  2:645594/150‑1 (MQ=255)
                                                  cAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTc                                                                                               >  2:446802/1‑150 (MQ=255)
                                                   aGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTcc                                                                                              <  1:325079/150‑1 (MQ=255)
                                                        gctgctTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGc                                                                                         >  2:1576686/1‑150 (MQ=255)
                                                         ctgctTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCg                                                                                        <  1:1117444/150‑1 (MQ=255)
                                                         ctgctTACCAGAACCCGGTAGTTGCGTTCGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCg                                                                                        <  2:1504/150‑1 (MQ=255)
                                                           gctTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGAc                                                                                      >  2:1629782/1‑150 (MQ=255)
                                                           gctTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGAc                                                                                      <  1:1411307/150‑1 (MQ=255)
                                                                ccAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTcc                                                                                 <  2:1533461/150‑1 (MQ=255)
                                                                 cAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCcg                                                                                >  2:258019/1‑150 (MQ=255)
                                                                 cAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCcg                                                                                >  2:958533/1‑150 (MQ=255)
                                                                               tGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCTCGGCGACATTCCGCGTTATATGGCggg                                                                  >  1:608872/1‑150 (MQ=255)
                                                                               tGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCggg                                                                  <  2:248857/150‑1 (MQ=255)
                                                                                gCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGTCGGGt                                                                 >  1:1204325/1‑150 (MQ=255)
                                                                                  gTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGa                                                               <  1:1629782/150‑1 (MQ=255)
                                                                                     tGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGAccc                                                            <  1:258019/150‑1 (MQ=255)
                                                                                      ggCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCg                                                           >  2:389702/1‑150 (MQ=255)
                                                                                      ggCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCg                                                           <  1:471030/150‑1 (MQ=255)
                                                                                        cATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGAc                                                         <  1:259338/150‑1 (MQ=255)
                                                                                            tACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCg                                                     >  2:1341980/1‑150 (MQ=255)
                                                                                              cGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCggg                                                   >  2:238769/1‑150 (MQ=255)
                                                                                                gTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTa                                                 <  1:1571856/150‑1 (MQ=255)
                                                                                                  ttCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTaaa                                               >  2:1450844/1‑150 (MQ=255)
                                                                                                   tCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAAc                                              <  2:522416/150‑1 (MQ=255)
                                                                                                     aTCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTg                                            >  2:689329/1‑150 (MQ=255)
                                                                                                      tCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGt                                           >  1:1538362/1‑150 (MQ=255)
                                                                                                        gAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCt                                         >  2:749226/1‑150 (MQ=255)
                                                                                                          aCGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCtgg                                       <  1:62071/150‑1 (MQ=255)
                                                                                                            gTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCtggtg                                     <  1:462067/150‑1 (MQ=255)
                                                                                                               gCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCt                                  >  2:465946/1‑150 (MQ=255)
                                                                                                                 cTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTc                                <  2:597414/150‑1 (MQ=255)
                                                                                                                    gTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTg                             <  1:465946/150‑1 (MQ=255)
                                                                                                                       ccGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAATCTGTCTGGTGGCTTCCTGTTc                          >  1:1504107/1‑150 (MQ=255)
                                                                                                                           ttGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCaaaa                      <  1:389702/150‑1 (MQ=255)
                                                                                                                               tCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTa                  <  1:991607/150‑1 (MQ=255)
                                                                                                                               tCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTa                  <  2:1538362/150‑1 (MQ=255)
                                                                                                                                 caatctgcaTCTGGAAAGGTACGTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACg                >  1:185477/8‑150 (MQ=255)
                                                                                                                                    accacATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTc             >  2:1552919/1‑150 (MQ=255)
                                                                                                                                     ctgcaTCTGGAAAGGTACGTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCt            >  1:365895/4‑150 (MQ=255)
                                                                                                                                      tgcaTCTGGAAAGGTACGTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTg           >  1:250066/3‑150 (MQ=255)
                                                                                                                                      cacaTCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTg           <  2:163729/150‑1 (MQ=255)
                                                                                                                                      cacaTCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGAGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTg           >  2:712914/1‑150 (MQ=255)
                                                                                                                                       gcaTCTGGAAAGGTACGTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGc          >  2:1160893/2‑150 (MQ=255)
                                                                                                                                             ggAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATg                    >  1:338235/1‑134 (MQ=255)
                                                                                                                                             ggAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATg                    <  2:338235/134‑1 (MQ=255)
                                                                                                                                               aaCGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGcc  >  2:191113/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TGCTGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCC  >  NC_000913/1158368‑1158660

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: