Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,068,929 A→C 15.2% intergenic (‑121/+18) argO ← / ← mscS arginine transporter/mechanosensitive channel protein, small conductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,068,9290AC15.2% 439.0 / 73.3 172intergenic (‑121/+18)argO/mscSarginine transporter/mechanosensitive channel protein, small conductance
Reads supporting (aligned to +/- strand):  ref base A (78/67);  new base C (15/11);  total (93/79)
Fisher's exact test for biased strand distribution p-value = 8.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.57e-01

Only 100 of 180 total aligned reads displayed.

CTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAAT  >  NC_000913/3068832‑3069028
                                                                                                 |                                                                                                   
cTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATc                                                                                                   >  1:3057259/1‑100 (MQ=255)
 tGATGCAGATACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGATAAGCTAATAAGAGAGAGGCCAGAGAGTATTATc                                                                                                   <  1:2499742/99‑1 (MQ=255)
 tGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCt                                                                                                  >  1:206821/1‑100 (MQ=255)
 tGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATAAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCt                                                                                                  >  1:2216075/1‑100 (MQ=255)
  gATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTg                                                                                                 >  2:3510675/1‑100 (MQ=255)
   aTGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTgg                                                                                                <  1:648762/100‑1 (MQ=255)
   aTGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTgc                                                                                                <  1:755203/100‑2 (MQ=255)
    tGCAGAGACTATATCGGGGGATTGTTATTAGCGGAAATTGATAATATTAAGTGGCAGTCAGAAAGACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGc                                                                                               <  1:143251/100‑1 (MQ=255)
    tGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTgg                                                                                                >  1:547974/1‑99 (MQ=255)
    tGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGc                                                                                               <  2:221131/100‑1 (MQ=255)
     gCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGATAATACTCTCTGGcc                                                                                              <  1:447131/100‑1 (MQ=255)
     gCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGcc                                                                                              <  2:2596419/100‑1 (MQ=255)
      cAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGATAATACTCTCTGGCCt                                                                                             <  2:1360046/100‑1 (MQ=255)
      cAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCt                                                                                             >  2:1420185/1‑100 (MQ=255)
      cAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCt                                                                                             >  1:748209/1‑100 (MQ=255)
       agagACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGATAATACTCTCTGGcc                                                                                              <  1:1292098/98‑1 (MQ=255)
       agagACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGATAATACTCTCTGGCCTc                                                                                            <  1:3665175/100‑1 (MQ=39)
       agagACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTc                                                                                            >  1:1221944/1‑100 (MQ=255)
        gagaCTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGATAATACTCTCTGGCCTct                                                                                           <  1:3399590/100‑2 (MQ=39)
          gaCTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGATAATACTCTCTGGCCTctct                                                                                         <  2:3043747/100‑4 (MQ=38)
          gaCTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGc                                                                                         >  2:1517565/1‑100 (MQ=255)
          gaCTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGc                                                                                         <  2:894355/100‑1 (MQ=255)
           aCTATATGGGGGGATTGTTCTTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGTCCAGATAATACTCTCTGGcc                                                                                              <  1:2608809/94‑1 (MQ=255)
           aCTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCg                                                                                        <  2:1572373/100‑1 (MQ=255)
             tataTCGGGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTg                                                                                     <  2:2966017/100‑1 (MQ=255)
             tataTCGCGGCATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGtt                                                                                      <  1:3283131/100‑1 (MQ=255)
              atatCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTg                                                                                     >  2:1930220/1‑100 (MQ=255)
                 tCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGAtt                                                                                  >  1:3010164/1‑100 (MQ=255)
                  cgcgGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTa                                                                                 <  2:3659030/100‑1 (MQ=255)
                   gcgGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTAc                                                                                <  2:2893941/100‑1 (MQ=255)
                   gcgGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTAc                                                                                >  1:281478/1‑100 (MQ=255)
                    cgGGATTGTTATTAGCGGAAAATGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATATGGCCTCAGCGTTGATTACg                                                                               >  1:1487186/1‑100 (MQ=255)
                      ggATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCa                                                                             <  2:1293603/100‑1 (MQ=255)
                        aTTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGc                                                                           <  1:1177001/100‑1 (MQ=255)
                           gTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCttt                                                                        >  1:1614625/1‑100 (MQ=255)
                           gTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCttt                                                                        >  2:2214787/1‑100 (MQ=255)
                             tattaGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGt                                                                      <  2:2542280/100‑1 (MQ=255)
                             tattaGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGt                                                                      >  1:1009573/1‑100 (MQ=255)
                                 aGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGtcttc                                                                  >  1:1600923/1‑100 (MQ=255)
                                 aGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGtcttc                                                                  >  1:414723/1‑100 (MQ=255)
                                     gAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTc                                                              <  2:597/100‑1 (MQ=255)
                                        aTTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCAcc                                                           <  1:1517556/100‑1 (MQ=255)
                                            ataataTTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGcn                                                       <  1:3286470/100‑2 (MQ=255)
                                              aataTTAAGTGTCAATCAGAAACACTAATAAGAGAGAGGCCACAGATTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTGCACCCGCt                                                       <  1:369198/98‑1 (MQ=255)
                                                tatTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTaaa                                                   >  1:100304/1‑100 (MQ=255)
                                                tatTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTaaa                                                   >  1:2405854/1‑100 (MQ=255)
                                                tatTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTaaa                                                   >  2:549038/1‑100 (MQ=255)
                                                tatTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTaaa                                                   >  1:3071699/1‑100 (MQ=255)
                                                 atTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAg                                                  >  2:3157306/1‑100 (MQ=255)
                                                 atTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAg                                                  <  2:1441555/100‑1 (MQ=255)
                                                  ttAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAg                                                  <  1:3557487/99‑1 (MQ=255)
                                                   tAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGtt                                                <  1:2063052/100‑1 (MQ=255)
                                                   tAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGtt                                                >  2:1674792/1‑100 (MQ=255)
                                                    aaTTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTc                                               >  2:1735857/1‑100 (MQ=255)
                                                     aTTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTca                                              >  1:364525/1‑100 (MQ=255)
                                                     aTTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTca                                              <  2:3383094/100‑1 (MQ=255)
                                                      ttGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTcac                                             <  2:772535/100‑1 (MQ=255)
                                                      ttGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTcac                                             >  1:346205/1‑100 (MQ=255)
                                                       tGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTcaca                                            >  2:812627/1‑100 (MQ=255)
                                                         gCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATc                                          >  1:2848787/1‑100 (MQ=255)
                                                         gCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATc                                          <  1:907500/100‑1 (MQ=255)
                                                             tCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTATTCTTTCACCCGCTTAAAGTTCACACCCa                                        >  1:1751191/1‑98 (MQ=255)
                                                              cAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCAttt                                     <  2:159633/100‑1 (MQ=255)
                                                               aGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTg                                    >  2:3762733/1‑100 (MQ=255)
                                                               aGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTg                                    >  2:708805/1‑100 (MQ=255)
                                                                gAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGc                                   >  2:3402982/1‑100 (MQ=255)
                                                                gAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGc                                   >  2:1099308/1‑100 (MQ=255)
                                                                 aaaaaCTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGc                                   <  1:1129311/99‑1 (MQ=255)
                                                                   aaaCTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCggg                                <  1:518020/100‑1 (MQ=255)
                                                                     aCTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTa                              <  1:2400967/100‑1 (MQ=255)
                                                                      cTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCATCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGTTa                              >  1:190753/1‑99 (MQ=255)
                                                                      cTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTAc                             >  1:3816711/1‑100 (MQ=255)
                                                                      cTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTAc                             >  2:2215816/1‑100 (MQ=255)
                                                                      cTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTAc                             >  1:1286622/1‑100 (MQ=255)
                                                                       taataaGAGAGAAGCCAGAGAGTATTATCTGGCCTCAGAGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACg                            <  1:846152/100‑1 (MQ=255)
                                                                        aataaGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACgg                           <  2:2490806/100‑1 (MQ=255)
                                                                        aataaGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACgg                           <  1:3787260/100‑1 (MQ=255)
                                                                          taaGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACgg                           >  2:2818028/1‑98 (MQ=255)
                                                                            acgctgagGCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAg                       >  2:2274209/6‑100 (MQ=37)
                                                                              gctgagGCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCt                     >  2:641397/4‑100 (MQ=38)
                                                                              gctgagGCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCt                     >  1:3681000/4‑100 (MQ=38)
                                                                              agagagGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCt                     <  1:2373865/100‑1 (MQ=255)
                                                                              agagagGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCt                     <  1:3209500/100‑1 (MQ=255)
                                                                               gagagGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTg                    >  1:2827431/1‑100 (MQ=255)
                                                                               ctgagGCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGACGCTg                    >  1:1790820/3‑100 (MQ=37)
                                                                               ctgagGCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTg                    >  2:1931988/3‑100 (MQ=38)
                                                                               ctgagGCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTg                    >  2:3565539/3‑100 (MQ=38)
                                                                                 gagGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGAt                  >  2:430495/1‑100 (MQ=255)
                                                                                   ggCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATAc                >  1:2096876/1‑100 (MQ=255)
                                                                                     ccAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACcg              >  1:2790521/1‑100 (MQ=255)
                                                                                     ccAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACcg              >  2:2320190/1‑100 (MQ=255)
                                                                                      cAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACTTCCCTTTGCGGGTACGGGAAGc                      <  1:3253888/91‑1 (MQ=255)
                                                                                        gagagTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACcggcg           <  2:3005761/100‑1 (MQ=255)
                                                                                        gagagTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACcggcg           <  2:748217/100‑1 (MQ=255)
                                                                                         agagTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACcggcgg          <  2:2696658/100‑1 (MQ=255)
                                                                                          gagTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCggcggc         <  2:207381/100‑1 (MQ=255)
                                                                                            gTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCAt       >  2:2456543/1‑100 (MQ=255)
                                                                                                 atCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAAt  <  1:2711869/100‑1 (MQ=255)
                                                                                                 atCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAAt  <  2:2441459/100‑1 (MQ=255)
                                                                                                 atCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAAt  >  2:1824217/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
CTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAAT  >  NC_000913/3068832‑3069028

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: