Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,112,175 A→G 7.3% pseudogene (767/849 nt) yghF ← pseudogene, secretion pathway protein, C‑type protein homology;putative transport; Not classified; putative general secretion pathway for protein export (GSP)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,112,1750AG 7.3% 207.1 / 3.8 80pseudogene (767/849 nt)yghFpseudogene, secretion pathway protein, C‑type protein homology;putative transport; Not classified; putative general secretion pathway for protein export (GSP)
Reads supporting (aligned to +/- strand):  ref base A (43/31);  new base G (3/3);  total (46/34)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGCCCCCGGCTTCGCTGCATAACCGACAATCCCCTCCTTACGCA  >  NC_000913/3112029‑3112319
                                                                                                                                                  |                                                                                                                                                
aaaTGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCGATTGAAGGTAACTGCCGCATc                                                                                                                                                >  2:38272/1‑149 (MQ=255)
   tGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCaga                                                                                                                                             >  1:164849/1‑149 (MQ=255)
       gatCATCGAAATCCTGGTTCAGCTCTTCTTCCAGTAAGCCGGGCAGAATTGTCGGTAGCGCACCGCGATGGAAATGTCGTGGCGCGCACCCTTCCGAAAAACCGTAAGTTGAATGCAATCCATAGAAGGTAACTGCCGCGTCAGAGCaa                                                                                                                                         >  1:434649/3‑149 (MQ=255)
       gTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCGTCAGAGCaa                                                                                                                                         >  1:438580/1‑149 (MQ=255)
        ttCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAACGCGATGGAAATGTCGTTGCTCGCACCCTTTCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGATCAAt                                                                                                                                        >  2:879211/1‑149 (MQ=255)
              gATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTGGAATGGAATCCAGTGAAGGTAACTGCCGCGTCAGAGCAATAATTGc                                                                                                                                  >  2:290793/1‑149 (MQ=255)
                tatCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCCGCGCGATGGAAATGTCGTGGCGCGCACCCGTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGCGCAATCATTGCTc                                                                                                                                >  1:368772/1‑149 (MQ=255)
                   ccTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGAAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATACATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTg                                                                                                                             >  2:926071/1‑149 (MQ=255)
                    cTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCGCCCGTGCGTTAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTgg                                                                                                                            >  2:343095/1‑149 (MQ=255)
                    cTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTgg                                                                                                                            >  1:168571/1‑149 (MQ=255)
                            gCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGa                                                                                                                    >  1:774819/1‑149 (MQ=255)
                            gCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGa                                                                                                                    >  1:157630/1‑149 (MQ=255)
                              tcttctTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGaaa                                                                                                                  >  1:691520/1‑149 (MQ=255)
                              gcatcttcCAGTAAGGCGGGCATAAGAGTCGGGAGAGCAGCGCGATGGAAATGGCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGaaa                                                                                                                  <  2:1015984/146‑1 (MQ=255)
                               cttcttCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAAt                                                                                                                 >  1:428944/1‑149 (MQ=255)
                               cttcttCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTGGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCTCATCAGAGCAATCATTGCTCGTGGATCAGTGAAAt                                                                                                                 >  2:24272/1‑149 (MQ=255)
                                ttcttcCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATc                                                                                                                <  2:957617/149‑1 (MQ=255)
                                ttcttcCAGGCAGGCGGGCAGAATAGTCGGTAGCGCCGTGCGATGGTCCTGTCGTTGCGCGCTCCCTTGCGTAAAACCGTTAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCCTTGCTCGTGGATCAGTGAAATc                                                                                                                <  2:927213/149‑1 (MQ=255)
                                    tcCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCctgc                                                                                                            >  1:94903/1‑149 (MQ=255)
                                     ccAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGTACCCTTGCGTAAATCCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCATAGCAATCATTGCTCGTGGATCAGTGAAATcctcct                                                                                                           >  2:559375/1‑149 (MQ=255)
                                       aGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGa                                                                                                         <  2:129411/149‑1 (MQ=255)
                                         tAAGGCGGGCAGAATAGTTGGAAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGAtt                                                                                                       <  1:70129/149‑1 (MQ=255)
                                         tAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGAtt                                                                                                       <  2:503316/149‑1 (MQ=255)
                                             gCGGGCAGAATATTCTGTATCGCAGCGCGATGGAAATGTCGTTGCGCGCACCCTTGCTTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGACTAATCAGAGCAATCATTTCTCGTGGATCATTGAAATCATGCTGATTTAgc                                                                                                   <  2:115526/149‑1 (MQ=255)
                                             gCGGGCAGAAGAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAgc                                                                                                   <  2:660107/149‑1 (MQ=255)
                                             gCGGCCAGCATAGTCGGTAGCGCAGCGCGATGGAACTGTCGTGCCGCGCACCCGTGCGGAAAACCGTCAGTTGAATGGAATCCATTGAAGGTAACTGCCGCGGCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCGGATTTAgc                                                                                                   <  2:578612/149‑1 (MQ=255)
                                              cGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAgcg                                                                                                  <  2:164849/149‑1 (MQ=255)
                                               gggCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGATTCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAgcgc                                                                                                 >  1:278964/1‑149 (MQ=255)
                                                   aGATTAGTCGGGAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCTTAAGTTGATTGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATg                                                                                             <  1:223011/149‑1 (MQ=255)
                                                   aGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATg                                                                                             <  2:691520/149‑1 (MQ=255)
                                                    gAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAGCCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATgg                                                                                            >  2:273071/1‑149 (MQ=255)
                                                     aaTAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGc                                                                                           >  1:386100/1‑149 (MQ=255)
                                                     aaTAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGAACGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGc                                                                                           >  2:959306/1‑149 (MQ=255)
                                                         gTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTACCTGCCGCATCAGAGCAATCATTGCTCGTCGATCAGTGAATTCCTGCTGATTTAGCGCAATGGCGata                                                                                       >  1:1012106/1‑149 (MQ=255)
                                                             gTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGc                                                                                   >  2:258625/1‑149 (MQ=255)
                                                                 cgcagcgcGATGGATATGTCGTGGCGCGCACGGTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAACTCCTGCGGATTTAGCGCAATGGCGATATCGccttc                                                                               <  2:645609/149‑1 (MQ=255)
                                                                  gcagcgcGATGGAAATGTCGGGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGAAATCGccttcc                                                                              >  1:452489/1‑149 (MQ=255)
                                                                   cagcgcGATGGAAATGTCGTTGCGCGCACACTTGAGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAAATGCCTCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGccttcct                                                                             <  1:883852/149‑1 (MQ=255)
                                                                   cagcgcGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGccttcct                                                                             >  1:855461/1‑149 (MQ=255)
                                                                    agcgcGATGGAAATGTCGTGGCGCGAACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCTCAATGGCGATATCGccttcctt                                                                            <  2:582473/149‑1 (MQ=255)
                                                                    agcgcGATGGAAATGTCGGGGCGCGCACCCTTGCGTAAACCCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCAATGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGccttcctt                                                                            <  2:631132/149‑1 (MQ=255)
                                                                    agcgcGATGGAAATGTCGAGGCGCGCACCCTTGCGTAAAACCGTAAGTTGCATGGAATCCAATGAAGGTAACAGCCGCATCAGAGCACTCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGccttcctt                                                                            <  2:678084/149‑1 (MQ=255)
                                                                     gcgcGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCACTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTg                                                                           >  2:525732/1‑149 (MQ=255)
                                                                       gcgATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGGAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGAATTAGCGCAATCGCGATATCGCCGTCCTTGaa                                                                         >  2:216304/1‑149 (MQ=255)
                                                                        cgATGGAAATGTCGTGGCGCGCACCCCTGCGAAAAACCGTCAGCGGACTCGAATCCATTGAAGGTCACTGCCCCATCAGAGCAATCATTCCGCGTGGATCTGTGAAATCCTGCTGATTTAGCGCAACGGCGATATCGCCTTCCTTGaaa                                                                        <  2:496207/149‑1 (MQ=255)
                                                                         gATGGAAATGTCGGGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAAc                                                                       <  1:677637/149‑1 (MQ=255)
                                                                            ggAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCt                                                                   >  1:649228/1‑149 (MQ=255)
                                                                               aaTGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAACTCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTgg                                                                 >  2:803061/1‑149 (MQ=255)
                                                                                 tGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCa                                                               >  1:463303/1‑149 (MQ=255)
                                                                                  gTCGGGGCGCGCACCCTTGCGTAAAACCGTACGTAGAAAGGCATCCCTTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCGTCCTTGAAACCGCTGGCAt                                                              <  2:582913/149‑1 (MQ=255)
                                                                                   tCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATc                                                             >  2:159150/1‑149 (MQ=255)
                                                                                      tGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGaa                                                          <  1:366525/149‑1 (MQ=255)
                                                                                      gggCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGACTCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGaa                                                          <  2:656235/148‑1 (MQ=255)
                                                                                      gggCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGaa                                                          <  2:94903/148‑1 (MQ=255)
                                                                                         cgcgcACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACag                                                       <  2:140185/149‑1 (MQ=255)
                                                                                          gagcACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATAGCTCGTGGATAAGTGATATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTACTTGAAACCGCTGGCCTCGAACaga                                                      <  2:69137/147‑1 (MQ=255)
                                                                                                ccTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGa                                                <  2:168571/149‑1 (MQ=255)
                                                                                                  ttGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATc                                              >  1:757247/1‑149 (MQ=255)
                                                                                                     cGTAAAACGGTAAGTATTATGGAATCGATTGTAGGTAACTGCCGCGTCAGTGGAATCATTGCTCGTGGATCAGTGTAATCCTGCTTATTTTGCGCAATGGCGATATCGCCTTCCTTGAATCCGCTGGCATCGAACAGAGAACGATCTGc                                           <  2:434649/149‑1 (MQ=255)
                                                                                                     cGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCGTCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGc                                           <  2:438580/149‑1 (MQ=255)
                                                                                                      gTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGcc                                          >  1:585396/1‑149 (MQ=255)
                                                                                                      gTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGcc                                          >  1:822533/1‑149 (MQ=255)
                                                                                                      gTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGcc                                          >  1:425410/1‑149 (MQ=255)
                                                                                                        aaaaCCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGcccc                                        <  2:309442/149‑1 (MQ=255)
                                                                                                         aaaCCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGccccc                                       >  2:416840/1‑149 (MQ=255)
                                                                                                           aCCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGCCCCCgg                                     <  1:86166/149‑1 (MQ=255)
                                                                                                             cGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGCCCCCGGCt                                   >  2:727342/1‑149 (MQ=255)
                                                                                                                aaGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGCCCCCGGCTTCg                                <  1:273071/149‑1 (MQ=255)
                                                                                                                  gTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATAGAACAGAGAACGATCTGCCCCCGGCTTCGCt                              <  1:926071/149‑1 (MQ=255)
                                                                                                                         ggAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGCCCCCGGCTTCGCTGCATAAc                       <  2:797328/149‑1 (MQ=255)
                                                                                                                         ggAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGCCCCCGGCTTCGCTGCATAAc                       >  1:71870/1‑149 (MQ=255)
                                                                                                                          gAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGCCCCCGGCTTCGCTGCATAAcc                      >  2:505686/1‑149 (MQ=255)
                                                                                                                               cATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGCCCCCGGCTTCGCTGCATAACCGACaa                 >  2:824696/1‑149 (MQ=255)
                                                                                                                                aTTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGCCCCCGGCTTCGCTGCATAACCGACAAt                >  2:529925/1‑149 (MQ=255)
                                                                                                                                  tGAAGGTAAATGCCGCATCAGAGCAATCATTGCTGGTGGATCTGTGAAATCCTGCTGATTTAGCGCAATCGCGAAATCGCCTTCCGTGAACCCGCTGGCATCGAGCAGAGACCGCTCTGCCCCCGCCTTCCCTGCATAACCGACAAGcc              >  1:686530/1‑149 (MQ=255)
                                                                                                                                     aGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGCCCCCGGCTTCGCTGCATAACCGACAATcccca           >  1:504108/1‑148 (MQ=255)
                                                                                                                                            tGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGCCCCCGGCTTCGATGCATAACAGACAATCCCCTCCTTAcc    >  1:272792/1‑148 (MQ=255)
                                                                                                                                             gCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGCCCCCGGCTTCGCTGCATAACCGACAATCCCCTCCTTACGc   >  2:585580/1‑149 (MQ=255)
                                                                                                                                              ccGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGCCCCCGGCTTCGCTGCATAACCGACAATCCCCTCCTTACGca  <  2:774819/149‑1 (MQ=255)
                                                                                                                                              ccGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGCCCCCGGCTTCGCTGCATAACCGACAATCCCCTCCTTACGca  >  2:420051/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
AAATGCAGTTCATCGATATCCTGGTTCAGCTCTTCTTCCAGTAAGGCGGGCAGAATAGTCGGTAGCGCAGCGCGATGGAAATGTCGTGGCGCGCACCCTTGCGTAAAACCGTAAGTTGAATGGAATCCATTGAAGGTAACTGCCGCATCAGAGCAATCATTGCTCGTGGATCAGTGAAATCCTGCTGATTTAGCGCAATGGCGATATCGCCTTCCTTGAAACCGCTGGCATCGAACAGAGAACGATCTGCCCCCGGCTTCGCTGCATAACCGACAATCCCCTCCTTACGCA  >  NC_000913/3112029‑3112319

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: