Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,173,155 A→G 5.5% pseudogene (208/442 nt) gatC ← pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,1550AG 5.5% 385.3 / 9.2 128pseudogene (208/442 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base A (67/53);  new base G (5/2);  total (72/56)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.65e-01

Only 100 of 132 total aligned reads displayed.

CCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGG  >  NC_000913/2173006‑2173304
                                                                                                                                                     |                                                                                                                                                     
ccGAAAATAATGAAACCGTGAATAATTAGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCTGCTTATAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATa                                                                                                                                                       <  1:234282/150‑1 (MQ=255)
  gATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGt                                                                                                                                                     >  1:1238460/1‑150 (MQ=255)
  gATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGt                                                                                                                                                     >  2:1709592/1‑150 (MQ=255)
    taataaTGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCGGCCAGTTGTGTGTGCAAACCAAGAGgtt                                                                                                                                                   >  2:208123/1‑150 (MQ=255)
     aataatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTg                                                                                                                                                  <  2:1076550/150‑1 (MQ=255)
        aatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGTGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGGTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTCTGCAAACCAATATTTTGCGt                                                                                                                                               >  2:37729/1‑150 (MQ=255)
        aatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACTATTGGTTTGcac                                                                                                                                               <  1:766767/150‑3 (MQ=255)
        aatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGt                                                                                                                                               >  1:1307046/1‑150 (MQ=255)
          tGAAACCGGGAAAATTTTGCGCTGAGAATACCTGAATAAGTAACCATGTTATGTGAGTACCGCCCTGATCCATTGAAGCCACCATACCCGCGGCTTTTAACGCCCCAGCATTATCCGCCAGTTGGGTGTGCAAACTATTGGTTTGcacac                                                                                                                                             <  1:1096105/150‑5 (MQ=255)
           gAAACCGTGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTTATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTcgc                                                                                                                                            <  2:283408/150‑1 (MQ=255)
           gAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCTAACCAATAGTTTGCGTcgc                                                                                                                                            >  2:171080/1‑150 (MQ=255)
           gAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTcgc                                                                                                                                            <  1:200723/150‑1 (MQ=255)
           gAAACCGGGAATATTTTGCGGGGAGAAAACCTGAAACAGTAACCAGGTAATGGGAGAACCTCCCTGATCAATTGAAGCCACCATACACCCGGCTTTTAACGCCCCAGCAATAGCCGCCAGTTGGCTGAGCAACCAAATAGTTTGAGTcgg                                                                                                                                            >  2:1344924/1‑149 (MQ=255)
              aCCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAAACGCCCTGATCCATTGACGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGGGTGCAAACCAATAGGTTGCGTCGCGAt                                                                                                                                         <  2:937831/150‑1 (MQ=255)
                cGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATcc                                                                                                                                       <  1:1638622/150‑1 (MQ=255)
                  ggAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCcac                                                                                                                                     >  2:1500229/1‑150 (MQ=255)
                       atTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGt                                                                                                                                <  1:1126358/150‑1 (MQ=255)
                         tttGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGa                                                                                                                              <  1:1466995/150‑1 (MQ=255)
                          ttGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGac                                                                                                                             >  2:78492/1‑149 (MQ=255)
                                ggAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCGATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCAt                                                                                                                       >  1:992969/1‑150 (MQ=255)
                                   gaAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAat                                                                                                                    <  1:825548/150‑1 (MQ=255)
                                   gaAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAat                                                                                                                    >  2:219167/1‑150 (MQ=255)
                                    aaatCCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATTCCCCCGGCTTTTAACGCCCCATCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGGGATGCTCATAatg                                                                                                                   <  2:1810228/150‑1 (MQ=255)
                                    aaaaCCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATaacg                                                                                                                   <  2:1655227/150‑3 (MQ=255)
                                     aaaCCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAACCCAATAGTTTGCGTCGCGATCCACAGGGTGCTGCTCATAatga                                                                                                                  >  2:1740583/1‑150 (MQ=255)
                                      aaCCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCACACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAatgat                                                                                                                 >  1:1614657/1‑150 (MQ=255)
                                      aaCCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAatgat                                                                                                                 <  2:520027/150‑1 (MQ=255)
                                       aCCTTAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAatgatg                                                                                                                <  2:1054874/150‑1 (MQ=255)
                                       aCCTGAATCAGTAACCATGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAatgatg                                                                                                                <  2:735219/150‑1 (MQ=255)
                                       aCCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAatgatg                                                                                                                <  1:647109/150‑1 (MQ=255)
                                        ccTGAATCAGTAAGCAGGTAATGGGAGAAGCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGGATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAAtgatga                                                                                                               <  2:1016460/150‑1 (MQ=255)
                                          tGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGacc                                                                                                             >  1:1597917/1‑149 (MQ=255)
                                          tGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGaca                                                                                                             >  1:1805268/1‑150 (MQ=255)
                                           gAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGacac                                                                                                            >  2:522549/1‑150 (MQ=255)
                                                aGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATCCCCCCGGCTTTTACCGCCCCGGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGGGTCGCGATCCCCAGGGTGATGCTCTTATTGATGCCAcccatg                                                                                                       >  2:870433/1‑147 (MQ=255)
                                                aGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTTCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCgag                                                                                                       >  1:431222/1‑150 (MQ=255)
                                                aGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCgag                                                                                                       <  1:764306/150‑1 (MQ=255)
                                                   aaCCAGGTAATGGGAGAACCGCCATGATACATTGAAGACACCATACCCCCGGCTTTTAACGCACCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTAGCGATCCACAGGGTGATGCTAATAATGATGACACCCGAGAtt                                                                                                    >  2:1134543/1‑150 (MQ=255)
                                                    aCCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGGGATGCTCATAATGATGACACCCGATATTa                                                                                                   <  2:1314480/150‑1 (MQ=255)
                                                     ccAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTaa                                                                                                  <  1:1512925/150‑1 (MQ=255)
                                                     ccAGGTAATGGGAGAACCGCCCTGATCCATCTAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTaa                                                                                                  >  2:1536081/1‑150 (MQ=255)
                                                       aGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAgg                                                                                                <  1:1593606/150‑1 (MQ=255)
                                                       aGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAgg                                                                                                <  1:942924/150‑1 (MQ=255)
                                                          tAATGGGAGAACCGCCCTGATCCATTGACGCCACCAAACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGGGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGc                                                                                             <  2:1620087/150‑1 (MQ=255)
                                                           taTGGGTGATTCGCCATGATCCATTGAAGCCACCATACCCCGGGCTTTTAACTCCCCAGTATTAGCCGCCAGTTGTTTGTGCAAACCATTTGTTTGCGTCGCGATCCACAGTGTGTTGCTCTTAATTATGACACCCGAGTTTAAGGTGCg                                                                                            <  2:1540943/149‑1 (MQ=255)
                                                           aaTGGGAGAACCGCCCTGATCCAATGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTAGGGTGTGCAAAGCAATAGTTTGCGTCGCGATGCACAGGGTGATGCTCATAATGATGACACCCGAGATTAACGTGCg                                                                                            >  2:127933/1‑150 (MQ=255)
                                                             tGGGAGAACCGCCCTGATCCATTGACGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGa                                                                                          >  1:108264/1‑150 (MQ=255)
                                                                gagaACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACa                                                                                       >  2:1443277/1‑150 (MQ=255)
                                                                    aCCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCACTAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGAtt                                                                                   <  2:913526/150‑1 (MQ=255)
                                                                    aCCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGAtt                                                                                   >  1:1732847/1‑150 (MQ=255)
                                                                      cGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTc                                                                                 <  1:1589403/150‑1 (MQ=255)
                                                                       gCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTcc                                                                                >  1:1186076/1‑150 (MQ=255)
                                                                         ccTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCTCCAGTGGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCAc                                                                              <  1:846259/150‑1 (MQ=255)
                                                                          cTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGCTGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACg                                                                             >  1:1414549/1‑150 (MQ=255)
                                                                            gATCCATTGAAGCCACCATACCCCCTGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGAt                                                                           <  1:219448/150‑1 (MQ=255)
                                                                              tCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGCTCCACCGGGTGTTGCTCATAATGATGACACCCGAGATTATGGAGCGGAAAAGATTTACACGATGc                                                                         >  1:901512/1‑150 (MQ=255)
                                                                              tCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGc                                                                         >  2:1032406/1‑150 (MQ=255)
                                                                               ccATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCa                                                                        >  1:1695719/1‑150 (MQ=255)
                                                                                cATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAAAAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGCCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCAc                                                                       >  1:742221/1‑150 (MQ=255)
                                                                                   tGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGc                                                                    >  2:136935/1‑150 (MQ=255)
                                                                                    gacgCCCCCATACCCCCGGCTTTTAACGCCCCAGCACTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGCGCGGAACAGATTTCCACGATGCACGGCg                                                                   <  2:959756/147‑1 (MQ=255)
                                                                                    gAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCATACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCg                                                                   >  1:753218/1‑150 (MQ=255)
                                                                                           ccaTACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGAGATGCTCATAATGATGACACCCGAGATTAAGGTGAGGAACAGATTTCCACGATGCACGGAGACCGCCa                                                            >  2:35872/1‑150 (MQ=255)
                                                                                              tACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCg                                                         >  1:603062/1‑150 (MQ=255)
                                                                                              aaCCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCg                                                         <  1:843269/149‑1 (MQ=255)
                                                                                                   ccGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACg                                                    >  1:1600713/1‑150 (MQ=255)
                                                                                                     ggCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACgaa                                                  <  2:1695719/150‑1 (MQ=255)
                                                                                                      gCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCTCCACgaag                                                 <  1:1295531/150‑1 (MQ=255)
                                                                                                        ttttAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCCCgaagca                                               >  2:166797/1‑148 (MQ=255)
                                                                                                          ttAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGc                                             >  2:511059/1‑150 (MQ=255)
                                                                                                           tAACGCCCCATAATTAGCCGCCAGTTGTGTGTGCAAACGAATAGTTTGCGTCGCGATCGACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGcc                                            <  2:866529/150‑1 (MQ=255)
                                                                                                           tAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCCCGAAGAAGcc                                            >  1:703863/1‑150 (MQ=255)
                                                                                                            aaCGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACCCCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCg                                           >  2:1312916/1‑150 (MQ=255)
                                                                                                             aCGCCCCAGCATTAGCCGCCAGTTGGGTGTGCATACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCAAAATTATGAAACACGAGATTAAGGTTCTGAACAGATTTCCCCGAGGCACGGCGACCGCCATCGCCACGAAGACGCCGa                                          >  2:879058/1‑150 (MQ=255)
                                                                                                              cGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGAt                                         >  1:497595/1‑150 (MQ=255)
                                                                                                                   cAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGAtggtgg                                    <  2:1132206/150‑1 (MQ=255)
                                                                                                                           gccgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGCTGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCg                            <  1:389369/150‑1 (MQ=255)
                                                                                                                           gccgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCg                            <  2:1251004/150‑1 (MQ=255)
                                                                                                                           gccgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCg                            >  2:610797/1‑150 (MQ=255)
                                                                                                                             cgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGcc                          >  1:738089/1‑150 (MQ=255)
                                                                                                                               ccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCaa                        <  1:1500229/150‑1 (MQ=255)
                                                                                                                                  gttgggttggCAAACCAATAGTTTGCGGCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGCTTAAGGTGCGGAACAGATTTCCACGATGCATGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACg                     <  1:1121341/150‑1 (MQ=255)
                                                                                                                                    tGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAATAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGgc                   >  1:1512329/1‑150 (MQ=255)
                                                                                                                                      ggTGTGCAAACCAATAGTTTGAGTCGCGATCAACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAAAAGATTTCCACGATGCACGGCGACCGCCATCGACACGAAGAAGCCGATGGTGGCAAGATAGCCAAACGgcag                 >  1:455195/1‑150 (MQ=255)
                                                                                                                                       cgacgCAAACTATTGGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGgcagc                >  1:1484842/5‑150 (MQ=255)
                                                                                                                                       cgacgCAAACTATTGGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGgcagc                >  2:763581/5‑150 (MQ=255)
                                                                                                                                       cgacgCAAACTATTGGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGgcagc                >  2:632156/5‑150 (MQ=255)
                                                                                                                                        gacgCAAACTATTGGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGgcagca               >  2:1266683/4‑150 (MQ=255)
                                                                                                                                        gacgCAAACTATTGGTTTGCGTCGCGATCCACAGGGTCATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGgcagcc               >  1:351055/4‑149 (MQ=255)
                                                                                                                                          tgCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGcccc             >  2:1097019/1‑150 (MQ=255)
                                                                                                                                           gCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCt            >  1:51936/1‑150 (MQ=255)
                                                                                                                                            cAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTg           <  1:892392/150‑1 (MQ=255)
                                                                                                                                            cAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTg           <  1:1420270/150‑1 (MQ=255)
                                                                                                                                             aaaCCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGa          <  2:1186076/150‑1 (MQ=255)
                                                                                                                                              aaCCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGAt         >  2:1030794/1‑150 (MQ=255)
                                                                                                                                              aaCCAATAGTTGGCGGCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGAt         <  1:1570738/150‑1 (MQ=255)
                                                                                                                                               aCCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGAtt        >  1:606442/1‑150 (MQ=255)
                                                                                                                                                     aGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGTCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCgg  <  2:1307046/150‑1 (MQ=255)
                                                                                                                                                     aGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGCTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCgg  >  1:246221/1‑150 (MQ=255)
                                                                                                                                                     aGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCgg  <  2:804995/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGG  >  NC_000913/2173006‑2173304

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