Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,966,720 C→T 100% R102C (CGC→TGC)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,966,7200CT97.6% 524.4 / 0.6 162R102C (CGC→TGC) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base T (76/82);  total (79/83)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 171 total aligned reads displayed.

CAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTG  >  NC_000913/3966623‑3966816
                                                                                                 |                                                                                                
 aGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTtgc                                                                                                >  2:301151/1‑99 (MQ=255)
  gCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTtgccc                                                                                              >  1:3498179/1‑100 (MQ=255)
   cTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTtacccg                                                                                             >  1:887709/1‑100 (MQ=255)
    tCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTtgcccgc                                                                                            >  1:3764829/1‑100 (MQ=255)
         ccTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGaaa                                                                                       >  1:1163812/1‑100 (MQ=255)
         ccTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGaaa                                                                                       >  1:891333/1‑100 (MQ=255)
         ccTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGaaa                                                                                       >  2:3766887/1‑100 (MQ=255)
          cTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAaag                                                                                      <  2:1325497/100‑1 (MQ=255)
          cTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGATGATACCATCTCTGGTCAGATTTGCccgccg                                                                                          >  2:3030296/1‑96 (MQ=255)
             gCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAaagaag                                                                                   <  2:1299980/100‑1 (MQ=255)
             gCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAaagaag                                                                                   >  2:3083887/1‑100 (MQ=255)
             gCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAaagaag                                                                                   <  1:2035832/100‑1 (MQ=255)
                   ccTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAAc                                                                             <  2:253398/100‑1 (MQ=255)
                    cTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACg                                                                            <  2:235569/100‑1 (MQ=255)
                    cTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACg                                                                            >  1:3277062/1‑100 (MQ=255)
                     tgatgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGc                                                                           <  2:493304/100‑1 (MQ=255)
                      gatgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCt                                                                          >  2:27778/1‑100 (MQ=255)
                        tgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTAt                                                                        <  1:301153/100‑1 (MQ=255)
                          aCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTAttt                                                                      <  2:3814672/100‑1 (MQ=255)
                           cATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTAtttt                                                                     >  2:1180582/1‑100 (MQ=255)
                           cATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTAtttt                                                                     >  2:1951078/1‑100 (MQ=255)
                             tCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTgc                                                                   <  1:3807423/100‑1 (MQ=255)
                             tCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTgc                                                                   <  2:2764214/100‑1 (MQ=255)
                              cTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTgcg                                                                  >  2:248201/1‑100 (MQ=255)
                              cTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTgcg                                                                  <  2:1467675/100‑1 (MQ=255)
                              cTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTgc                                                                   <  1:1995265/99‑1 (MQ=255)
                               tACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTgcgc                                                                 <  1:2748667/100‑1 (MQ=255)
                               tACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTgcgc                                                                 <  2:1782979/100‑1 (MQ=255)
                                aCGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgct                                                                >  2:767217/1‑100 (MQ=255)
                                 cGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctg                                                               <  1:2413292/100‑1 (MQ=255)
                                 cGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctg                                                               >  1:205561/1‑100 (MQ=255)
                                  gTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctgc                                                              >  2:3385271/1‑100 (MQ=255)
                                   tttCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctgct                                                             <  2:1477417/100‑1 (MQ=255)
                                   tttCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctgct                                                             >  2:2229410/1‑100 (MQ=255)
                                   tttCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctgct                                                             <  1:1754291/100‑1 (MQ=255)
                                   ttaccccTAGCCAAATCCGCCGTGTCAACCACCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctgc                                                              <  1:1321135/96‑1 (MQ=255)
                                     tCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGa                                                           >  1:2292844/1‑100 (MQ=255)
                                      ccccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaa                                                          >  2:70385/1‑100 (MQ=255)
                                       cccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaaa                                                         >  2:3702205/1‑100 (MQ=255)
                                       cccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaaa                                                         <  1:1768870/100‑1 (MQ=255)
                                       acatAGCCAAATCCGCCGTTTCAACCTCCGCACTGGAGATAAAATCTCTGGTAAAATTTGACCGACGAAAGAAGGAGATCGCTATTTGGCTCTGCTGaaa                                                         >  2:3120503/4‑100 (MQ=255)
                                        cgtAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAg                                                        <  1:2895320/98‑1 (MQ=255)
                                        ccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAg                                                        <  1:1910923/100‑1 (MQ=255)
                                         cTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGt                                                       <  1:3498691/100‑1 (MQ=255)
                                          tAGCCAAATCCGCCGTTTCAACCTCCCCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGtt                                                      <  1:767225/100‑1 (MQ=255)
                                           aGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTa                                                     <  1:3059992/100‑1 (MQ=255)
                                             ccAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAAc                                                   <  1:3784968/100‑1 (MQ=255)
                                             ccAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAAc                                                   <  1:499742/100‑1 (MQ=255)
                                              cAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACg                                                  <  1:2102843/100‑1 (MQ=255)
                                              cAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACg                                                  >  2:3076793/1‑100 (MQ=255)
                                               aaaTCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGTTATTTTGCGCTGCTGAAAGTTAACGa                                                 <  1:2322366/100‑1 (MQ=255)
                                                 aTCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAg                                               <  2:1731939/100‑1 (MQ=255)
                                                   ccgccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGtt                                             <  1:1464067/100‑1 (MQ=255)
                                                   ccgccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGtt                                             <  2:3155318/100‑1 (MQ=255)
                                                   ccgccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGtt                                             <  2:459083/100‑1 (MQ=255)
                                                    ngccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTa                                            >  2:3499599/2‑100 (MQ=255)
                                                    cgccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTa                                            >  1:949196/1‑100 (MQ=255)
                                                    cgccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTa                                            >  1:27166/1‑100 (MQ=255)
                                                    cgccgTTTCAACCTCCGCACTGATGATACCATCTATGGTAAGATTTGCCCGCCGAAAGAAGGTGAACTCTATTTTGCGCTGCTGAAAGTTAACgtcggca                                            >  2:431501/1‑94 (MQ=255)
                                                      ccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAAc                                          <  2:2828015/100‑1 (MQ=255)
                                                        gTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACtt                                        <  2:699609/100‑1 (MQ=255)
                                                         tttCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTc                                       <  1:1131697/100‑1 (MQ=255)
                                                          ttCACCCTCCGCACTGATGATACTTTCTCTGGTAAGTTTTGACCGCCAAAATAAGGTGAACGTTATTTTGCGCTGCTGAAAGTTAACTTCGTTAACtttc                                      <  2:3677670/100‑3 (MQ=255)
                                                          ttCAACCTCCGCACTGGTGATACCATCTCTGGTAATATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCTCTGCTGAAAGTTAACGAAGTTAACTTCg                                      <  2:3199551/100‑1 (MQ=255)
                                                            cAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGAc                                    <  1:1070151/100‑1 (MQ=255)
                                                              aCCTCCGGACTGGTGATACCATCTCTGGTAAGATTTGACCGCCGAAAGAAGGTGAACGCTATTTTGCGCTTCTGAAAGTTAACTTCGTTAACtttcag                                    <  2:225879/98‑5 (MQ=255)
                                                               ccTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACaaa                                 >  2:2241992/1‑100 (MQ=255)
                                                               ccTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACaaa                                 >  2:1932581/1‑100 (MQ=255)
                                                                cTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAc                                >  2:3624485/1‑100 (MQ=255)
                                                                 tCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAcc                               <  1:2839539/100‑1 (MQ=255)
                                                                  ccGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCt                              >  1:600522/1‑100 (MQ=255)
                                                                     cACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaa                           >  1:2761926/1‑100 (MQ=255)
                                                                      aCTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaaa                          <  1:3477837/100‑1 (MQ=255)
                                                                         ggTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGATAACTTCGACAAACCTGCAAAc                        >  2:3091541/1‑99 (MQ=255)
                                                                           tGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGcc                     <  1:1038103/100‑1 (MQ=255)
                                                                            gATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGccc                    >  1:3279432/1‑100 (MQ=255)
                                                                            gATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGcc                     >  2:3535958/1‑99 (MQ=255)
                                                                              tACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGc                  <  2:3060387/100‑1 (MQ=255)
                                                                              tACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGc                  >  1:603313/1‑100 (MQ=255)
                                                                              tACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGc                  >  2:3482178/1‑100 (MQ=255)
                                                                                ccATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaa                <  2:1800113/100‑1 (MQ=255)
                                                                                ccATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaa                >  1:710623/1‑100 (MQ=255)
                                                                                ccATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaa                >  2:2573743/1‑100 (MQ=255)
                                                                                 cATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaac               <  2:2942316/100‑1 (MQ=255)
                                                                                 cATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaac               >  1:3186084/1‑100 (MQ=255)
                                                                                 cATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaac               >  1:2877636/1‑100 (MQ=255)
                                                                                  aTCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaaca              <  2:1188289/100‑1 (MQ=255)
                                                                                     tctGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACaaaa           <  1:1220863/100‑1 (MQ=255)
                                                                                     tctGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACaaa            >  1:2742218/1‑99 (MQ=255)
                                                                                      ctGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAAt          <  1:1585460/100‑1 (MQ=255)
                                                                                        ggTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTGACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc        <  1:2698061/100‑1 (MQ=255)
                                                                                        ggTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc        <  1:1557525/100‑1 (MQ=255)
                                                                                        ggTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc        <  1:1264360/100‑1 (MQ=255)
                                                                                        ggTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc        >  1:1208746/1‑100 (MQ=255)
                                                                                        ggTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc        <  1:1040288/100‑1 (MQ=255)
                                                                                          tAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc        <  2:2452266/98‑1 (MQ=255)
                                                                                          tAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctc      >  2:655765/1‑100 (MQ=255)
                                                                                             gATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCttt   >  1:1696252/1‑100 (MQ=255)
                                                                                             gATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCttt   >  2:2726032/1‑100 (MQ=255)
                                                                                             gATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCttt   >  1:3607685/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTG  >  NC_000913/3966623‑3966816

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: