Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,647,4810TC43.0% 16.1 / 110.8 87C395R (TGC→CGC) gorglutathione oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (26/23);  new base C (21/16);  total (48/39)
Fisher's exact test for biased strand distribution p-value = 8.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.86e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCAC  >  NC_000913/3647345‑3647615
                                                                                                                                        |                                                                                                                                      
tGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGc                                                                                                                                      <  2:468490/139‑1 (MQ=255)
  gTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGt                                                                                                                                    <  1:77695/139‑1 (MQ=255)
      tGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGt                                                                                                                                    >  2:693371/1‑135 (MQ=255)
      tGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGt                                                                                                                                    <  1:693371/135‑1 (MQ=255)
       gTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGAt                                                                                                                               <  1:505316/139‑1 (MQ=255)
        ttGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCTGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATc                                                                                                                              <  1:763396/139‑1 (MQ=255)
        ttGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATc                                                                                                                              >  1:191036/1‑139 (MQ=255)
         tGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCCCTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCg                                                                                                                             <  1:713497/139‑1 (MQ=255)
            tttAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGAt                                                                                                                               >  2:93649/1‑134 (MQ=255)
            tttAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGAt                                                                                                                               <  1:93649/134‑1 (MQ=255)
             ttAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGgaag                                                                                                                         >  1:452698/1‑139 (MQ=255)
                aCGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGgaagaga                                                                                                                      >  1:174107/1‑139 (MQ=255)
                aCGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGgaagaga                                                                                                                      <  2:191036/139‑1 (MQ=255)
                     aCCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGAtt                                                                                                                 >  1:335649/1‑139 (MQ=255)
                                      aGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCAt                                                                                                <  2:369708/139‑1 (MQ=255)
                                      aGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCAt                                                                                                <  1:608662/139‑1 (MQ=255)
                                       gTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCAtt                                                                                               <  2:174107/139‑1 (MQ=255)
                                             cgacgaTCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGga                                                                                                                           >  2:213732/1‑105 (MQ=255)
                                             cgacgaTCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGga                                                                                                                           <  1:213732/105‑1 (MQ=255)
                                                 gaTCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTa                                                                                     >  2:536678/1‑139 (MQ=255)
                                                  aTCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGAt                                                                                                                               <  1:669532/96‑1 (MQ=255)
                                                  aTCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATc                                                                                                                              >  2:696574/1‑97 (MQ=255)
                                                  aTCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATc                                                                                                                              <  1:696574/97‑1 (MQ=255)
                                                  aTCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGGGCGTTGGAt                                                                                                                               >  2:669532/1‑96 (MQ=255)
                                                  aTCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTAt                                                                                    >  1:20797/1‑139 (MQ=255)
                                                     aGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGa                                                                                 <  1:536678/139‑1 (MQ=255)
                                                               gtATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTg                                                                       >  1:846529/1‑139 (MQ=255)
                                                                 ataAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCAtt                                                                                               >  1:942495/1‑113 (MQ=255)
                                                                 ataAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCAtt                                                                                               <  2:942495/113‑1 (MQ=255)
                                                                 ataAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGGAAGAGAAGATTGTCGGTATTc                                                                                                       <  2:83811/105‑1 (MQ=255)
                                                                 ataAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGGAAGAGAAGATTGTCGGTATTc                                                                                                       >  1:83811/1‑105 (MQ=255)
                                                                     aTCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGc                                                                 >  2:661773/1‑139 (MQ=255)
                                                                     aTCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGc                                                                 >  2:378123/1‑139 (MQ=255)
                                                                     aTCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGc                                                                 >  2:461374/1‑139 (MQ=255)
                                                                       ccTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCtt                                                               >  1:185180/1‑139 (MQ=255)
                                                                        ctctTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCAc                                                                                                     >  1:301560/1‑100 (MQ=255)
                                                                        ctctTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCAc                                                                                                     <  2:301560/100‑1 (MQ=255)
                                                                          ctTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGAt                                                                                                                  <  1:392884/85‑1 (MQ=255)
                                                                          ctTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGAt                                                                                                                  >  2:392884/1‑85 (MQ=255)
                                                                          ctTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTcgc                                                            >  2:580255/1‑139 (MQ=255)
                                                                          ctTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTcgc                                                            >  2:675110/1‑139 (MQ=255)
                                                                            ttCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCgg                                                          <  2:20797/139‑1 (MQ=255)
                                                                                   cgATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCt                                                   >  2:845202/1‑139 (MQ=255)
                                                                                   cgATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCt                                                   <  1:452525/139‑1 (MQ=255)
                                                                                    gATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTg                                                  <  1:461374/139‑1 (MQ=255)
                                                                                     aTGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGa                                                 >  2:389088/1‑139 (MQ=255)
                                                                                        tataCCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTc                                                                                                       >  2:694035/1‑82 (MQ=255)
                                                                                        tataCCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTc                                                                                                       <  1:694035/82‑1 (MQ=255)
                                                                                        tataCCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGa                                              >  2:820524/1‑139 (MQ=255)
                                                                                             cgccgTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCa                                                                                                      >  1:861759/1‑78 (MQ=255)
                                                                                             cgccgTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCa                                                                                                      <  2:861759/78‑1 (MQ=255)
                                                                                             cgccgTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGa                                                                                 >  1:560460/1‑99 (MQ=255)
                                                                                             cgccgTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGa                                                                                 <  2:560460/99‑1 (MQ=255)
                                                                                               ccgTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTgg                                                       >  1:386704/1‑123 (MQ=255)
                                                                                               ccgTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTgg                                                       <  2:386704/123‑1 (MQ=255)
                                                                                                cgTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCa                                      >  2:95821/1‑139 (MQ=255)
                                                                                                cgTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGACGATGGGGGCa                                      >  1:144209/1‑139 (MQ=255)
                                                                                                cgTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCa                                      >  2:296044/1‑139 (MQ=255)
                                                                                                 gTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCaa                                     <  1:378123/139‑1 (MQ=255)
                                                                                                  tCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGgaagagaaga                                                                                                                   <  1:941786/60‑1 (MQ=255)
                                                                                                  tCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGgaagagaaga                                                                                                                   >  2:941786/1‑60 (MQ=255)
                                                                                                  tCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAAc                                    >  2:334356/1‑139 (MQ=255)
                                                                                                  tCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGGAAGAGAAGATTGACGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAAc                                    >  2:612107/1‑139 (MQ=255)
                                                                                                    accacTCACCGCCAGCCGTGCCGCATGAAGCTGGTGCGCGATGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCa                                  >  2:615284/1‑139 (MQ=255)
                                                                                                       acTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGGTGAAGATGGGGGCAACCaaaa                               >  1:179470/1‑139 (MQ=255)
                                                                                                         tCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGAt                                                                                                                  <  1:939104/54‑1 (MQ=255)
                                                                                                         tCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGAt                                                                                                                  >  2:939104/1‑54 (MQ=255)
                                                                                                         tCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCaaaaaa                             >  2:891212/1‑139 (MQ=255)
                                                                                                            ccGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGAc                          <  2:335649/139‑1 (MQ=255)
                                                                                                             cGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGga                                                                                                                           >  2:981186/1‑41 (MQ=255)
                                                                                                             cGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGga                                                                                                                           <  1:981186/41‑1 (MQ=255)
                                                                                                             cGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCt                                                                                           <  1:153604/73‑1 (MQ=255)
                                                                                                             cGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCt                                                                                           >  2:153604/1‑73 (MQ=255)
                                                                                                             cGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACt                         >  2:534178/1‑139 (MQ=255)
                                                                                                             cGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGga                                                                                                                           <  2:726832/41‑1 (MQ=255)
                                                                                                             cGCCAGCCGTGCCGCATGAAGCTGGTGCGCGTTGGATCGga                                                                                                                           >  1:726832/1‑41 (MQ=255)
                                                                                                                   ccGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCt                                                                                           <  2:242380/67‑1 (MQ=255)
                                                                                                                   ccGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCt                                                                                           >  1:242380/1‑67 (MQ=255)
                                                                                                                   ccGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACa                   <  2:452698/139‑1 (MQ=255)
                                                                                                                   ccGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACa                   <  2:846529/139‑1 (MQ=255)
                                                                                                                      tGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGaaa                                                                            <  1:57472/79‑1 (MQ=255)
                                                                                                                      tGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGaaa                                                                            >  2:57472/1‑79 (MQ=255)
                                                                                                                      tGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATa                >  1:488990/1‑139 (MQ=255)
                                                                                                                       gCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATAc               >  2:746758/1‑139 (MQ=255)
                                                                                                                           cATGAAGCTGGTGCGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCtt                                                               <  1:96404/87‑1 (MQ=255)
                                                                                                                           cATGAAGCTGGTGCGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCtt                                                               >  2:96404/1‑87 (MQ=255)
                                                                                                                                    ggTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCAc  <  1:95821/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                      
TGGTACTGTTGGTTTAACGGAACCGCAGGCGCGCGAGCAGTATGGCGACGATCAGGTGAAAGTGTATAAATCCTCTTTCACCGCGATGTATACCGCCGTCACCACTCACCGCCAGCCGTGCCGCATGAAGCTGGTGTGCGTTGGATCGGAAGAGAAGATTGTCGGTATTCACGGCATTGGCTTTGGTATGGACGAAATGTTGCAGGGCTTCGCGGTGGCGCTGAAGATGGGGGCAACCAAAAAAGACTTCGACAATACCGTCGCCATTCAC  >  NC_000913/3647345‑3647615

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: