Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 383,477 G→A 26.8% G153G (GGC→GGT yaiP ← putative family 2 glycosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913383,4770GA26.8% 269.3 / 101.3 157G153G (GGC→GGTyaiPputative family 2 glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (57/57);  new base A (27/15);  total (85/72)
Fisher's exact test for biased strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 165 total aligned reads displayed.

GATAATAAACGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGAACATAG  >  NC_000913/383378‑383576
                                                                                                   |                                                                                                   
gATAATAAACGGTGCGCCACCCAGGAGCTGTTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCcgcg                                                                                                     <  2:648513/100‑1 (MQ=255)
 ataataAACGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCgcgc                                                                                                    <  1:1272310/100‑1 (MQ=255)
    ataaACGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCtt                                                                                                 <  2:284253/100‑1 (MQ=255)
    ataaACGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCtt                                                                                                 <  1:1759758/100‑1 (MQ=255)
     taaACGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCttt                                                                                                <  1:1087598/100‑1 (MQ=255)
     taaACGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCttt                                                                                                >  2:58122/1‑100 (MQ=255)
     taaACGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCttt                                                                                                >  1:654949/1‑100 (MQ=255)
     taaACGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCttt                                                                                                >  1:1318163/1‑100 (MQ=255)
      aaaCGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTc                                                                                               <  1:27741/100‑1 (MQ=255)
       aaCGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCa                                                                                              >  1:3749421/1‑100 (MQ=255)
        aCGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCaa                                                                                             >  2:2650877/1‑100 (MQ=255)
        aCGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCa                                                                                              <  2:1577829/99‑1 (MQ=255)
          ggTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACg                                                                                           <  1:1536920/100‑1 (MQ=255)
           gTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGc                                                                                          <  1:649821/100‑1 (MQ=255)
           gTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGc                                                                                          <  2:1981714/100‑1 (MQ=255)
           gTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGc                                                                                          >  1:2852972/1‑100 (MQ=255)
           gTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGc                                                                                          >  2:1866371/1‑100 (MQ=255)
             gtgcCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAg                                                                                        <  2:2303675/98‑1 (MQ=255)
              cgcCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGt                                                                                       >  1:3156868/1‑100 (MQ=255)
               gcCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTa                                                                                      >  1:197195/1‑100 (MQ=255)
               gcCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTa                                                                                      >  1:3665500/1‑100 (MQ=255)
                ccACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTag                                                                                     >  2:3205671/1‑100 (MQ=255)
                ccACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTag                                                                                     >  2:330921/1‑100 (MQ=255)
                ccACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTag                                                                                     >  1:2388876/1‑100 (MQ=255)
                 cACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTaga                                                                                    >  2:1071651/1‑100 (MQ=255)
                 cACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTaga                                                                                    >  2:2454004/1‑100 (MQ=255)
                  aCCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTagag                                                                                   >  1:1481460/1‑100 (MQ=255)
                    ccAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAggg                                                                                 >  2:3753762/1‑100 (MQ=255)
                     cAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGa                                                                                <  1:3809028/100‑1 (MQ=255)
                      aGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGaa                                                                               >  1:1680862/1‑100 (MQ=255)
                      aGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGaa                                                                               >  1:378450/1‑100 (MQ=255)
                            tgctgTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATgccgc                                                                         <  2:1928710/100‑1 (MQ=255)
                            tgctgTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATgccgc                                                                         <  1:1913424/100‑1 (MQ=255)
                            tgctgTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGCATgccgc                                                                         >  1:219409/1‑100 (MQ=255)
                             gctgTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATgccgcc                                                                        >  1:2181274/1‑100 (MQ=255)
                             gctgTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATgccgc                                                                         <  2:759950/99‑1 (MQ=255)
                                 tAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACg                                                                    <  1:3190567/100‑1 (MQ=255)
                                 tAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACg                                                                    <  1:1230249/100‑1 (MQ=255)
                                  aGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGTATGCCGCCTAccg                                                                   <  1:2782320/100‑3 (MQ=255)
                                   gCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGc                                                                  >  2:1262965/1‑100 (MQ=255)
                                   gCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGc                                                                  >  2:2974352/1‑100 (MQ=255)
                                   gCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGc                                                                  <  1:2502424/100‑1 (MQ=255)
                                    cGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCa                                                                 <  1:3303870/100‑1 (MQ=255)
                                      tGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATc                                                               <  1:1578872/100‑1 (MQ=255)
                                       gcgcTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCg                                                              >  2:319594/1‑100 (MQ=255)
                                        cgcTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCgg                                                             >  1:2268242/1‑100 (MQ=255)
                                        cgcTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCgg                                                             >  2:1491379/1‑100 (MQ=255)
                                          cTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCa                                                           >  1:2079094/1‑100 (MQ=255)
                                           ttCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCAc                                                          <  1:4137996/100‑1 (MQ=255)
                                            tCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCAc                                                          <  2:3433961/99‑1 (MQ=255)
                                               tAACAATCATCGGGAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACcgcg                                                      >  2:2560914/1‑100 (MQ=255)
                                                aaCAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCgcgc                                                     <  1:854988/100‑1 (MQ=255)
                                                aaCAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCgcgc                                                     >  1:3394897/1‑100 (MQ=255)
                                                  cAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTc                                                   >  2:951632/1‑100 (MQ=255)
                                                   aaTCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCa                                                  >  1:295395/1‑100 (MQ=255)
                                                     tcatcGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAAt                                                <  2:3119699/100‑1 (MQ=255)
                                                     tcatcGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAAt                                                >  2:1321821/1‑100 (MQ=255)
                                                        tcGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTc                                             >  1:1189172/1‑100 (MQ=255)
                                                        tcGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTc                                             <  1:2102171/100‑1 (MQ=255)
                                                         cGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCt                                            >  1:485884/1‑100 (MQ=255)
                                                          ggCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTg                                           >  1:3094123/1‑100 (MQ=255)
                                                          ggCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTg                                           <  1:725073/100‑1 (MQ=255)
                                                          ggCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTg                                           <  2:1853011/100‑1 (MQ=255)
                                                          ggCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTg                                           >  1:3694041/1‑100 (MQ=255)
                                                          ggCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTg                                           <  2:2360674/100‑1 (MQ=255)
                                                            cAACGTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGcc                                         <  2:3778797/100‑1 (MQ=255)
                                                              aCTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAg                                       <  1:104401/100‑1 (MQ=255)
                                                              aCTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAg                                       <  2:3622413/100‑1 (MQ=255)
                                                              aCTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAg                                       >  2:1170497/1‑100 (MQ=255)
                                                                tttACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCa                                     <  1:958953/100‑1 (MQ=255)
                                                                  tACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCata                                   <  1:1327326/100‑1 (MQ=255)
                                                                   aCGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTCTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCata                                   >  2:3222774/1‑99 (MQ=255)
                                                                     ggTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAg                                >  1:3994983/1‑100 (MQ=255)
                                                                     ggTCGCGCGGATGTGCGGGAACATACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCAtata                                 <  2:1208052/99‑1 (MQ=255)
                                                                       tCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCAtata                                 <  2:3614383/97‑1 (MQ=255)
                                                                        cgcgcgGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCCCCGCGCTCAATTTCGGCCAGCATATAGccc                             >  1:4067890/1‑100 (MQ=255)
                                                                        cgcgcgGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGccc                             <  2:1680850/100‑1 (MQ=255)
                                                                        cgcgcgGATGGGCGGTAACAGACCCGCTCCTTTCAACGCAGTGGAGGTAATTCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGc                               >  2:3293170/1‑98 (MQ=255)
                                                                          cgcgGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCAt                           >  2:2963467/1‑100 (MQ=255)
                                                                           gcgGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCAtt                          <  1:698727/100‑1 (MQ=255)
                                                                              gATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCg                       <  1:3041579/100‑1 (MQ=255)
                                                                                 gtgCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCt                    <  2:1577791/100‑1 (MQ=255)
                                                                                   gCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTg                  >  2:753062/1‑100 (MQ=255)
                                                                                    cGGTAACAGACCCGCACCTTTCAAAGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGACAGCATATAGCCCATTCCGTCTTGa                 >  2:2741685/1‑100 (MQ=255)
                                                                                     ggTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGAt                >  2:2082532/1‑100 (MQ=255)
                                                                                      gTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATc               <  1:3388527/100‑1 (MQ=255)
                                                                                        aaCAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCgg             >  2:1380254/1‑100 (MQ=255)
                                                                                          cAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCggga           >  2:1966534/1‑99 (MQ=255)
                                                                                          cAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGc           <  1:3953131/100‑1 (MQ=255)
                                                                                          cAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGc           <  2:3720180/100‑1 (MQ=255)
                                                                                          cAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGc           >  2:520802/1‑100 (MQ=255)
                                                                                           aGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCg          >  1:426117/1‑100 (MQ=255)
                                                                                              cccGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCggca       >  1:3938959/1‑98 (MQ=255)
                                                                                              cccGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGaa       <  1:3824800/100‑1 (MQ=255)
                                                                                              cccGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGaa       >  1:3994815/1‑100 (MQ=255)
                                                                                                 gcgcCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATCTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGAACAt    <  2:399216/100‑1 (MQ=255)
                                                                                                 gcgcCTTTCAACGCAGTAGAAGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGAACAt    >  2:2527268/1‑100 (MQ=255)
                                                                                                  cgcCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGAACATa   <  1:2454018/100‑1 (MQ=255)
                                                                                                   gcCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGAACATAg  <  1:3483691/100‑1 (MQ=255)
                                                                                                   gcCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGAACATAg  <  1:951638/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GATAATAAACGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGAACATAG  >  NC_000913/383378‑383576

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