Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,673,783 G→A 16.1% intergenic (+15/+33) yphA → / ← yphB DoxX family inner membrane protein/mutarotase superfamily protein, YphB family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,673,7830GA16.1% 465.0 / 65.3 186intergenic (+15/+33)yphA/yphBDoxX family inner membrane protein/mutarotase superfamily protein, YphB family
Reads supporting (aligned to +/- strand):  ref base G (86/69);  new base A (14/16);  total (100/86)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

Only 100 of 207 total aligned reads displayed.

TTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGA  >  NC_000913/2673684‑2673880
                                                                                                   |                                                                                                 
 ttCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATgc                                                                                                  >  2:2314209/1‑100 (MQ=255)
   cTGGAAGAATGTCTGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATgcgg                                                                                                >  1:2719023/1‑100 (MQ=255)
   cTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCGCTCGAGCGGCGTTAGGAAAGATGCCGGATgcgg                                                                                                <  1:4135950/100‑1 (MQ=255)
    tGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATgcggc                                                                                               >  1:2376398/1‑100 (MQ=255)
    tGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATgcggc                                                                                               >  2:370234/1‑100 (MQ=255)
    tGGAAAGATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATgcggc                                                                                               <  1:3781221/100‑1 (MQ=255)
     ggAAGATTGGCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACAGGGAAGAGGGAAATTTCCCACGATCGACGATAGGAAACATGCCGGATgcggc                                                                                               >  1:1022347/1‑99 (MQ=255)
     ggAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATgcggcg                                                                                              <  2:678925/100‑1 (MQ=255)
        agaaTGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGCAAGATGCCGGATGCGGAGTGa                                                                                           >  1:687296/1‑100 (MQ=255)
           aTGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACg                                                                                        <  2:1273953/100‑1 (MQ=255)
                aGTATCGCTGTCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTCAGTAAAGATGCCGGATGCGGCGTGAACGCCTTa                                                                                   <  1:2896412/100‑1 (MQ=255)
                aGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTa                                                                                   >  2:3204961/1‑100 (MQ=255)
                aGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTa                                                                                   <  1:2193618/100‑1 (MQ=255)
                  tATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATc                                                                                 >  1:3173217/1‑100 (MQ=255)
                   aTCGCTGGCGCGTTCTTGCTCTTGGCAATTACAGGGCCGGGGGAAATTTCTCTCGATCGGCGTTAGGAAAGATGGCGGATGCTGCTTGAACGCCTTATcc                                                                                <  1:3018546/100‑1 (MQ=255)
                      gCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAGCGCCTTATCCgg                                                                              >  1:4243406/1‑99 (MQ=255)
                       cTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGcn                                                                            <  1:283473/100‑2 (MQ=255)
                       cTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                                            <  2:2803393/100‑1 (MQ=255)
                       cTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                                            >  2:847756/1‑100 (MQ=255)
                       cTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                            <  2:3563861/100‑1 (MQ=25)
                        ttgcgcGTTCTTGATAGTGGCAATTACAGGGCCGGCGGCGATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAATGCGTTCACGCCTCATCCGGCAt                                                                           <  1:1460661/98‑1 (MQ=255)
                         ggCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGCAtc                                                                          <  2:2057726/100‑2 (MQ=25)
                         ggCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGCAtc                                                                          <  1:3286701/100‑2 (MQ=25)
                          gcggGTTCTTGCTTTTGGGAATTATCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGTATCCGGCa                                                                            <  1:426706/97‑1 (MQ=25)
                          gcgcGTTCTTGCTATTGGCAATTACCGTGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGCAtct                                                                         <  2:309215/100‑3 (MQ=21)
                          gcgcGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCAttt                                                                         >  2:2161270/1‑100 (MQ=255)
                           cgcgTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGCAtctt                                                                        >  2:1574746/1‑97 (MQ=25)
                            gcgTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaa                                                                       >  1:148506/1‑100 (MQ=255)
                              gTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGCAtctttcc                                                                     <  2:1978388/100‑7 (MQ=21)
                                tCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaataa                                                                    >  1:1726730/1‑99 (MQ=255)
                                tCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATaaa                                                                   >  1:495802/1‑100 (MQ=255)
                                    gCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTAc                                                               >  1:304871/1‑100 (MQ=255)
                                    gCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTAc                                                               >  2:367927/1‑100 (MQ=255)
                                    gCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTAc                                                               >  2:1897008/1‑100 (MQ=255)
                                    gCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAACATGCCGGATGGGGCGTGAACGCCTTATCCGGCATTTAATAAATTAc                                                               >  1:3682167/1‑100 (MQ=255)
                                     cTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACa                                                              <  2:1618589/100‑1 (MQ=255)
                                      tATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAg                                                             <  2:2459503/100‑1 (MQ=255)
                                       aTTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGc                                                            <  2:2095106/100‑1 (MQ=255)
                                       aTTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGc                                                            >  1:3244817/1‑100 (MQ=255)
                                        ttGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGcc                                                           >  2:1126916/1‑100 (MQ=255)
                                          ngCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCAc                                                         >  1:1266968/2‑100 (MQ=255)
                                          ggCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCAc                                                         >  2:633381/1‑100 (MQ=255)
                                           gCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCAct                                                        <  2:1136124/100‑1 (MQ=255)
                                            cAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCActc                                                       >  1:1346444/1‑100 (MQ=255)
                                               ttACGGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAc                                                    <  2:1942070/100‑1 (MQ=255)
                                                tACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAcc                                                   <  2:48410/100‑1 (MQ=255)
                                                tACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAcc                                                   <  1:1856894/100‑1 (MQ=255)
                                                  ccgggccgggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCg                                                 >  1:3410040/1‑100 (MQ=255)
                                                   cgggccgggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGc                                                >  2:2488953/1‑100 (MQ=255)
                                                    gggccgggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTGAATAAATTACAGCCACTCTACCCGCa                                               >  1:3319811/1‑100 (MQ=255)
                                                       ccgggGGCAATTTCTCTTGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACg                                            <  1:2420010/100‑1 (MQ=255)
                                                         gggggCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGAc                                          >  2:1807673/1‑100 (MQ=255)
                                                         gggggCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGAc                                          <  2:2858488/100‑1 (MQ=255)
                                                          ggggCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGaga                                         >  1:28097/1‑98 (MQ=255)
                                                          ggggCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACa                                         >  1:3862612/1‑100 (MQ=255)
                                                            ggCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACAtc                                       <  1:2594430/100‑1 (MQ=255)
                                                              cAATTTCTCTCTATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACCCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACAtctc                                     <  1:355236/100‑1 (MQ=255)
                                                              cAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACAtctc                                     <  2:4317907/100‑1 (MQ=255)
                                                               aaTTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATctct                                    >  2:3543625/1‑100 (MQ=255)
                                                                aTTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTg                                   <  1:1565181/100‑1 (MQ=255)
                                                                aTTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTg                                   <  2:179391/100‑1 (MQ=255)
                                                                  ttCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGaa                                 <  2:2533062/100‑1 (MQ=255)
                                                                   tctctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAg                                >  1:467001/1‑100 (MQ=255)
                                                                    ctctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAg                                >  1:2485557/1‑99 (MQ=255)
                                                                     tctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgtt                              >  1:2241018/1‑100 (MQ=255)
                                                                     tctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgtt                              >  2:2516255/1‑100 (MQ=255)
                                                                      ctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgttg                             >  2:1733682/1‑100 (MQ=255)
                                                                      ctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgttg                             >  2:4257281/1‑100 (MQ=255)
                                                                      ctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgttg                             >  1:2145603/1‑100 (MQ=255)
                                                                        cGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAgttgtt                           >  2:140441/1‑100 (MQ=255)
                                                                         gATCGGCGTTAGGAAAGATTCCTGATTCGGCGTGAACGACTTATCCGGCATTTCATAAATTACAGCCGCGCTACCCGCAACGACATCTCTGAAGttgttg                          <  1:1303552/100‑1 (MQ=255)
                                                                          aTCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACATCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGa                         <  1:757004/100‑1 (MQ=255)
                                                                          aTCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTCCCCGCAACGACATCTCTGAAGTTGTTGa                         <  1:4060921/100‑1 (MQ=255)
                                                                          aTCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGa                         <  2:1822266/100‑1 (MQ=255)
                                                                           tCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTCCCCGCAACGACATCTCTGAAGTTGTTGAt                        <  1:1978765/100‑1 (MQ=255)
                                                                           tCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGAt                        >  2:2110048/1‑100 (MQ=255)
                                                                             ggCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTc                      <  2:317437/100‑1 (MQ=255)
                                                                              gCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTcc                     <  1:313005/100‑1 (MQ=255)
                                                                                 ttAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCt                  <  2:265300/100‑1 (MQ=255)
                                                                                  tAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTg                 >  1:3499048/1‑100 (MQ=255)
                                                                                  tAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTg                 <  1:3735049/100‑1 (MQ=255)
                                                                                   aGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTgg                >  1:1617868/1‑100 (MQ=255)
                                                                                   aGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTgg                >  1:1513684/1‑100 (MQ=255)
                                                                                    ggAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTgct               >  1:4335010/1‑98 (MQ=255)
                                                                                    ggAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGt               <  2:3363186/100‑1 (MQ=255)
                                                                                    ggAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGt               >  2:687690/1‑100 (MQ=255)
                                                                                     gAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTg              >  1:1149160/1‑100 (MQ=255)
                                                                                     attaaaTGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTg              >  2:3409025/6‑100 (MQ=21)
                                                                                          aTGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAgc         >  1:959852/1‑100 (MQ=255)
                                                                                          aTGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAgc         >  2:2207745/1‑100 (MQ=255)
                                                                                          aTGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAGTAAATTACAGCCACTCTACCCGCAACGCTATCTCTGAAGTTGTTGATTCCCCCGGTGCTAgc         >  1:4203737/1‑100 (MQ=21)
                                                                                          aTGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAgc         >  1:1677203/1‑100 (MQ=25)
                                                                                          aTGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAgc         >  2:3762957/1‑100 (MQ=25)
                                                                                          aTGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAgc         >  2:4017169/1‑100 (MQ=25)
                                                                                           tGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAgcg        >  1:1363017/1‑100 (MQ=255)
                                                                                             ccGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCa      >  1:315853/1‑100 (MQ=37)
                                                                                                gATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATg   <  2:3244835/100‑1 (MQ=255)
                                                                                                gATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATg   <  1:2058450/100‑1 (MQ=255)
                                                                                                 aTGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATga  <  2:3857229/100‑1 (MQ=255)
                                                                                                 aTGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACAGCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATga  <  2:2368077/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
TTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGA  >  NC_000913/2673684‑2673880

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