Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,781,051 C→G 6.9% intergenic (+34/‑164) lldD → / → trmL L‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,0510CG 6.9% 547.9 / 24.1 158intergenic (+34/‑164)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (75/72);  new base G (3/8);  total (78/80)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.27e-01

Only 100 of 159 total aligned reads displayed.

CGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCC  >  NC_000913/3780905‑3781198
                                                                                                                                                  |                                                                                                                                                    
cGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCa                                                                                                                                                     <  3:173855/148‑1 (MQ=255)
cGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAg                                                                                                                                                    <  2:68925/149‑1 (MQ=255)
 gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCGCCCTAACCCTCTCCCTCAgg                                                                                                                                                   <  2:171368/149‑1 (MQ=255)
 gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCGCCCTAACCCTCTCCCTCAgg                                                                                                                                                   >  1:171368/1‑149 (MQ=255)
 gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAgg                                                                                                                                                   <  2:185107/149‑1 (MQ=255)
      cTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAGGgaga                                                                                                                                               <  4:48629/148‑1 (MQ=255)
       tGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAGGGAGAgg                                                                                                                                             <  2:74384/149‑1 (MQ=255)
       tGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGagag                                                                                                                                              <  3:355672/148‑1 (MQ=255)
        ggCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAggg                                                                                                                                            <  2:359677/149‑1 (MQ=255)
           gcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGa                                                                                                                                          <  3:117412/148‑1 (MQ=255)
             gAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAgg                                                                                                                                                   >  1:112899/1‑137 (MQ=255)
             gAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAgg                                                                                                                                                   <  2:112899/137‑1 (MQ=255)
              aaaTCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑g                                                                                                                                      >  3:174840/1‑148 (MQ=255)
                aTCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTc                                                                                                                                   <  2:391871/149‑1 (MQ=255)
                     tCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGgcgc                                                                                                                              >  2:53778/1‑149 (MQ=255)
                      cAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCt                                                                                                                             <  2:267769/149‑1 (MQ=255)
                        gCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAgggg                                                                                                                                           <  1:15508/134‑1 (MQ=255)
                        gCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAgggg                                                                                                                                           >  2:15508/1‑134 (MQ=255)
                        gCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCtgt                                                                                                                           >  1:86235/1‑149 (MQ=255)
                        gCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCtg                                                                                                                            >  4:324879/1‑148 (MQ=255)
                          gAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCtgtat                                                                                                                         <  1:245075/149‑1 (MQ=255)
                          gAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCtgta                                                                                                                          <  3:324879/148‑1 (MQ=255)
                              ttACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCtgtatgta                                                                                                                      <  4:399162/148‑1 (MQ=255)
                                 cGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTcc                                                                                                                  <  2:354170/149‑1 (MQ=255)
                                  gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGgc                                                                                                                                >  1:182901/1‑134 (MQ=255)
                                    aaGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTc                                                                                                               >  2:324965/1‑149 (MQ=255)
                                     aGATTCGCTGGTGAAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTc                                                                                                               <  4:256419/148‑1 (MQ=255)
                                         tCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCActc                                                                                                           >  4:19673/1‑148 (MQ=255)
                                             tGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCCCTCTGaaa                                                                                                      <  2:239412/149‑1 (MQ=255)
                                                  cAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGaca                                                                                                 <  2:195053/149‑1 (MQ=255)
                                                    ggggCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTc                                                                                                                   <  3:349397/129‑1 (MQ=255)
                                                    ggggCAGAGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGacac                                                                                                >  4:246140/1‑148 (MQ=255)
                                                          gggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCtgtat                                                                                                                         <  2:64104/117‑1 (MQ=255)
                                                          gggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCAct                                                                                          >  4:1890/1‑148 (MQ=255)
                                                            gTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTTTGAAACGACACCGCActct                                                                                        <  4:76159/148‑1 (MQ=255)
                                                             tAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAgg                                                                                                                                             <  2:391593/95‑1 (MQ=255)
                                                             tAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttt                                                                                      >  2:168868/1‑149 (MQ=255)
                                                             tAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttt                                                                                      <  2:71864/149‑1 (MQ=255)
                                                             tAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCtt                                                                                       <  4:100716/148‑1 (MQ=255)
                                                             tAAAGAGTTGCCTGCGGCACTGGCGCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCATCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACTCCGCACTCtt                                                                                       <  3:267063/148‑1 (MQ=255)
                                                              aaaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttt                                                                                      >  3:181253/1‑148 (MQ=255)
                                                              aaaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttt                                                                                      <  4:282877/148‑1 (MQ=255)
                                                               aaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCtttt                                                                                     >  4:235887/1‑148 (MQ=255)
                                                                agagTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttttt                                                                                    <  4:352492/148‑1 (MQ=255)
                                                                 gagTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCtttttt                                                                                   <  3:228831/148‑1 (MQ=255)
                                                                  agTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttttttt                                                                                  >  4:292620/1‑148 (MQ=255)
                                                                  agTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtc                                                                                 <  1:36475/149‑1 (MQ=255)
                                                                     tGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCAct                                                                                          >  3:324566/1‑137 (MQ=255)
                                                                        cTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCt                                                                            <  4:23720/148‑1 (MQ=255)
                                                                        cTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTc                                                                           >  1:71056/1‑149 (MQ=255)
                                                                           cGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttt                                                                                      <  1:126655/135‑1 (MQ=255)
                                                                           cGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttt                                                                                      >  2:126655/1‑135 (MQ=255)
                                                                                 tGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTcc                                                                   >  3:400279/1‑148 (MQ=255)
                                                                                 tGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCg                                                                  >  1:381886/1‑149 (MQ=255)
                                                                                   gCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgg                                                                 <  4:201025/148‑1 (MQ=255)
                                                                                     tCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtctc                                                                               <  2:147254/132‑1 (MQ=255)
                                                                                     tCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtctc                                                                               >  1:147254/1‑132 (MQ=255)
                                                                                       ccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgaga                                                            <  1:113312/149‑1 (MQ=255)
                                                                                       ccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgag                                                             <  4:181253/148‑1 (MQ=255)
                                                                                        cATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgaga                                                            >  3:366172/1‑148 (MQ=255)
                                                                                        cATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTGTTCTCCCTCGCCCCTCCGGGgaga                                                            >  4:203934/1‑148 (MQ=255)
                                                                                          tGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACt                                                                                                                    >  3:167654/1‑90 (MQ=255)
                                                                                          tGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACt                                                                                                                    <  4:167654/90‑1 (MQ=255)
                                                                                          tGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAgg                                                          >  4:353209/1‑148 (MQ=255)
                                                                                          tGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAgg                                                          >  4:133711/1‑148 (MQ=255)
                                                                                           ggCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑Gt                                                                                                                                     <  4:102973/72‑1 (MQ=255)
                                                                                           ggCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑Gt                                                                                                                                     >  3:102973/1‑72 (MQ=255)
                                                                                            gCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGcc                                                       <  1:335874/149‑1 (MQ=255)
                                                                                              gAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCtgtat                                                                                                                         <  1:290373/81‑1 (MQ=255)
                                                                                              gAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCtgtat                                                                                                                         >  2:290373/1‑81 (MQ=255)
                                                                                               aaaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtct                                                                                <  1:310546/121‑1 (MQ=255)
                                                                                               aaaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtct                                                                                >  2:310546/1‑121 (MQ=255)
                                                                                               aaaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                                    <  2:354739/149‑1 (MQ=255)
                                                                                               aaaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgg                                                     <  4:145496/148‑1 (MQ=255)
                                                                                                 aGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTATCCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttttt                                                                                    <  1:355822/115‑1 (MQ=255)
                                                                                                 aGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgggg                                                   >  3:87554/1‑148 (MQ=255)
                                                                                                     aaTGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                              <  2:345601/149‑1 (MQ=255)
                                                                                                      aTGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAggg                                             <  1:53778/149‑1 (MQ=255)
                                                                                                           gCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                        <  1:66090/149‑1 (MQ=255)
                                                                                                           gCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTATTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaat                                        <  2:346236/149‑2 (MQ=255)
                                                                                                              tAGTCGTTTGCGCCCCCTCACCCAACCCCTCCCCCTCGGGGAGGGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTGTGAAAGGACCCCGCACTCTTTTTTTCTCCCTCGCCCCGCCGGGGAGAGGGCGGGGGTGAGGGAAAAAGGc                                     <  2:134996/149‑1 (MQ=255)
                                                                                                                 tCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCtt                                                                                       <  1:132645/96‑1 (MQ=255)
                                                                                                                  cGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGc                                  >  4:226070/1‑148 (MQ=255)
                                                                                                                               tCACCCTAACCCTCTCCATCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCg                                                                  >  4:186096/1‑103 (MQ=255)
                                                                                                                                tccccTAACCCCCCCCCCCGGGGGGGGGGGGCCGGTTGGGCGTTGTATCTCCCCCCCCCCTCGGAAACGCCCCCGCACTTTTTTTTCCCCCCCCCCCCCCCGGGGGAGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAAtt                    <  2:114045/147‑1 (MQ=255)
                                                                                                                                cACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtt                    >  3:217699/1‑148 (MQ=255)
                                                                                                                                   ccTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                >  1:211241/1‑149 (MQ=255)
                                                                                                                                    cTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAgg                                                          >  3:197698/1‑106 (MQ=255)
                                                                                                                                    cTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                >  3:96649/1‑148 (MQ=255)
                                                                                                                                    cTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg               >  1:194342/1‑149 (MQ=255)
                                                                                                                                      tccccTCTCCCTGAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTCTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCt             >  1:393151/3‑149 (MQ=255)
                                                                                                                                      aaCCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGTCACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGTGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGc              >  3:356582/1‑148 (MQ=255)
                                                                                                                                        cccTCCCCCCGAGGGAGGGGGCCCC‑GTCCGGCGTGGTTTTTCCCCCCCCCCTCGAACCCGACCCCCCCCTTTTTTTTCCCCCCCCCCCCCGCCGGGGAGGGGGCGGGGGGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTAt           <  2:381886/149‑1 (MQ=255)
                                                                                                                                         ccTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTAt           >  3:182698/1‑148 (MQ=255)
                                                                                                                                         ccTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc          >  2:370831/1‑149 (MQ=255)
                                                                                                                                          ctctCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc          >  4:140655/1‑148 (MQ=255)
                                                                                                                                            ctcCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCt        >  3:367256/1‑148 (MQ=255)
                                                                                                                                            ctcCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCt        >  3:271938/1‑148 (MQ=255)
                                                                                                                                                cTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcccc   >  1:75934/1‑149 (MQ=255)
                                                                                                                                                cTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcccc   >  2:344436/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCC  >  NC_000913/3780905‑3781198

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑