Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 443,023 T→C 17.3% intergenic (+26/+28) thiI → / ← yajL tRNA s(4)U8 sulfurtransferase/oxidative‑stress‑resistance chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913443,0230TC17.3% 381.8 / 71.7 179intergenic (+26/+28)thiI/yajLtRNA s(4)U8 sulfurtransferase/oxidative‑stress‑resistance chaperone
Reads supporting (aligned to +/- strand):  ref base T (68/79);  new base C (18/13);  total (87/92)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.24e-02

Only 100 of 179 total aligned reads displayed.

GGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCC  >  NC_000913/442925‑443121
                                                                                                  |                                                                                                  
 gTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGATCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCa                                                                                                  <  1:3511420/100‑1 (MQ=255)
 gTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCa                                                                                                  <  1:929499/100‑1 (MQ=255)
 gTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCa                                                                                                  <  1:502992/100‑1 (MQ=255)
   gatgaGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAg                                                                                                 >  2:3257347/1‑99 (MQ=255)
      gaGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGc                                                                                             <  2:3563912/100‑1 (MQ=255)
      gaGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGCTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATg                                                                                              <  2:3822248/99‑1 (MQ=255)
       aGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCg                                                                                            >  1:2371112/1‑100 (MQ=255)
       aGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCg                                                                                            <  1:3686629/100‑1 (MQ=255)
         ccGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGACCAGATGCGAc                                                                                          <  2:1218595/100‑1 (MQ=255)
         ccGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGACCAGATGCGAc                                                                                          <  1:1020341/100‑1 (MQ=255)
          cGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCaagatcggaa                                                                                         >  1:3546680/1‑91 (MQ=255)
          cGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGACCAGATGCGACg                                                                                         <  1:617488/100‑1 (MQ=255)
          cGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGACCAGATGCGACg                                                                                         <  1:777299/100‑1 (MQ=255)
            tCTGCAGGCGCTCTGTCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGtt                                                                                       <  1:394008/100‑1 (MQ=255)
               gCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGc                                                                                    >  1:1220266/1‑100 (MQ=255)
                cAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGcc                                                                                   >  2:2570287/1‑100 (MQ=255)
                cAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGcc                                                                                   <  1:1423913/100‑1 (MQ=255)
                cAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGcc                                                                                   <  1:2570306/100‑1 (MQ=255)
                cAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGACCAGATGCGACGtaaaaa                                                                                   <  1:3088267/100‑6 (MQ=38)
                 aGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCg                                                                                  <  1:2086391/100‑1 (MQ=255)
                 aGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCg                                                                                  <  1:2196215/100‑1 (MQ=255)
                 aGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCg                                                                                  <  2:65251/100‑1 (MQ=255)
                   gcgcTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCa                                                                                <  1:1078707/100‑1 (MQ=255)
                      cactATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATcc                                                                             >  1:2571365/3‑100 (MQ=255)
                       tctATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGACGCATCTGGCCAGATGCGACGTTTGCCGCATcc                                                                             >  2:3311357/1‑99 (MQ=255)
                         tATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTCTATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGac                                                                          <  2:1954312/100‑1 (MQ=255)
                         tATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGac                                                                          <  1:486841/100‑1 (MQ=255)
                         tATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGac                                                                          <  1:3686400/100‑1 (MQ=255)
                          aTCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGaca                                                                         >  1:3205377/1‑100 (MQ=255)
                          aTCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGaca                                                                         >  1:3624200/1‑100 (MQ=255)
                          aTCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGGCAGATGCGACGTTTGCCGCATCCGaca                                                                         >  2:1320418/1‑100 (MQ=255)
                           tCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGacac                                                                        <  2:3189989/100‑1 (MQ=255)
                              gcgcGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACactac                                                                     >  2:1989887/1‑100 (MQ=255)
                                gcgAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTc                                                                   >  2:3771071/1‑100 (MQ=255)
                                 cgAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCg                                                                  >  1:1529995/1‑100 (MQ=255)
                                   aGCAGTGCTTTAACAATGTGAAGGTAGATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGta                                                                <  1:492502/100‑1 (MQ=255)
                                   aGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGta                                                                <  2:2701606/100‑1 (MQ=255)
                                    gCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGtaa                                                               >  1:2875661/1‑100 (MQ=255)
                                     cAGGGCTTTAACAATGTGAAGGTATATCGCACGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGtaat                                                              <  2:2162970/100‑1 (MQ=255)
                                       gggCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGtaataa                                                            <  1:277121/100‑1 (MQ=255)
                                        ggCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTaataat                                                           <  1:1414650/100‑1 (MQ=255)
                                         gcgttAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGAGACGTTTGCCGCATCCGCCACTACTCGTaataat                                                           >  2:491979/4‑99 (MQ=255)
                                         gCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCATAATAAtt                                                          >  2:2764365/1‑100 (MQ=255)
                                          cTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATta                                                         >  1:2006779/1‑100 (MQ=255)
                                           tttAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCTACGTTTGCCGCATCCGACACTACTCGTAATAATtat                                                        <  1:3789856/100‑1 (MQ=255)
                                            ttAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATtata                                                       <  2:3378712/100‑1 (MQ=255)
                                            ttAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATtata                                                       <  1:1674306/100‑1 (MQ=255)
                                            ttAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGCAATAATtata                                                       >  1:391298/1‑100 (MQ=255)
                                             tAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATaa                                                      <  1:3459295/100‑1 (MQ=255)
                                              aaCAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATaaa                                                     <  1:3150097/100‑1 (MQ=255)
                                              aaCAATGTGAAGGTATATCGACCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATaaa                                                     >  2:2327667/1‑100 (MQ=255)
                                               aCAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAAt                                                    <  2:488627/100‑1 (MQ=255)
                                                cAATGTGAAGGTATATCGGCCGTAGTGTGTGGTTTTTACGCCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGAAACTACTCGTAATAATTATAAATc                                                   >  2:1514239/1‑100 (MQ=255)
                                                cAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATc                                                   <  1:2300034/100‑1 (MQ=255)
                                                 aaTGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATcc                                                  <  1:264541/100‑1 (MQ=255)
                                                  aTGTGAAGGTATATCGCCCGTAGTTTGTGGTTTTTACGACGAATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAAGTATAAATccc                                                 >  2:1266579/1‑100 (MQ=255)
                                                  aTGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATccc                                                 <  2:956884/100‑1 (MQ=255)
                                                   tgtgAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCt                                                >  1:431128/1‑100 (MQ=255)
                                                    gtgAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTg                                               <  1:2531983/100‑1 (MQ=255)
                                                     tgAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTgc                                              >  1:2490543/1‑100 (MQ=255)
                                                      gAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTgcc                                             >  2:2884476/1‑100 (MQ=255)
                                                            atatCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATc                                       <  1:2814982/100‑1 (MQ=255)
                                                             tatCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCa                                      <  1:1896334/100‑1 (MQ=255)
                                                              atCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCAc                                     <  2:1007885/100‑1 (MQ=255)
                                                              atCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCAc                                     >  2:748548/1‑100 (MQ=255)
                                                              atCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCAc                                     >  1:3364137/1‑100 (MQ=255)
                                                               tCGCCGGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCAcc                                    <  2:3675599/100‑1 (MQ=255)
                                                                 gCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAg                                  <  1:835436/100‑1 (MQ=255)
                                                                    cGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTtg                               <  1:372282/100‑1 (MQ=255)
                                                                    cGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTtg                               >  2:2662041/1‑100 (MQ=255)
                                                                    cGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTtg                               <  1:3641110/100‑1 (MQ=255)
                                                                         tttGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATg                          <  2:3704906/100‑1 (MQ=255)
                                                                         tttGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATg                          <  2:3028374/100‑1 (MQ=255)
                                                                           tgtgGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGcc                        >  1:2434957/1‑100 (MQ=255)
                                                                            gtgGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCa                       >  1:1321109/1‑100 (MQ=255)
                                                                              ggTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACt                     >  1:822769/1‑100 (MQ=255)
                                                                              ggTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACAACTGGTGATggcggc                     <  2:2497079/100‑6 (MQ=255)
                                                                               cggcaaaCGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCAc                      >  2:1401299/7‑98 (MQ=38)
                                                                                tttttACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTc                   >  1:2808691/1‑100 (MQ=255)
                                                                                tttttACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTc                   <  2:1442968/100‑1 (MQ=255)
                                                                                  tttACGTCGCATCAGGTCAGATGCGACGATTGCCGCATCCGACACTACTCGTGATACTTATAACTCCCTGCCACCATCACCAATTGTGATGCCACtt                    >  2:1587297/1‑97 (MQ=255)
                                                                                   ttACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATg                <  1:2996748/100‑1 (MQ=255)
                                                                                    tACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATgg               <  1:3573415/100‑1 (MQ=255)
                                                                                     aCGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATggg              >  2:2883881/1‑100 (MQ=255)
                                                                                     aCGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATggg              >  2:2417939/1‑100 (MQ=255)
                                                                                      cGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATTGGc             >  2:2025771/1‑100 (MQ=255)
                                                                                      cGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATCATTATAAATCCATGACAACATCACCAGTTGTGATGCCACTTCATGGGc             >  2:286642/1‑100 (MQ=255)
                                                                                      cGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGc             >  1:683874/1‑100 (MQ=255)
                                                                                      cGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGc             >  2:541825/1‑100 (MQ=255)
                                                                                      cGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGc             >  2:1453339/1‑100 (MQ=255)
                                                                                      cGTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGc             >  1:3106545/1‑100 (MQ=255)
                                                                                       gTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCt            >  1:505233/1‑100 (MQ=255)
                                                                                       gTCGCATCTGACCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCt            >  1:162399/1‑100 (MQ=255)
                                                                                          gCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCtttt         <  1:2294986/100‑1 (MQ=255)
                                                                                          gCAGCTGGTCAGATGCGACGTGTGCCGCATCCGACACTACTCTTAATAATTATAAAGCCCGGCCGCCATCACCAGTTGTGATGTCACTTCATGGTCt            <  2:3010523/97‑1 (MQ=255)
                                                                                            aTCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTc       >  2:3519797/1‑100 (MQ=255)
                                                                                                ggTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATATTTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGc   <  2:2673899/100‑1 (MQ=255)
                                                                                                ggTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGc   >  2:2461561/1‑100 (MQ=255)
                                                                                                ggTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGc   >  2:1107080/1‑100 (MQ=255)
                                                                                                 gTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGcc  >  2:1618046/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
GGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCATGGGCTTTTTCACGCC  >  NC_000913/442925‑443121

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: