Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 691,390 G→C 5.7% A132G (GCC→GGC)  ybeX ← CorC‑HlyC family protein YbeX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913691,3900GC 5.7% 293.4 / 7.3 87A132G (GCC→GGC) ybeXCorC‑HlyC family protein YbeX
Reads supporting (aligned to +/- strand):  ref base G (40/42);  new base C (3/2);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.07e-01

AATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGTTTC  >  NC_000913/691242‑691539
                                                                                                                                                    |                                                                                                                                                     
aaTGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTCGTCCATGCTGAAGGc                                                                                                                                                      <  1:1791336/150‑1 (MQ=255)
   ggTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGTAACGACAACCGCCTGACGTAACACTTTGTCCATGCTTAAGGCTTc                                                                                                                                                   <  2:1989456/150‑1 (MQ=255)
    gTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCa                                                                                                                                                  >  1:1365870/1‑150 (MQ=255)
      cACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGc                                                                                                                                                <  1:1487106/150‑1 (MQ=255)
        ccAGACCGGAAACCCCACCGAATTCGTCAATAACCATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCAt                                                                                                                                              >  2:1214055/1‑150 (MQ=255)
         cAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCTTTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAtg                                                                                                                                             <  2:924289/150‑2 (MQ=255)
         cAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATc                                                                                                                                             >  2:870190/1‑150 (MQ=255)
           gACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGc                                                                                                                                           >  2:468295/1‑150 (MQ=255)
           gACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAtgga                                                                                                                                           <  2:1243089/150‑4 (MQ=255)
             ccGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTg                                                                                                                                         >  1:1312557/1‑150 (MQ=255)
             ccGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTAAAGGCTTCAGCATCGCTg                                                                                                                                         >  1:1987461/1‑150 (MQ=255)
              cGGAAACCCCACCGAATTCGTCTATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTgc                                                                                                                                        <  1:1357521/150‑1 (MQ=255)
              cGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCAGGCTGAAGGCTTCAGCATCGCTgc                                                                                                                                        <  1:584075/150‑1 (MQ=255)
               ggAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTgcg                                                                                                                                       <  2:942727/150‑1 (MQ=255)
               ggAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTgcg                                                                                                                                       <  1:1872641/150‑1 (MQ=255)
               ggAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGGAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTgcg                                                                                                                                       >  1:1047542/1‑150 (MQ=255)
                gAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTgcgc                                                                                                                                      >  1:96561/1‑150 (MQ=255)
                gAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTgcgc                                                                                                                                      <  2:1369792/150‑1 (MQ=255)
                 aaaCCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATTGCTGCGCa                                                                                                                                     <  1:1552880/150‑1 (MQ=255)
                 aaaCCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCa                                                                                                                                     <  1:842276/150‑1 (MQ=255)
                                       aaCGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGc                                                                                                               <  1:1098256/150‑1 (MQ=255)
                                           aTCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATc                                                                                                           >  2:910190/1‑150 (MQ=255)
                                                  tgtgGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATcc                                                                                                    >  2:365061/1‑150 (MQ=255)
                                                      gTAACGCTGAGAGCGAAACTCTTTTAGCATCCGGTCTACGAGCTTACTTTCAGGAACGACAACCGCCTGACTTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTc                                                                                                <  1:869562/150‑1 (MQ=255)
                                                       tAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCa                                                                                               >  2:1669537/1‑150 (MQ=255)
                                                       tAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCa                                                                                               >  1:515009/1‑150 (MQ=255)
                                                          cGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCGTCAAtg                                                                                            >  2:1753774/1‑150 (MQ=255)
                                                                gagCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCATCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCt                                                                                      <  2:698830/150‑1 (MQ=255)
                                                                 agCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCtt                                                                                     >  1:12006/1‑150 (MQ=255)
                                                                  gCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCttt                                                                                    >  2:335056/1‑150 (MQ=255)
                                                                   cGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTg                                                                                   >  2:914596/1‑150 (MQ=255)
                                                                   cGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTg                                                                                   <  2:1047542/150‑1 (MQ=255)
                                                                   cGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCCTTg                                                                                   <  1:870190/150‑1 (MQ=255)
                                                                     aaaCTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTc                                                                                 <  2:1312557/150‑1 (MQ=255)
                                                                     aaaCTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTc                                                                                 >  2:880556/1‑150 (MQ=255)
                                                                     aaaCTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCCGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGa                                                                                         <  2:331915/142‑1 (MQ=255)
                                                                     aaaCTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCCGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGa                                                                                         >  1:331915/1‑142 (MQ=255)
                                                                         tctTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACGCTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCg                                                                             >  2:1948506/1‑150 (MQ=255)
                                                                         tctTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTAAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCg                                                                             <  2:1987461/150‑1 (MQ=255)
                                                                          ctTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGc                                                                            <  1:382392/150‑1 (MQ=255)
                                                                          ctTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGc                                                                            >  2:1219412/1‑150 (MQ=255)
                                                                           tttCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCt                                                                           >  1:776431/1‑150 (MQ=255)
                                                                              cAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAAt                                                                        >  1:201065/1‑150 (MQ=255)
                                                                               aGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATc                                                                       >  1:1913424/1‑150 (MQ=255)
                                                                               aGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATc                                                                       <  1:791930/150‑1 (MQ=255)
                                                                                gCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAAtg                                                                                            <  1:215235/128‑1 (MQ=255)
                                                                                gCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAAtg                                                                                            >  2:215235/1‑128 (MQ=255)
                                                                                gCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCa                                                                      <  1:884299/150‑1 (MQ=255)
                                                                                 cATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCAc                                                                     >  2:1957489/1‑150 (MQ=255)
                                                                                    ccGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAATGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCgg                                                                  >  2:1758756/1‑150 (MQ=255)
                                                                                    ccGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCgg                                                                  <  1:914596/150‑1 (MQ=255)
                                                                                     cGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCggg                                                                 <  2:1624004/150‑1 (MQ=255)
                                                                                      ggTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGa                                                                >  2:1353634/1‑150 (MQ=255)
                                                                                        tCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGaaa                                                              <  2:1063133/150‑1 (MQ=255)
                                                                                         cTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAAc                                                             >  2:1774172/1‑150 (MQ=255)
                                                                                            cgcgCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTg                                                          <  1:1787169/150‑1 (MQ=255)
                                                                                               gcTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAAc                                                             >  1:685739/1‑144 (MQ=255)
                                                                                               gcTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAAc                                                             <  2:685739/144‑1 (MQ=255)
                                                                                               gcTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGt                                                       <  1:1161250/150‑1 (MQ=255)
                                                                                                  tACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTggg                                                    <  2:1367493/150‑1 (MQ=255)
                                                                                                   aCTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGc                                                   <  1:1214055/150‑1 (MQ=255)
                                                                                                     tttCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCgg                                                 <  2:515009/150‑1 (MQ=255)
                                                                                                      ttCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGa                                                <  2:1365870/150‑1 (MQ=255)
                                                                                                         aGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTc                                             >  1:509964/1‑150 (MQ=255)
                                                                                                          ggAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCg                                            <  2:293354/150‑1 (MQ=255)
                                                                                                             aCGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCgatg                                         <  1:365061/150‑1 (MQ=255)
                                                                                                             aCGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCgatg                                         <  2:1814398/150‑1 (MQ=255)
                                                                                                              cGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCgatga                                        <  1:1219412/150‑1 (MQ=255)
                                                                                                              cGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCgatga                                        <  2:857385/150‑1 (MQ=255)
                                                                                                                  aaCCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAATGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGAc                                    <  1:1758756/150‑1 (MQ=255)
                                                                                                                    ccGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACAt                                  >  1:1441738/1‑150 (MQ=255)
                                                                                                                     cGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATc                                 <  1:910190/150‑1 (MQ=255)
                                                                                                                         tGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGa                             <  1:1353634/150‑1 (MQ=255)
                                                                                                                           aCGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGAGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTTAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGAGGACTAGATGATGACATCAAGaaa                           >  1:1303807/1‑148 (MQ=255)
                                                                                                                               aaCACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCg                       >  2:415850/1‑150 (MQ=255)
                                                                                                                               aaCACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCg                       >  1:1305904/1‑150 (MQ=255)
                                                                                                                                 cacTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTc                     >  1:202501/1‑150 (MQ=255)
                                                                                                                                 cacTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTc                     <  1:880556/150‑1 (MQ=255)
                                                                                                                                  acTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTcc                    >  1:1992829/1‑150 (MQ=255)
                                                                                                                                  acTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTcc                    >  2:785894/1‑150 (MQ=255)
                                                                                                                                    cgcagcgaTGCTGAAGCCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGTTCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAg                  >  1:1867748/8‑150 (MQ=255)
                                                                                                                                        gcgaTGCTGAAGCCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACAGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTc              >  1:1924317/4‑150 (MQ=255)
                                                                                                                                         ccATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCt             >  1:1788348/1‑150 (MQ=255)
                                                                                                                                          gaTGCTGAAGCCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTg            >  1:698929/2‑150 (MQ=255)
                                                                                                                                           aTGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTgg           <  1:894101/150‑1 (MQ=255)
                                                                                                                                            tGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGt          >  2:954068/1‑150 (MQ=255)
                                                                                                                                                gAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCg      >  2:954727/1‑150 (MQ=255)
                                                                                                                                                    gCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGTTTc  >  1:1636441/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGTTTC  >  NC_000913/691242‑691539

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