Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,966,720 C→T 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,7200CT96.2% 429.6 / 3.5 132R102C (CGC→TGC) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base C (4/1);  new base T (70/57);  total (74/58)
Fisher's exact test for biased strand distribution p-value = 3.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 137 total aligned reads displayed.

CAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTG  >  NC_000913/3966623‑3966816
                                                                                                 |                                                                                                
cAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTtgc                                                                                                <  2:2780054/100‑1 (MQ=255)
 aGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTtgcc                                                                                               <  2:2423045/100‑1 (MQ=255)
   cTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAACTCCGCCGTTTCAACCTCCCCACTGGTGATACCATCTCTGGTAAGATTtgcccg                                                                                             >  1:2460637/1‑100 (MQ=255)
    tCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTtgcccgc                                                                                            >  2:1813611/1‑100 (MQ=255)
     ccTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTtgcccgcc                                                                                           <  1:2040877/100‑1 (MQ=255)
       tACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGa                                                                                         >  2:3000611/1‑100 (MQ=255)
        aCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGaa                                                                                        >  1:2844476/1‑100 (MQ=255)
         ccTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGaaa                                                                                       <  1:2709954/100‑1 (MQ=255)
          cTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAaag                                                                                      <  2:2760880/100‑1 (MQ=255)
           tCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAaaga                                                                                     <  1:2900737/100‑1 (MQ=255)
           tCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAaaga                                                                                     >  2:2369931/1‑100 (MQ=255)
            cGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAaagaa                                                                                    >  1:2575001/1‑100 (MQ=255)
            cGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAaagaa                                                                                    <  1:2345687/100‑1 (MQ=255)
            cGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAaagaa                                                                                    >  1:843869/1‑100 (MQ=255)
              ccGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAgg                                                                                  <  1:587202/100‑1 (MQ=255)
                ggTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTg                                                                                <  2:846752/100‑1 (MQ=255)
                ggTCCTCATGACATCTACGTCTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTg                                                                                <  1:2841510/100‑1 (MQ=255)
                 gTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGa                                                                               >  1:2248821/1‑100 (MQ=255)
                 gTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCAACTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGa                                                                               <  1:1073070/100‑1 (MQ=255)
                    cTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACg                                                                            <  1:2681964/100‑1 (MQ=255)
                     tgatgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGc                                                                           >  2:2557410/1‑100 (MQ=255)
                     tgatgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGc                                                                           >  2:2098719/1‑100 (MQ=255)
                      gatgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCt                                                                          >  1:1164081/1‑100 (MQ=255)
                           cATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGTTTTGCCCGCCGAAAGAAGGTGAACGCTAtttt                                                                     >  1:2009870/1‑100 (MQ=255)
                           cATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTAtttt                                                                     >  2:2242346/1‑100 (MQ=255)
                           cATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTAtttt                                                                     >  2:301983/1‑100 (MQ=255)
                           cATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTAtttt                                                                     <  2:1140020/100‑1 (MQ=255)
                            aTCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTg                                                                    >  1:1736926/1‑100 (MQ=255)
                             ttccgatctCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTgc                                                                   <  2:317671/92‑1 (MQ=255)
                             tCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTgc                                                                   >  2:1719947/1‑100 (MQ=255)
                               tACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTgcgc                                                                 <  2:2594801/100‑1 (MQ=255)
                               tACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTgcgc                                                                 >  2:1702324/1‑100 (MQ=255)
                               tACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTgcgc                                                                 <  1:828795/100‑1 (MQ=255)
                               tACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTgcgc                                                                 <  1:2403262/100‑1 (MQ=255)
                                 cGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctg                                                               >  2:2154580/1‑100 (MQ=255)
                                 cGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctg                                                               >  1:1576705/1‑100 (MQ=255)
                                  gTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctgc                                                              <  2:3041184/100‑1 (MQ=255)
                                   tttCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctgct                                                             >  2:1388384/1‑100 (MQ=255)
                                      ccccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTgcagatcggaa                                                          >  1:317672/1‑91 (MQ=255)
                                      ccccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaa                                                          >  2:3150153/1‑100 (MQ=255)
                                           aGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTa                                                     <  2:1359028/100‑1 (MQ=255)
                                            gCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTaa                                                    >  1:144692/1‑100 (MQ=255)
                                            gCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTaa                                                    >  2:1996399/1‑100 (MQ=255)
                                             ccAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAAc                                                   <  1:1621947/100‑1 (MQ=255)
                                             ccAAATCCGCCGTTTCAACCGCCGCACTGGTTATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAAc                                                   <  1:1956616/100‑1 (MQ=255)
                                              cgaaTCCGCCGTTACAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACg                                                  >  1:2275082/3‑100 (MQ=255)
                                              cAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACg                                                  <  1:2640216/100‑1 (MQ=255)
                                              cAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGCAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACg                                                  >  1:1984725/1‑100 (MQ=255)
                                                   ccgccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGtt                                             >  2:2055868/1‑100 (MQ=255)
                                                    cgccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTa                                            >  2:2407357/1‑100 (MQ=255)
                                                    cgccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTa                                            <  1:1036020/100‑1 (MQ=255)
                                                     gccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTaa                                           >  2:1044850/1‑100 (MQ=255)
                                                     gccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGAGTTGCCCGCCGAAAGAAGGTGCACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTaa                                           >  2:336330/1‑100 (MQ=255)
                                                      ccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAAc                                          <  2:2244014/100‑1 (MQ=255)
                                                      ccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGCAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAAc                                          >  1:66735/1‑100 (MQ=255)
                                                        gTTTCAAGCTCCGCACTGGTGATATCATCTCTTGTGAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAActtcgtta                                            <  2:2690641/96‑8 (MQ=255)
                                                        gTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAACGTTCACTTCGTTAACtt                                        <  2:151804/100‑1 (MQ=255)
                                                        ccgatctccTCCGCACTGGTGATACCATCTCTGGTAAGATTTGTCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACTTCGTTAACtt                                        <  1:2752443/93‑1 (MQ=255)
                                                          ttCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCg                                      <  1:2843562/100‑1 (MQ=255)
                                                          ttCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCg                                      >  2:2516395/1‑100 (MQ=255)
                                                          ttCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCg                                      <  1:301984/100‑1 (MQ=255)
                                                           tgaaCCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGACGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACTTCGTTAACtttca                                     <  2:160105/98‑4 (MQ=255)
                                                            cAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGACAGTTAACGAAGTTAACTTCGAc                                    >  2:532054/1‑100 (MQ=255)
                                                            cAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGAc                                    >  1:1195214/1‑100 (MQ=255)
                                                             aaCCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACa                                   >  1:1486838/1‑100 (MQ=255)
                                                             aaCCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACa                                   >  1:506663/1‑100 (MQ=255)
                                                               ccTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAACGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACaaa                                 >  1:1938237/1‑100 (MQ=255)
                                                               ccTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACTTCGTTAACtttcagcag                                 <  1:1601754/100‑8 (MQ=38)
                                                               ccTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACaaa                                 >  2:251370/1‑100 (MQ=255)
                                                               ccTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACaaa                                 >  2:2128713/1‑100 (MQ=255)
                                                                cTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAc                                >  2:1360024/1‑100 (MQ=255)
                                                                cTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAc                                >  2:836504/1‑100 (MQ=255)
                                                                 tCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAcc                               <  1:2453070/100‑1 (MQ=255)
                                                                  ccGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCt                              >  1:548595/1‑100 (MQ=255)
                                                                   cGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTg                             >  1:914554/1‑100 (MQ=255)
                                                                   cGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTg                             <  2:495889/100‑1 (MQ=255)
                                                                     cACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaa                           >  1:1008702/1‑100 (MQ=255)
                                                                      aCTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaaa                          <  1:3000008/100‑1 (MQ=255)
                                                                      aCTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaaa                          <  1:1405516/100‑1 (MQ=255)
                                                                         ggTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACg                       >  1:2966004/1‑100 (MQ=255)
                                                                             aTACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCg                   <  1:395119/100‑1 (MQ=255)
                                                                             aTACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCg                   <  2:551125/100‑1 (MQ=255)
                                                                              tACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGc                  <  2:2189663/100‑1 (MQ=255)
                                                                              tACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGc                  <  1:3165537/100‑1 (MQ=255)
                                                                              tACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGc                  <  2:1164078/100‑1 (MQ=255)
                                                                                ccATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaa                >  2:3272921/1‑100 (MQ=255)
                                                                                   tctcTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaacaa             <  2:2749485/100‑1 (MQ=255)
                                                                                   tctcTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaacaa             <  1:380373/100‑1 (MQ=255)
                                                                                     tctGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACaaaa           <  2:2232965/100‑1 (MQ=255)
                                                                                      ctGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAAt          <  1:1320292/100‑1 (MQ=255)
                                                                                      ctGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAAt          <  1:3162699/100‑1 (MQ=255)
                                                                                      ctGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAAt          >  2:1180280/1‑100 (MQ=255)
                                                                                       tGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATc         >  1:2920791/1‑100 (MQ=255)
                                                                                        ggTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCCACAAAATcc        >  1:1996366/1‑100 (MQ=255)
                                                                                        ggTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc        >  1:13252/1‑100 (MQ=255)
                                                                                        ggTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc        <  1:1254257/100‑1 (MQ=255)
                                                                                           aaGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctct     >  1:1933691/1‑100 (MQ=255)
                                                                                            aGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCtt    <  2:3023023/100‑1 (MQ=255)
                                                                                             gATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCttt   >  2:1932996/1‑100 (MQ=255)
                                                                                              aTTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTg  <  2:1426926/100‑1 (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 < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑