Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,966,720 C→T 54.7% 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,7200CT54.7% 24.7 / 95.8 77R102C (CGC→TGC) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base C (20/15);  new base T (25/17);  total (45/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGA  >  NC_000913/3966624‑3966819
                                                                                                |                                                                                                   
aGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTtgcc                                                                                                  <  1:797118/100‑1 (MQ=255)
aGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcc                                                                                                  >  1:2048046/1‑100 (MQ=255)
 gCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTtgccc                                                                                                 >  1:349999/1‑100 (MQ=255)
 gCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcc                                                                                                  >  1:812529/1‑99 (MQ=255)
 gCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcaccc                                                                                                 >  1:257739/1‑100 (MQ=255)
   tCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcccgc                                                                                               >  2:2045206/1‑100 (MQ=255)
    ccTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTtgcccgcc                                                                                              >  2:1944123/1‑100 (MQ=255)
    ccTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTtgcccgcc                                                                                              >  2:1835716/1‑100 (MQ=255)
    ccTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATCtgcccgcc                                                                                              >  2:1830044/1‑100 (MQ=255)
      tACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGa                                                                                            >  2:473593/1‑100 (MQ=255)
       aCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGaa                                                                                           >  1:2209220/1‑100 (MQ=255)
       aCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGaa                                                                                           <  2:865802/100‑1 (MQ=255)
       aCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGaa                                                                                           >  2:1603468/1‑100 (MQ=255)
        ccTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGaaa                                                                                          <  1:1951166/100‑1 (MQ=255)
        ccTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGaaa                                                                                          <  2:961854/100‑1 (MQ=255)
         cTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAaag                                                                                         >  2:1349874/1‑100 (MQ=255)
            gCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAaagaag                                                                                      <  1:688370/100‑1 (MQ=255)
            gCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAaagaag                                                                                      <  2:2234223/100‑1 (MQ=255)
                gTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGa                                                                                  <  2:2081676/100‑1 (MQ=255)
                gTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGa                                                                                  <  1:1572197/100‑1 (MQ=255)
                  ccTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAAc                                                                                >  2:1049639/1‑100 (MQ=255)
                   cTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACTATCTCTGGTAAGATTTGCCCGCCGTAAGAAGGTGAACg                                                                               <  2:1454814/100‑1 (MQ=255)
                    tgatgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGc                                                                              >  1:871121/1‑100 (MQ=255)
                     gatgaCATCTACGTTTCCCCTAGCCAAATCCGCTGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCt                                                                             >  2:459866/1‑100 (MQ=255)
                     catgaCCTCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACAGGTGATACCATCCCTGATAATATTTTCCCGCCGAAAGAAGGTGAACGCt                                                                             >  2:1639718/2‑100 (MQ=255)
                      atgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGAGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTa                                                                            >  1:2308038/1‑100 (MQ=255)
                          cATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTAtttt                                                                        >  2:1227086/1‑100 (MQ=255)
                           aTCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTg                                                                       <  1:231765/100‑1 (MQ=255)
                            tCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTgc                                                                      <  1:76258/100‑1 (MQ=255)
                            tCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTgc                                                                      >  2:624526/1‑100 (MQ=255)
                             cTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTgcg                                                                     >  2:843455/1‑100 (MQ=255)
                              tACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTgcgc                                                                    <  1:1881900/100‑1 (MQ=255)
                                 gTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctgc                                                                 >  1:2246326/1‑100 (MQ=255)
                                 gTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctgc                                                                 <  1:1407705/100‑1 (MQ=255)
                                     ccccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaa                                                             >  1:1395894/1‑100 (MQ=255)
                                     ccccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaa                                                             >  1:1635189/1‑100 (MQ=255)
                                     ccccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGa                                                              <  2:800593/99‑1 (MQ=255)
                                        cTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGt                                                          <  1:764865/100‑1 (MQ=255)
                                           gCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTaa                                                       >  1:1551767/1‑100 (MQ=255)
                                             cAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACg                                                     >  1:1937883/1‑100 (MQ=255)
                                                aTCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAg                                                  >  1:1998596/1‑100 (MQ=255)
                                                 tCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGt                                                 <  1:145702/100‑1 (MQ=255)
                                                 tCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGt                                                 <  1:696130/100‑1 (MQ=255)
                                                   cgccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGGTTTGCCCGCCGAAAGAAGGTGAACGCTATTTCGCGCTGCTGATAGTTAACGGAGtn                                                <  2:317550/99‑2 (MQ=255)
                                                   cgccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTa                                               <  2:1991246/100‑1 (MQ=255)
                                                   cgccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTa                                               >  1:1636192/1‑100 (MQ=255)
                                                   cgccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTa                                               >  2:322116/1‑100 (MQ=255)
                                                    gccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTaa                                              >  1:427558/1‑100 (MQ=255)
                                                         ttCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACtt                                           >  2:1852330/1‑98 (MQ=255)
                                                         ttCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCg                                         >  1:1774655/1‑100 (MQ=255)
                                                         ttCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAGTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCg                                         >  2:609040/1‑100 (MQ=255)
                                                          tCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGa                                        >  1:936707/1‑100 (MQ=255)
                                                          tCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGa                                        >  2:1727895/1‑100 (MQ=255)
                                                           cAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGAc                                       >  1:1070729/1‑100 (MQ=255)
                                                            aaCCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACa                                      >  2:2351999/1‑100 (MQ=255)
                                                            aaCCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACa                                      >  2:169215/1‑100 (MQ=255)
                                                            aaCCTCCGCACTGGTGAGACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGACCGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACa                                      <  1:1319152/100‑1 (MQ=255)
                                                              ccTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACaaa                                    >  1:1805291/1‑100 (MQ=255)
                                                              ccTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACaaa                                    >  1:2245725/1‑100 (MQ=255)
                                                              ccTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACaac                                    >  1:1823814/1‑99 (MQ=255)
                                                                tCCGCACTGGTGATACCATCTCTGGTGAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAcc                                  >  1:1429081/1‑100 (MQ=255)
                                                                tCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAcc                                  <  1:2181890/100‑1 (MQ=255)
                                                                tCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAcc                                  >  2:1054443/1‑100 (MQ=255)
                                                                tCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAcc                                  <  1:735358/100‑1 (MQ=255)
                                                                 ccGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCt                                 >  2:898586/1‑100 (MQ=255)
                                                                  cGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTg                                >  1:584668/1‑100 (MQ=255)
                                                                      cTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaaaa                            >  2:748658/1‑100 (MQ=255)
                                                                       tGGTGATACCATCGCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaaaa                            <  2:1751944/99‑1 (MQ=255)
                                                                        ggTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAcc                                  <  2:2252991/92‑1 (MQ=255)
                                                                        ggTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACg                          <  2:1780731/100‑1 (MQ=255)
                                                                            aTACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCg                      <  2:1119623/100‑1 (MQ=255)
                                                                             tACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGc                     >  1:1172644/1‑100 (MQ=255)
                                                                              aCCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGca                    <  1:843459/100‑1 (MQ=255)
                                                                               ccAACTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaa                   <  2:287772/100‑1 (MQ=255)
                                                                                   ctctGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACaaa               >  1:2154883/1‑100 (MQ=255)
                                                                                     ctGGTAAGATTTGCCCACCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAAt             <  1:1349648/100‑1 (MQ=255)
                                                                                     ctGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAAt             <  2:257738/100‑1 (MQ=255)
                                                                                      tGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACaaaa              >  1:819377/1‑98 (MQ=255)
                                                                                       ggTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc           <  1:185745/100‑1 (MQ=255)
                                                                                       ggTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc           >  2:979237/1‑100 (MQ=255)
                                                                                       ggTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc           >  1:2024058/1‑100 (MQ=255)
                                                                                        gTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCct          >  1:1327062/1‑100 (MQ=255)
                                                                                        gTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCct          >  2:676976/1‑100 (MQ=255)
                                                                                         tAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctc         <  2:816107/100‑1 (MQ=255)
                                                                                          aaGATTTGTCCGCCGAAAGAAGTTTAACTTTATTTTGGTCTGCTGAAAGTTAACGTAGTTAACTTCGACACACCTGAAAACGCCCGCAACAAAATCctct        <  2:2043095/100‑1 (MQ=255)
                                                                                           aGATTTGCTCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCtt       >  1:552707/1‑100 (MQ=255)
                                                                                                cgcccgcgGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGATGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTgaga  >  2:1170965/1‑100 (MQ=255)
                                                                                                cgcccgccGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTgaga  <  2:1774642/100‑1 (MQ=255)
                                                                                                |                                                                                                   
AGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGA  >  NC_000913/3966624‑3966819

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