Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,173,155 A→G 8.2% pseudogene (1124/1358 nt) gatC ← galactitol‑specific PTS enzyme IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,1550AG 8.2% 241.0 / 10.2 85pseudogene (1124/1358 nt)gatCgalactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base A (38/40);  new base G (3/4);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

CGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCC  >  NC_000913/2173007‑2173302
                                                                                                                                                    |                                                                                                                                                   
cGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAg                                                                                                                                                    <  1:1056460/150‑1 (MQ=255)
cGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAg                                                                                                                                                    >  1:131903/1‑150 (MQ=255)
   taataaTGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGttt                                                                                                                                                 <  2:1216652/150‑1 (MQ=255)
     ataatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGATTCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACTATTGGTTTGc                                                                                                                                               <  1:813080/150‑1 (MQ=255)
     ataatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCGTTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGc                                                                                                                                               <  1:1156032/150‑1 (MQ=255)
      taatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCg                                                                                                                                              >  2:909217/1‑150 (MQ=255)
       aatGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACTATTGGTTTGcac                                                                                                                                             <  2:611170/150‑3 (MQ=255)
       aatGAAACCGGGAATATTTTGCGGGGAGAAAACCTCAATCAGTAACAAGGTAATGGGAGAACCGCCCTGATCAATTGAAGCCACCATACCCCCGGCTTATAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAAAAGTTTGCGt                                                                                                                                             >  2:668243/1‑150 (MQ=255)
        atGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCCCCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTc                                                                                                                                            >  2:1115024/1‑150 (MQ=255)
         tGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACTATTGGTTTGcacac                                                                                                                                           <  1:328707/150‑5 (MQ=255)
          gAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACTATTGGTTTGcacacc                                                                                                                                          <  2:215445/150‑6 (MQ=255)
          gAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTcgc                                                                                                                                          >  2:62865/1‑150 (MQ=255)
               cGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATcc                                                                                                                                     <  2:1129336/150‑1 (MQ=255)
                       ttttGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTg                                                                                                                             <  2:196492/150‑1 (MQ=255)
                       ttttGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTg                                                                                                                             <  2:1238996/150‑1 (MQ=255)
                        tttGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGa                                                                                                                            <  1:627662/150‑1 (MQ=255)
                           gCGGGGAGAAAACCTGAATCAGTACCCAGGTAATAGGAGAACCGCCCTGATACATTGAAGCCACCATACACCCGGATTTTAACGACCAAGCATAAGCCGCCAGTTGAGTGAGCAAACCAAAAGATTGCGTCCCCATCCACAGGGAGATGc                                                                                                                         >  2:604817/1‑150 (MQ=255)
                            cGGGGAGAAAACCTGAATCAGTCAGCAGGTAATGGGAGCACCGCCCTGATCCATTGACGCCCCCATGCCCCCGGCTGTTAACGCCCCAGCATTAGCCGCGAGTTGGGTGTGCACCCCAATAGGTTGCGGCGCGATCCCCAGGGTGAGGcg                                                                                                                        >  2:480461/1‑149 (MQ=255)
                                gagaAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATa                                                                                                                    >  1:1473127/1‑150 (MQ=255)
                                    aaaCCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAatga                                                                                                                >  2:13994/1‑150 (MQ=255)
                                     aaCCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAatgat                                                                                                               <  1:759619/150‑1 (MQ=255)
                                      aCCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAatgatg                                                                                                              <  1:1058503/150‑1 (MQ=255)
                                      aCCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAatgatg                                                                                                              <  2:404051/150‑1 (MQ=255)
                                         tGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGaca                                                                                                           <  1:670558/150‑1 (MQ=255)
                                              cAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCga                                                                                                      >  1:1194159/1‑150 (MQ=255)
                                               aGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAGCCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCgag                                                                                                     >  1:804860/1‑150 (MQ=255)
                                                 tAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGAt                                                                                                   <  1:705034/150‑1 (MQ=255)
                                                    ccAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTaa                                                                                                >  2:417966/1‑150 (MQ=255)
                                                      aGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCGGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAgg                                                                                              <  2:789893/150‑1 (MQ=255)
                                                      aGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAgg                                                                                              <  2:155040/150‑1 (MQ=255)
                                                         tAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGc                                                                                           >  2:685247/1‑150 (MQ=255)
                                                                agaACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAg                                                                                    <  2:900978/150‑1 (MQ=255)
                                                                   aCCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGAtt                                                                                 >  2:736831/1‑150 (MQ=255)
                                                                   aCCGCCATGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGAtt                                                                                 >  1:404031/1‑150 (MQ=255)
                                                                    ccGCCCTGCTCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCAGAATGATGCCACCCGCGATTAAGGTGCCGAACCGAttt                                                                                >  1:279021/1‑150 (MQ=255)
                                                                    ccGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGAttt                                                                                >  1:1486445/1‑150 (MQ=255)
                                                                       cccTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCa                                                                             >  1:725812/1‑150 (MQ=255)
                                                                        ccTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCAc                                                                            <  1:417966/150‑1 (MQ=255)
                                                                         cTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACg                                                                           >  1:952302/1‑150 (MQ=255)
                                                                           taTCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGAt                                                                         <  2:1087289/149‑1 (MQ=255)
                                                                            aTCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATg                                                                        >  1:219359/1‑150 (MQ=255)
                                                                             tCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGc                                                                       <  2:1194159/150‑1 (MQ=255)
                                                                             tCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGc                                                                       <  2:351175/150‑1 (MQ=255)
                                                                              ccATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCa                                                                      <  2:1473127/150‑1 (MQ=255)
                                                                              ccATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCa                                                                      <  1:909217/150‑1 (MQ=255)
                                                                              ccATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCa                                                                      <  1:719712/150‑1 (MQ=255)
                                                                                  tGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCTATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGc                                                                  <  2:1145306/150‑1 (MQ=255)
                                                                                  tGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGAGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGc                                                                  <  2:51495/150‑1 (MQ=255)
                                                                                        caccaTACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGc                                                            >  1:438700/1‑150 (MQ=255)
                                                                                         accaTACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGcc                                                           >  1:303508/1‑150 (MQ=255)
                                                                                            aTACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATc                                                        >  2:890956/1‑150 (MQ=255)
                                                                                            aTACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATc                                                        <  1:62865/150‑1 (MQ=255)
                                                                                               cccccGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAatgatg                                                                                                              <  2:973975/93‑1 (MQ=255)
                                                                                               cccccGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAatgatg                                                                                                              >  1:973975/1‑93 (MQ=255)
                                                                                                ccccGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCa                                                    <  2:1446154/150‑1 (MQ=255)
                                                                                                  ccGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACg                                                  <  2:1105160/150‑1 (MQ=255)
                                                                                                        tttAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGaagaag                                            <  2:952302/150‑1 (MQ=255)
                                                                                                        tttAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGaagaag                                            >  2:1278156/1‑150 (MQ=255)
                                                                                                        tttAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGaagaag                                            >  2:1187059/1‑150 (MQ=255)
                                                                                                          tAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACgaaga                                              <  1:1389108/146‑1 (MQ=255)
                                                                                                          tAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACgaaga                                              >  2:1389108/1‑146 (MQ=255)
                                                                                                             cGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAACCCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGAt                                       >  1:1101417/1‑150 (MQ=255)
                                                                                                             cGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGACTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGAt                                       <  1:480461/150‑1 (MQ=255)
                                                                                                              gCCCCAGCATTAGCCGCCTGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCAGAGGGTGATGCTCATAATGATGACAGCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCAGGGCGACCGGCATCGCCACGATGATGCCGAtg                                      <  2:860793/150‑1 (MQ=255)
                                                                                                               ccccAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGAtgg                                     >  2:719318/1‑150 (MQ=255)
                                                                                                                 ccAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGAtggtg                                   <  1:736831/150‑1 (MQ=255)
                                                                                                                        tAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGAt                            >  2:809586/1‑150 (MQ=255)
                                                                                                                            cgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGcc                        >  1:820972/1‑150 (MQ=255)
                                                                                                                              ccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCaa                      <  1:685247/150‑1 (MQ=255)
                                                                                                                              ccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCaa                      <  2:774238/150‑1 (MQ=255)
                                                                                                                               cAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCaaa                     <  1:772399/150‑1 (MQ=255)
                                                                                                                                   tGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGgc                 <  2:219359/150‑1 (MQ=255)
                                                                                                                                    gggTGTGCAATCCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGgca                <  2:804860/150‑1 (MQ=255)
                                                                                                                                    gggTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGgca                >  2:181568/1‑150 (MQ=255)
                                                                                                                                    cgcgacgCAAACTATTGGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGgca                >  2:1078843/7‑150 (MQ=255)
                                                                                                                                    cgcgacgCAAACTATTGGTTTGAGTCGCGATCCACAGGGTGATGCTCATAATGCTGACACCCGAGATTAAGGTGCGGAACAGATTTACACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGCTGGCAAGATCGCCAAACGgcc                >  2:221616/7‑149 (MQ=255)
                                                                                                                                      gtgtGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGgcagc              >  1:775164/1‑150 (MQ=255)
                                                                                                                                      cgacgCAAACTATTGGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGgcagc              >  1:128110/5‑150 (MQ=255)
                                                                                                                                          gCATACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCATGATCGCCAAACGGCAGCACCt          <  2:775164/150‑1 (MQ=255)
                                                                                                                                              aCCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGAtt      >  1:113171/1‑150 (MQ=255)
                                                                                                                                               ccAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTa     <  1:1154503/150‑1 (MQ=255)
                                                                                                                                               ccAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTa     >  2:873177/1‑150 (MQ=255)
                                                                                                                                                cAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTAc    >  1:1225952/1‑150 (MQ=255)
                                                                                                                                                 aaTAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTAcc   >  2:1472068/1‑150 (MQ=255)
                                                                                                                                                  aTAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTAccc  >  2:457506/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGATAATAATGAAACCGGGAATATTTTGCGGGGAGAAAACCTGAATCAGTAACCAGGTAATGGGAGAACCGCCCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCC  >  NC_000913/2173007‑2173302

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