Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,068,437 T→A 5.1% intergenic (+183/+74) ymdF → / ← rutG KGG family protein/pyrimidine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,068,4370TA 5.1% 290.5 / 3.8 98intergenic (+183/+74)ymdF/rutGKGG family protein/pyrimidine permease
Reads supporting (aligned to +/- strand):  ref base T (49/44);  new base A (3/2);  total (52/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

Only 100 of 101 total aligned reads displayed.

GTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATC  >  NC_000913/1068288‑1068586
                                                                                                                                                     |                                                                                                                                                     
gTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGTACTGTTTTGTAGGCCGGATAAGGCGttt                                                                                                                                                       <  1:1513616/150‑1 (MQ=255)
  aGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTAc                                                                                                                                                     >  2:366960/1‑150 (MQ=255)
   gCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTAcg                                                                                                                                                    >  1:418313/1‑150 (MQ=255)
    cTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCATAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                                   >  2:1710087/1‑150 (MQ=255)
    cTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                                   <  2:650553/150‑1 (MQ=255)
      ggTCAAAAAGCCATCACCGCCGCTTACCCGGTGGTTGTTACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGGAATATATTGAATTTGCACTGTTTTGTAGGCCGGTTAAGGCGTTTAcgccg                                                                                                                                                 <  2:1550970/149‑1 (MQ=255)
      cggcAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                                                                 <  2:771144/147‑1 (MQ=255)
        tCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTAACAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCa                                                                                                                                               <  2:798616/150‑1 (MQ=255)
           aaaaGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATcc                                                                                                                                            >  1:213111/1‑150 (MQ=255)
              aGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                         >  1:1041661/1‑150 (MQ=255)
               gCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                        <  1:580738/150‑2 (MQ=37)
               gCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                        <  2:1636963/150‑1 (MQ=255)
               gCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                        >  2:876760/1‑150 (MQ=255)
                ccATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcct                                                                                                                                       <  2:4936/150‑3 (MQ=37)
                 cATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATg                                                                                                                                      >  1:894666/1‑150 (MQ=255)
                 cATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATg                                                                                                                                      >  2:632726/1‑150 (MQ=255)
                    cACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAc                                                                                                                                   >  2:1172842/1‑150 (MQ=255)
                     aCCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCTCATCCGGCATGAACa                                                                                                                                  >  2:329660/1‑150 (MQ=255)
                      ccgccgGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACaa                                                                                                                                 >  1:1676009/1‑150 (MQ=255)
                         ccgGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGc                                                                                                                              >  2:635397/1‑150 (MQ=255)
                          cgGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCg                                                                                                                             <  1:1191310/150‑1 (MQ=255)
                           ggTTACCCGGAGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAAACTATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGt                                                                                                                            <  1:136748/150‑1 (MQ=255)
                               aCCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGt                                                                                                                        <  2:1000785/150‑1 (MQ=255)
                                cccGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGtt                                                                                                                       >  1:1683105/1‑150 (MQ=255)
                                      ggtTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTcaac                                                                                                                 >  1:500707/1‑150 (MQ=255)
                                      ggtTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTcaac                                                                                                                 >  2:1244674/1‑150 (MQ=255)
                                          gATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATc                                                                                                             >  2:1539691/1‑150 (MQ=255)
                                             aCTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGc                                                                                                          <  1:1559201/150‑1 (MQ=255)
                                              cTGATGCCAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCa                                                                                                         >  1:102350/1‑150 (MQ=255)
                                              cTGATGACAAATGTAAGCTTGCCTGATGCGCGCTGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACACTGCGTACGTTGTCAACAATCTGCa                                                                                                         >  1:106964/1‑150 (MQ=255)
                                              cTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCa                                                                                                         >  1:155622/1‑150 (MQ=255)
                                                  tgaCAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCAccgc                                                                                                     <  1:1527748/150‑1 (MQ=255)
                                                    aCAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCAccgccg                                                                                                   >  1:502159/1‑150 (MQ=255)
                                                     cAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCgg                                                                                                  >  1:605212/1‑150 (MQ=255)
                                                      aaaTGTAAGCTGGCCTGATGCGCGATGCTTATCAGTCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGt                                                                                                 <  2:180936/150‑1 (MQ=255)
                                                       aaTGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCAGTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGGCTCCGGCATGAACTATGCGTACGTTGTCATCAATCTGCTCCGCCGGTa                                                                                                >  2:236040/1‑150 (MQ=21)
                                                       aaTGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTa                                                                                                <  1:401850/150‑1 (MQ=255)
                                                           tAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAAccc                                                                                            >  1:1224944/1‑150 (MQ=255)
                                                            aaGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAAcccc                                                                                           >  1:209398/1‑150 (MQ=255)
                                                            aaGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAAcccc                                                                                           >  2:1276537/1‑150 (MQ=255)
                                                             aGCTTGCCTGATGCTCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCcg                                                                                          <  1:771815/150‑1 (MQ=255)
                                                             aGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCcg                                                                                          >  1:1476667/1‑150 (MQ=255)
                                                                  gCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGt                                                                                     <  1:790758/150‑1 (MQ=255)
                                                                   ccTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGtt                                                                                    <  2:373350/150‑1 (MQ=255)
                                                                    cTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGttt                                                                                   <  1:1379580/150‑1 (MQ=255)
                                                                        tGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTg                                                                               <  1:876760/150‑1 (MQ=255)
                                                                         gcgcGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGt                                                                              >  2:1400338/1‑150 (MQ=255)
                                                                              aTGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATg                                                                         <  1:635397/150‑1 (MQ=255)
                                                                              aTGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATg                                                                         <  2:1725091/150‑1 (MQ=255)
                                                                              aTGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATg                                                                         <  2:894666/150‑1 (MQ=255)
                                                                                gCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGc                                                                       >  1:1099292/1‑150 (MQ=255)
                                                                                 cTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCt                                                                      >  2:1626829/1‑150 (MQ=255)
                                                                                   tATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTcc                                                                    >  1:401643/1‑150 (MQ=255)
                                                                                   tATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTcc                                                                    >  1:741928/1‑150 (MQ=255)
                                                                                    aTCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCt                                                                   >  2:7097/1‑150 (MQ=255)
                                                                                    aTCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCt                                                                   >  2:112553/1‑150 (MQ=255)
                                                                                    aTCAGGCCTACCAGAAGATTGCAAAATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCt                                                                   <  1:258997/150‑1 (MQ=255)
                                                                                       aGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTCACCCCCGCGGTTTTCTGTTTATGGATCCtgat                                                                >  2:330300/1‑150 (MQ=255)
                                                                                       aGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCtgat                                                                <  1:1539691/150‑1 (MQ=255)
                                                                                       aGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCtgat                                                                >  1:1683507/1‑150 (MQ=255)
                                                                                       aGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCtgat                                                                >  2:175243/1‑150 (MQ=255)
                                                                                        ggCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCtgatg                                                               <  1:667348/150‑1 (MQ=255)
                                                                                        ggCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCtgatg                                                               <  2:241042/150‑1 (MQ=255)
                                                                                         gCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCtgatga                                                              <  2:1131989/150‑1 (MQ=255)
                                                                                          ccTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGaa                                                             >  2:599849/1‑150 (MQ=255)
                                                                                            tACCAGAAGATTGCAATATTTTGAATTTGCACTGTTTTGTATGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATTCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGaaca                                                           <  2:212734/150‑1 (MQ=255)
                                                                                              ccAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGaacaac                                                         >  2:907979/1‑150 (MQ=255)
                                                                                                   agaTTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTg                                                    <  1:1042752/150‑1 (MQ=255)
                                                                                                     aTTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGc                                                  >  2:46114/1‑150 (MQ=255)
                                                                                                      ttGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCg                                                 <  1:165619/150‑1 (MQ=255)
                                                                                                       tGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTATGTTTATTGGTACAGATGAACAACTTATGGCgg                                                >  1:1428987/1‑150 (MQ=25)
                                                                                                        gCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGt                                               >  2:1677387/1‑150 (MQ=255)
                                                                                                        gCAATATATTGAATTTGCACTGGTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGt                                               <  2:1213781/150‑1 (MQ=255)
                                                                                                          aaTATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTgg                                             <  1:1031131/150‑1 (MQ=255)
                                                                                                          aaTATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTgg                                             >  2:1683699/1‑150 (MQ=255)
                                                                                                              tatTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCGGGCGGTGGAACg                                         >  1:1667814/1‑150 (MQ=255)
                                                                                                                 tGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTca                                      <  2:1659127/150‑1 (MQ=255)
                                                                                                                 tGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTca                                      <  2:1598550/150‑1 (MQ=255)
                                                                                                                  taaTTTGCACTGTTTTTTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTTATGAACAACTTCTGGCGGTGGAACGTcaa                                     <  2:634636/149‑1 (MQ=255)
                                                                                                                       tGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTcaaccaac                                <  2:826261/150‑1 (MQ=255)
                                                                                                                       tGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTcaaccaac                                <  1:1712676/150‑1 (MQ=255)
                                                                                                                             gTTTTGTAGGCCGGATAAGGCGTTTCCGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGGTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGACCGTCATCCAACGTTCTg                          >  1:836004/1‑150 (MQ=38)
                                                                                                                                    atgCCGGATGCGGCGTAAACGGCTTATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTaaca                   >  1:236272/3‑150 (MQ=255)
                                                                                                                                     ggCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTaacaa                  >  2:410207/1‑150 (MQ=255)
                                                                                                                                     ggCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTaacaa                  >  2:1192764/1‑150 (MQ=255)
                                                                                                                                      gacGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTaacaac                 >  1:1727113/3‑150 (MQ=255)
                                                                                                                                      gCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTaacaac                 >  1:1727117/1‑150 (MQ=255)
                                                                                                                                       ccGGATGCGGCGTAAACGCCTTATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGATCAACGTCTGGCGGTGGCACGTCCCCCAACTTTCTGCTTTACAACg                >  1:1414258/1‑150 (MQ=255)
                                                                                                                                       ccGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTAGGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACg                <  2:102350/150‑1 (MQ=255)
                                                                                                                                       ccGGAGAAGGCGTTTACGCCGCATCCGGCATGAACAATTCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACg                <  2:106964/150‑1 (MQ=255)
                                                                                                                                        cGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTGTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGc               <  1:51412/150‑1 (MQ=255)
                                                                                                                                         ggATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCa              <  1:632726/150‑1 (MQ=255)
                                                                                                                                            tAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTg           <  1:112553/150‑1 (MQ=255)
                                                                                                                                              aGGCGTTTACGCCGCATCCTGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAg         <  1:329660/150‑1 (MQ=255)
                                                                                                                                              aGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAg         <  2:418313/150‑1 (MQ=255)
                                                                                                                                               ggCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGt        <  2:1676009/150‑1 (MQ=255)
                                                                                                                                               ggCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGt        <  1:175243/150‑1 (MQ=255)
                                                                                                                                                 cGTTTACGCCGCATCCGGCATGAACAATGCGTTCGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGGTTTCTGTTTATGGCTCCTGAGGAACAACTTCTGGCGGTGGATCGTCACCCACCTGTCTGCTTAACAACGCAGTGAGTaa      >  2:929745/1‑150 (MQ=255)
                                                                                                                                                     tACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTACCCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACCACTTCTGGCGGTGGCACGTCCACCAACTTTCTGCTTAACAACGCATTGAGTAAAATc  >  1:450515/1‑150 (MQ=255)
                                                                                                                                                     tACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACATCAGGCGGAGGAACGTCCACCAACTTTGTGCTTAACAACGCATTGAGTAAAATc  >  1:1718996/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAACAATGCGTACGTTGTCAACAATCTGCACCGCCGGTAACCCCGGCGGTTTTCTGTTTATGGCTCCTGATGAACAACTTCTGGCGGTGGAACGTCAACCAACTTTCTGCTTAACAACGCATTGAGTAAAATC  >  NC_000913/1068288‑1068586

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