Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,967,461 A→C I349L (ATC→CTC)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,967,4610AC88.9% 340.9 / 27.1 127I349L (ATC→CTC) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base A (9/5);  new base C (59/53);  total (68/59)
Fisher's exact test for biased strand distribution p-value = 5.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 127 total aligned reads displayed.

GCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGT  >  NC_000913/3967365‑3967557
                                                                                                 |                                                                                                
gcgACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCg                                                                                               >  1:1431023/1‑100 (MQ=255)
 agACGGCGCTTATCGATACCGAGACTAACAGAAACGAAGTCATCTCCGATG‑GTTTAACAGGTACACGCAACATGGAACGGCACCTCTCTCGTAAGCTCGc                                                                                              <  2:1951751/99‑1 (MQ=255)
      gcgcTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGaaa                                                                                         <  2:3826006/100‑1 (MQ=255)
       cgcTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGaaaa                                                                                        <  2:3166222/100‑1 (MQ=255)
         cTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAAc                                                                                      <  1:567297/100‑1 (MQ=255)
         attATCGATACCGGCGCTAACATGGACGAAGTTATCTACGAAGAGCTTAA‑AGGTACAGGCACAATGAAACTACACCTCTCTCGTAAGATCGCTGaacaac                                                                                      >  2:1306463/2‑100 (MQ=255)
          ttATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAAcg                                                                                     <  1:3508812/100‑1 (MQ=255)
          ttATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCCGAAAAAcg                                                                                     >  2:1683897/1‑100 (MQ=255)
            aTCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAAcgcg                                                                                   <  1:3150714/100‑1 (MQ=255)
             tCGATACGGCTTCAAAAATGGACGAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCGACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGt                                                                                  >  2:3296359/1‑100 (MQ=255)
              cGATACCGGTTCTATAATGAACGAAGTTATCTACTAAGAATTTAA‑AGGTACAGGCAACAAAGAACTGCACCTCTCTCGTATGACCGCTGAAAAACGCGTc                                                                                 <  2:810799/100‑1 (MQ=255)
               gATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCt                                                                                >  2:3409228/1‑100 (MQ=255)
               gATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCt                                                                                <  1:168826/100‑1 (MQ=255)
                aTACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACCTGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCtt                                                                               <  1:2100555/100‑1 (MQ=255)
                  aCCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTcc                                                                             >  2:3236746/1‑100 (MQ=255)
                    cGGGTCTAAAAAGGACGAAGTGATCTACGAAGACGATAA‑AGTTAAAGGCAACAGGTAACTGCACCTCTCGTTTAAGCTCGCTGAAAAACGCGTCTTCCCg                                                                           <  2:3909478/100‑1 (MQ=255)
                     ggTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCgg                                                                          <  1:2675322/100‑1 (MQ=255)
                     ggTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCgg                                                                          <  2:3317283/100‑1 (MQ=255)
                      gTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGc                                                                         <  1:99507/100‑1 (MQ=255)
                        tCTAAAATGGACGAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTa                                                                       >  1:1246691/1‑100 (MQ=255)
                        tCTAAAATGGACGAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTa                                                                       >  2:3926188/1‑100 (MQ=255)
                         cTAAAATGGACGAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTAt                                                                      >  1:3446888/1‑100 (MQ=255)
                           aaaaGGGACGAAGTTATCTACGAAGAGTTTAA‑AGCTACTGGCAACATGAAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCGTCCCGGCTATCg                                                                    <  2:3854642/100‑1 (MQ=255)
                            aaaTGGACGAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGa                                                                   >  1:1886676/1‑100 (MQ=255)
                            aaaTGGACGAAGTTATCGACGAAGAGTTTAA‑AGGTACAGGCAACATGTAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGa                                                                   >  2:1621872/1‑100 (MQ=255)
                            aaaTGGACGAAGTTATCGACGAAGAGTATAA‑AGGTACAGGCAACATGGAACTGCACCTCCCTCGTAAGCTCGCTGGAAAAGGGGTCTTCCCGGCTATCGa                                                                   >  2:1932617/1‑100 (MQ=255)
                              aTGGACGAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACt                                                                 <  1:3296760/100‑1 (MQ=255)
                                ggACGAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTAc                                                               >  1:417972/1‑100 (MQ=255)
                                 gACGAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACa                                                              <  1:2899057/100‑1 (MQ=255)
                                 gACGAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACa                                                              >  2:2492080/1‑100 (MQ=255)
                                  aCGAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACaa                                                             >  1:2999265/1‑100 (MQ=255)
                                    gAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAAcc                                                           >  1:3829951/1‑100 (MQ=255)
                                    gAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAAcc                                                           <  2:3647541/100‑1 (MQ=255)
                                    gAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAAcc                                                           >  1:1930226/1‑100 (MQ=255)
                                    gAAGGTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTGCCCGGCTATCGACTACAAcc                                                           <  2:941249/100‑1 (MQ=255)
                                         tATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCt                                                      <  1:1881221/100‑1 (MQ=255)
                                          aTCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTg                                                     >  1:2850395/1‑100 (MQ=255)
                                          aTCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTg                                                     <  2:1131900/100‑1 (MQ=255)
                                              aCGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCCGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTAc                                                 >  1:167283/1‑100 (MQ=255)
                                                gAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTAccc                                               <  1:2530980/100‑1 (MQ=255)
                                                gAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTAccc                                               >  2:1198819/1‑100 (MQ=255)
                                                 aaGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCg                                              >  1:3573682/1‑100 (MQ=255)
                                                   gagTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTa                                            <  1:3481743/100‑1 (MQ=255)
                                                     gTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTaaa                                          >  2:1332829/1‑100 (MQ=255)
                                                     gTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTaaa                                          >  2:1671547/1‑100 (MQ=255)
                                                      tttAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAaag                                         >  2:1835089/1‑100 (MQ=255)
                                                       ttAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAaaga                                        <  1:2674254/100‑1 (MQ=255)
                                                        tAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAaagaa                                       <  2:309755/100‑1 (MQ=255)
                                                            aGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAagag                                    >  2:2576382/1‑100 (MQ=255)
                                                             ggTACAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCACGGCTAGCGACTACAACCGTTCTGGTAGCCGTAAAGAAGAgc                                   >  2:3264177/1‑100 (MQ=255)
                                                                 cAGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAgctgct                               <  2:1405837/100‑1 (MQ=255)
                                                                 cAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAgctgct                               >  2:483692/1‑100 (MQ=255)
                                                                  aGGCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTc                              <  1:1102501/100‑1 (MQ=255)
                                                                   ggCAACATGGAACTGCGCCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTCCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCa                             >  1:736532/1‑100 (MQ=255)
                                                                   ggCAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCa                             >  1:3511745/1‑100 (MQ=255)
                                                                    tcAACATGGAACTGCACCTCTCTCGTAAGTTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCTTAAAGAAGAGCTGCTCAc                            <  2:1888312/99‑1 (MQ=255)
                                                                     cAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACg                           <  2:2948823/100‑1 (MQ=255)
                                                                     cAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACg                           >  1:2088452/1‑100 (MQ=255)
                                                                     cAACATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACg                           <  2:3131938/100‑1 (MQ=255)
                                                                      aaCATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGa                          <  2:1713763/100‑1 (MQ=255)
                                                                      aaCATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGa                          >  2:623883/1‑100 (MQ=255)
                                                                       aCATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGAc                         >  1:2917749/1‑100 (MQ=255)
                                                                       aCATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGAc                         >  1:2917000/1‑100 (MQ=255)
                                                                        cATGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACt                        >  2:2426575/1‑100 (MQ=255)
                                                                          tGGAACTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCa                      >  1:562174/1‑100 (MQ=255)
                                                                          tGGAACTGCACCTCTCTCGAAAAAAAAAAAAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCa                      <  2:3565798/100‑1 (MQ=255)
                                                                            gAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAgg                    <  1:1609244/100‑1 (MQ=255)
                                                                             aaCTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGga                   <  2:705439/100‑1 (MQ=255)
                                                                              aCTGCACCTCTCTCGTAAGCTCGCTGAAAACCGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGgaa                  <  1:900736/100‑1 (MQ=255)
                                                                              aCTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGgaa                  <  2:3358368/100‑1 (MQ=255)
                                                                              aCTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGgaa                  <  2:3359946/100‑1 (MQ=255)
                                                                              aCTGCACCCCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGgaa                  >  1:3570582/1‑100 (MQ=255)
                                                                               cTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGgaag                 <  1:304750/100‑1 (MQ=255)
                                                                               cTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGgaag                 >  2:494889/1‑100 (MQ=255)
                                                                               cTGCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGgaag                 <  1:1911259/100‑1 (MQ=255)
                                                                                tGCACCTCTCTCGTAAGCTCGCTGAATTACGCGTCTTCCCGGTTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGCGCTGCTCACGACTCAGgaaga                <  2:454519/100‑1 (MQ=255)
                                                                                tGCACCTCGCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTTCCCGTAAAAAAGAGCTGCTCACGACTCAGgaaga                >  2:737723/1‑100 (MQ=255)
                                                                                 gCACCTCTCTCGTAAGCTCGCTGAAAAACGCTTCTTCCCGGCTATCGACTACAACCGTTCTGTTACGCGTAAAGAAGAGCTGCTCACGACTCAGggagaa               >  2:577816/1‑100 (MQ=255)
                                                                                 gCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGgaagaa               >  1:2554495/1‑100 (MQ=255)
                                                                                 gCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGgaagaa               <  1:2118022/100‑1 (MQ=255)
                                                                                 gCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCATAAAGAAGAGCTGCTCACGACTCAGgaagaa               >  1:3201333/1‑100 (MQ=255)
                                                                                 gCACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTCCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGgaagaa               >  2:1653854/1‑100 (MQ=255)
                                                                                  cACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAAc              <  2:702228/100‑1 (MQ=255)
                                                                                  cACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAAc              <  2:1022245/100‑1 (MQ=255)
                                                                                  cACCTCTCTCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAAc              >  1:3712436/1‑100 (MQ=255)
                                                                                     ctctctCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGc           >  2:1775217/1‑100 (MQ=255)
                                                                                     ctctctCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGc           <  1:4009131/100‑1 (MQ=255)
                                                                                     ctctctCGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGc           <  1:477649/100‑1 (MQ=255)
                                                                                     atctctcGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGc           <  2:2354098/99‑1 (MQ=255)
                                                                                       ctctcGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAg         >  2:508408/1‑100 (MQ=255)
                                                                                         ctcGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGat       >  1:2390816/1‑99 (MQ=255)
                                                                                         ctcGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGaa       >  2:2005848/1‑100 (MQ=255)
                                                                                         ctcGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGaa       >  2:2125444/1‑100 (MQ=255)
                                                                                          tcGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGaaa      >  2:2399325/1‑100 (MQ=255)
                                                                                           cGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGaaaa     >  1:2326238/1‑100 (MQ=255)
                                                                                           cGTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGaaaa     >  2:836255/1‑100 (MQ=255)
                                                                                            gTAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAAt    <  1:1236978/100‑1 (MQ=255)
                                                                                             tAAGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAAtg   <  2:2624738/100‑1 (MQ=255)
                                                                                             tAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGatcgga   >  2:3259043/1‑95 (MQ=255)
                                                                                              aaGCTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAAtgt  <  1:1537632/100‑1 (MQ=255)
                                                                                                 |                                                                                                
GCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAA‑AGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGT  >  NC_000913/3967365‑3967557

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