Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,068,922 G→T 15.2% intergenic (‑114/+25) 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,9220GT15.2% 469.8 / 61.2 170intergenic (‑114/+25)argO/mscSarginine transporter/mechanosensitive channel protein, small conductance
Reads supporting (aligned to +/- strand):  ref base G (80/63);  new base T (15/11);  total (96/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 181 total aligned reads displayed.

ATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGC  >  NC_000913/3068823‑3069021
                                                                                                   |                                                                                                   
aTTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCagag                                                                                                     <  1:1289660/100‑1 (MQ=255)
 ttAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCagaga                                                                                                    <  1:328367/100‑1 (MQ=255)
     gTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAACTAATAAGAGAGAGGCCAGAGAGta                                                                                                 >  1:2966002/1‑98 (MQ=255)
      tATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGtatt                                                                                               >  2:467263/1‑100 (MQ=255)
      tATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGtatt                                                                                               >  2:2682169/1‑100 (MQ=255)
      tATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGtatt                                                                                               <  1:1629992/100‑1 (MQ=255)
         cTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATc                                                                                            >  1:1919927/1‑100 (MQ=255)
         cTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAAAGtattat                                                                                             >  1:3539623/1‑99 (MQ=255)
          tGATGCAGATACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGATAAGCTAATAAGAGAGAGGCCAGAGAGTATTATc                                                                                            <  1:2499742/99‑1 (MQ=255)
          tGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCt                                                                                           >  1:206821/1‑100 (MQ=255)
          tGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATAAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCt                                                                                           >  1:2216075/1‑100 (MQ=255)
           gATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCt                                                                                           >  1:850940/1‑99 (MQ=255)
             tGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGc                                                                                        <  2:1690466/100‑1 (MQ=255)
              gCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGcc                                                                                       <  2:2596419/100‑1 (MQ=255)
               cAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCt                                                                                      >  1:748209/1‑100 (MQ=255)
                agagACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGATAATACTCTCTGGCCTc                                                                                     <  1:3665175/100‑1 (MQ=39)
                agagACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCt                                                                                      <  1:493559/99‑1 (MQ=255)
                 gagaCTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGATAATACTCTCTGGCCTct                                                                                    <  2:1514197/100‑2 (MQ=38)
                  agaCTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCa                                                                                    <  1:606210/99‑1 (MQ=255)
                   gaCTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGc                                                                                  <  2:894355/100‑1 (MQ=255)
                      tataTCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGtt                                                                               <  2:3436777/100‑1 (MQ=255)
                      tataTCGCGGCATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGtt                                                                               <  1:3283131/100‑1 (MQ=255)
                       atatCGCGGGATTGTTATGAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTg                                                                              <  1:3024140/100‑1 (MQ=255)
                        tatCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGATAATACTCTCTGGCCTctctcttat                                                                             <  2:1551739/100‑9 (MQ=25)
                          tCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGAtt                                                                           >  2:3035313/1‑100 (MQ=255)
                          tCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGAtt                                                                           >  1:3010164/1‑100 (MQ=255)
                          tCGCGGGATTGTGATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGATAATACTCTCTGGCCTctctcttatta                                                                           <  2:1751163/100‑11 (MQ=25)
                            gcgGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTAc                                                                         <  2:2893941/100‑1 (MQ=255)
                            gcgGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTAc                                                                         >  1:281478/1‑100 (MQ=255)
                            gcgGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCAAGAGAGTATTATCTGGCCTCAGCGTTGATTAc                                                                         >  1:159630/1‑100 (MQ=255)
                             cgGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAACCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACg                                                                        >  1:759794/1‑100 (MQ=255)
                               ggATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCa                                                                      <  2:1293603/100‑1 (MQ=255)
                                 aTTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGc                                                                    <  1:1177001/100‑1 (MQ=255)
                                      tattaGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGt                                                               >  1:1009573/1‑100 (MQ=255)
                                       attaGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGtc                                                              <  1:1354072/100‑1 (MQ=255)
                                          aGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGtcttc                                                           <  1:145841/100‑1 (MQ=255)
                                          aGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGtcttc                                                           >  1:1600923/1‑100 (MQ=255)
                                          aGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGtcttc                                                           >  1:414723/1‑100 (MQ=255)
                                          aGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGtcttc                                                           >  1:414199/1‑100 (MQ=255)
                                            cGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTcttctt                                                         <  2:1055337/100‑1 (MQ=255)
                                              gAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTc                                                       <  2:597/100‑1 (MQ=255)
                                                 aTTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCAcc                                                    <  1:1222705/100‑1 (MQ=255)
                                                    gATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGc                                                 >  2:2400979/1‑100 (MQ=255)
                                                         tatTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTaaa                                            >  1:100304/1‑100 (MQ=255)
                                                         tatTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTaaa                                            >  2:549038/1‑100 (MQ=255)
                                                         tatTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTaaa                                            >  1:2405854/1‑100 (MQ=255)
                                                         tatTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTaaa                                            >  1:3071699/1‑100 (MQ=255)
                                                             acatGGCAATAAGCAAATCTACTAAGAGAGGGGCCACAGAGTATTTTCTGGCGTCAGCGTTGAGTACGCAGCTTTGTCTTCTTTCGCCAGCTTAAAGTTc                                        <  1:680597/97‑1 (MQ=255)
                                                              aTTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTca                                       <  2:3383094/100‑1 (MQ=255)
                                                              aTTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTca                                       <  2:3348441/100‑1 (MQ=255)
                                                               ttGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTcac                                      <  2:772535/100‑1 (MQ=255)
                                                                tGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTcaca                                     >  2:812627/1‑100 (MQ=255)
                                                                tGGCAATCAGAAAAACTAATAACAGAGAGTCCAGAGAGTATTATCTTGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTcaca                                     <  1:2587127/100‑1 (MQ=255)
                                                                  gCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATc                                   >  1:2848787/1‑100 (MQ=255)
                                                                   cAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGCTGAATCCGCGGATGTGTCTTCCTTCCCCCGCGCAAAGCTCACATcc                                  >  1:3428211/1‑100 (MQ=255)
                                                                      tCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTATTCTTTCACCCGCTTAAAGTTCACACCCa                                 >  1:1751191/1‑98 (MQ=255)
                                                                       cAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTCAAGTTCACATCCAttt                              >  2:1437500/1‑100 (MQ=255)
                                                                       cAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCAttt                              <  2:159633/100‑1 (MQ=255)
                                                                        aGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTg                             >  2:22142/1‑100 (MQ=255)
                                                                        aGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTg                             >  2:3762733/1‑100 (MQ=255)
                                                                         gAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGc                            >  1:3701019/1‑100 (MQ=255)
                                                                         gAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGc                            >  2:3640794/1‑100 (MQ=255)
                                                                          aaaaaCTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGc                            <  1:1129311/99‑1 (MQ=255)
                                                                          aaaaaCTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCg                           >  2:2187050/1‑100 (MQ=255)
                                                                           aaaaCTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCgg                          >  2:119582/1‑100 (MQ=255)
                                                                           aaaaCTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCgg                          <  1:3659048/100‑1 (MQ=255)
                                                                              aCTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTa                       <  1:2400967/100‑1 (MQ=255)
                                                                               cTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCATCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGTTa                       >  1:190753/1‑99 (MQ=255)
                                                                               cTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTAc                      >  1:1286622/1‑100 (MQ=255)
                                                                               cTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTAc                      >  2:2215816/1‑100 (MQ=255)
                                                                               cTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTAc                      >  2:351148/1‑100 (MQ=255)
                                                                                 aataaGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACgg                    <  2:2490806/100‑1 (MQ=255)
                                                                                   taaGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACgg                    >  2:2818028/1‑98 (MQ=255)
                                                                                   taaGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGa                  <  2:6852/100‑1 (MQ=255)
                                                                                     acgctgagGCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAg                >  2:141000/6‑100 (MQ=37)
                                                                                     acgctgagGCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAg                >  2:3119669/6‑100 (MQ=37)
                                                                                     acgctgagGCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAg                >  2:2274209/6‑100 (MQ=37)
                                                                                      gagagagGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGc               >  1:3296446/1‑100 (MQ=255)
                                                                                       gctgagGCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCt              >  2:2733989/4‑100 (MQ=38)
                                                                                       gctgagGCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCt              >  2:641397/4‑100 (MQ=38)
                                                                                       gctgagGCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCt              >  1:3686587/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)
                                                                                        ctgcggCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAGCGTTCCCATCCATTTCCGGGTACGGGCAGCt              >  1:1510910/5‑99 (MQ=25)
                                                                                        ctgagGCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGACGCTg             >  1:1790820/3‑100 (MQ=37)
                                                                                        ctgagGCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTg             >  2:1931988/3‑100 (MQ=38)
                                                                                        ctgagGCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTg             >  2:1598424/3‑100 (MQ=38)
                                                                                        ctgagGCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTg             >  1:519312/3‑100 (MQ=38)
                                                                                           agGCCAGAGAGTATTATCTTGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATa          <  1:3598405/100‑1 (MQ=255)
                                                                                            ggCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATAc         >  1:2096876/1‑100 (MQ=255)
                                                                                            ggCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATAc         >  1:3810961/1‑100 (MQ=255)
                                                                                             gCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATAc         <  2:3701042/99‑1 (MQ=255)
                                                                                              ccAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACcg       >  2:2320190/1‑100 (MQ=255)
                                                                                              ccAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACcg       >  1:2790521/1‑100 (MQ=255)
                                                                                               cAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACTTCCCTTTGCGGGTACGGGAAGc               <  1:3253888/91‑1 (MQ=255)
                                                                                                 gagagTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACcggcg    <  2:748217/100‑1 (MQ=255)
                                                                                                  agagTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACcggcgg   <  2:3000016/100‑1 (MQ=255)
                                                                                                   gagTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCggcggc  <  2:207381/100‑1 (MQ=255)
                                                                                                   gagTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCggcggc  <  1:2615289/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
ATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGC  >  NC_000913/3068823‑3069021

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: