Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,781,051 C→G 5.1% 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 5.1% 8615.4 / inf 2432intergenic (+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 (1196/1107);  new base G (54/71);  total (1252/1180)
Fisher's exact test for biased strand distribution p-value = 6.58e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.55e-03

Only 100 of 2446 total aligned reads displayed.

GACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCG  >  NC_000913/3780903‑3781199
                                                                                                                                                    |                                                                                                                                                      
 aCGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCa                                                                                                                                                       <  2:6984293/149‑1 (MQ=255)
  cGCTGACTGGTGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAg                                                                                                                                                      <  2:4768002/149‑1 (MQ=255)
  cGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAg                                                                                                                                                      <  2:2008178/149‑1 (MQ=255)
   gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGTGCCCCCTCACCCTAACCCTCTCCCTCAgg                                                                                                                                                     <  4:5303782/149‑1 (MQ=255)
       aCTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGgaga                                                                                                                                                 <  4:5069399/149‑1 (MQ=255)
        cTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGagag                                                                                                                                                <  3:3476179/149‑1 (MQ=255)
        cTGGCGCGAAATCGATCAGCGAAATTACGCAAGATACGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGagag                                                                                                                                                <  4:5786147/149‑1 (MQ=255)
         tGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAGGGAGAgg                                                                                                                                               <  4:4835316/149‑1 (MQ=255)
          ggCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAGGGAGAggg                                                                                                                                              <  3:4623241/149‑1 (MQ=255)
          ggCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAggg                                                                                                                                              <  4:4604551/149‑1 (MQ=255)
          ggCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAggg                                                                                                                                              <  1:80204/149‑1 (MQ=255)
          ggCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAggg                                                                                                                                              <  2:4011303/149‑1 (MQ=255)
                aaaTCGATCAGCGAAATTACGCGAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGt                                                                                                                                        >  1:4418374/1‑149 (MQ=255)
                 aaTCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGtt                                                                                                                                       >  3:1686302/1‑149 (MQ=255)
                     gATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAgg                                                                                                                                                     <  1:1989072/131‑1 (MQ=255)
                      aTCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcg                                                                                                                                  >  3:6267740/1‑149 (MQ=255)
                       tCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcgc                                                                                                                                 <  4:3636489/149‑1 (MQ=255)
                        cAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCt                                                                                                                                <  1:2719589/149‑1 (MQ=255)
                        cAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCt                                                                                                                                <  3:5826096/149‑1 (MQ=255)
                         aGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCGCCCTCAGGGAGAGGTGACCGTTCGGCGAtg                                                                                                                               <  2:2811325/149‑1 (MQ=255)
                          gCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgt                                                                                                                              >  3:6515756/1‑149 (MQ=255)
                          gCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtg                                                                                                                               <  3:4827921/148‑1 (MQ=255)
                           cGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgta                                                                                                                             >  1:5049504/1‑149 (MQ=255)
                           cGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgta                                                                                                                             <  1:3661989/149‑1 (MQ=255)
                            gAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtat                                                                                                                            <  2:4374367/149‑1 (MQ=255)
                            gAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtat                                                                                                                            <  3:6990016/149‑1 (MQ=255)
                                 tACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACt                                                                                                                       <  3:3248705/149‑1 (MQ=255)
                                   cGCAAGATTCGCTGGTGCAGGGGGTGGGTAAAGAGTTGCCTGCGGCACCGGCCCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTcc                                                                                                                     <  2:3030575/149‑1 (MQ=255)
                                    gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTccc                                                                                                                    >  2:3734771/1‑149 (MQ=255)
                                     cAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCt                                                                                                                   <  3:6026783/149‑1 (MQ=255)
                                     cAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCt                                                                                                                   <  2:6721207/149‑1 (MQ=255)
                                      aaGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAGGGAGAgggttagg                                                                                                                                         >  3:5299192/1‑121 (MQ=255)
                                      aaGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGagag                                                                                                                                                >  1:2956016/1‑119 (MQ=255)
                                       aGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAgg                                                                                                                                               <  4:380549/119‑1 (MQ=255)
                                       aGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCa                                                                                                                 <  2:5026899/149‑1 (MQ=255)
                                       aGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCa                                                                                                                 <  1:316765/149‑1 (MQ=255)
                                        gATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCAc                                                                                                                >  4:1512322/1‑149 (MQ=255)
                                         aTTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCAct                                                                                                               <  4:1866614/149‑1 (MQ=255)
                                             gCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAgg                                                                                                                                                     >  3:5122895/1‑107 (MQ=255)
                                             gCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGa                                                                                                          <  1:621274/149‑1 (MQ=255)
                                               tGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGaaa                                                                                                        <  1:5750539/149‑1 (MQ=255)
                                                ggTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAAc                                                                                                       <  2:2552304/149‑1 (MQ=255)
                                                ggTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAAc                                                                                                       <  4:2366693/149‑1 (MQ=255)
                                                ggTGCAGGGGCTGGGTAAAGAGTTGCCTGCGCCACTGGCTCTCATGGCGAAAGGGAATGCGGAATAGTCGCTTGCGCCCCCTCACCATAACCCTCTCCCTCAGGGAGAGCGGACCGTGCGGCGCTGTATGTACTCCCTCACTC‑TGAAAc                                                                                                       <  3:4401921/149‑1 (MQ=255)
                                                   gCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGac                                                                                                    <  3:1232824/149‑1 (MQ=255)
                                                     aGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGacac                                                                                                  >  4:3341621/1‑149 (MQ=255)
                                                          cTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCAc                                                                                             >  2:535615/1‑149 (MQ=255)
                                                            gggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCActc                                                                                           >  3:4088483/1‑149 (MQ=255)
                                                             ggTAAAGAGTTGCCTGCGGCACTGGCTCCCGTGGCGAAAGGGAATGCGGCATAGTAGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCActct                                                                                          >  2:1709416/1‑149 (MQ=255)
                                                              gTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCtt                                                                                         <  1:6817033/149‑1 (MQ=255)
                                                               tAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCttt                                                                                        >  1:1022653/1‑149 (MQ=255)
                                                                 aaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCttttt                                                                                      >  1:5308048/1‑149 (MQ=255)
                                                                  agagTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCtttttt                                                                                     >  4:2617454/1‑149 (MQ=255)
                                                                    agTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTtc                                                                                   <  1:7082105/149‑1 (MQ=255)
                                                                        gCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGAc                                                                                                                                           <  4:4232694/90‑1 (MQ=255)
                                                                        gCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCg                                                                                                                                     >  3:229913/1‑96 (MQ=255)
                                                                         ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGacac                                                                                                  <  4:4037676/129‑1 (MQ=255)
                                                                          cTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTc                                                                             >  1:1507842/1‑149 (MQ=255)
                                                                                  cTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAAc                                                                                                       >  2:1124675/1‑115 (MQ=255)
                                                                                     gCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑cccc                                                                       >  3:910247/1‑143 (MQ=255)
                                                                                     gCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCggg                                                                 <  1:6003486/149‑1 (MQ=255)
                                                                                     gCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCggg                                                                 <  4:3038926/149‑1 (MQ=255)
                                                                                     gCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTGTTCTGCCTCG‑CCCCTCCggg                                                                 <  2:3268221/149‑1 (MQ=255)
                                                                                       tCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTcc                                                                                <  1:3389269/133‑1 (MQ=255)
                                                                                        cccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTC‑TTTTTTCTCCCTCG‑CCCCTCCGGGgaga                                                             <  4:6780174/149‑1 (MQ=255)
                                                                                         ccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGgaga                                                             <  2:3539928/149‑1 (MQ=255)
                                                                                         ccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGgaga                                                             <  4:5511810/149‑1 (MQ=255)
                                                                                         ccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGgaga                                                             >  2:4203499/1‑149 (MQ=255)
                                                                                         ccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGgaga                                                             >  3:2785146/1‑149 (MQ=255)
                                                                                         ccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCCCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCGCTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGgaga                                                             <  1:4763068/149‑1 (MQ=255)
                                                                                          cATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGGagag                                                            <  2:7204290/149‑1 (MQ=255)
                                                                                           aTGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCttttttt                                                                                    <  3:376707/125‑1 (MQ=255)
                                                                                                gAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgt                                                                                                                              >  4:6285375/1‑79 (MQ=255)
                                                                                                   aGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGGAGAGGGCCGGGGt                                                   <  1:1737988/149‑1 (MQ=255)
                                                                                                       aaTGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGGAGAGGGCCGGGGTGAgg                                               <  2:968533/149‑1 (MQ=255)
                                                                                                       aaTGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGGAGAGGGCCGGGGTGAgg                                               <  3:2572650/149‑1 (MQ=255)
                                                                                                       aaTGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGGAGAGGGCCGGGGTGAgg                                               <  3:5518825/149‑1 (MQ=255)
                                                                                                       aaTGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGGAGAGGGCCGGGGTGAgg                                               <  1:2492049/149‑1 (MQ=255)
                                                                                                           cGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaa                                           >  2:6198613/1‑149 (MQ=255)
                                                                                                             gCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                         <  2:1475682/149‑1 (MQ=255)
                                                                                                              cATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTGCTCCCTCACTC‑TGAAACGACACCGCACTCttttt                                                                                      >  4:6008266/1‑104 (MQ=255)
                                                                                                              cATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAg                                        <  3:6721920/149‑1 (MQ=255)
                                                                                                                 aGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTc                                                                             >  3:203627/1‑110 (MQ=255)
                                                                                                                  gTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCg                                    >  4:2167158/1‑149 (MQ=255)
                                                                                                                   tCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGc                                   >  1:4083524/1‑149 (MQ=255)
                                                                                                                          cgcCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTc                            >  2:5613198/1‑149 (MQ=255)
                                                                                                                           gcCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGGAGAGGGCCGGGGTGAgg                                               <  1:6487245/129‑1 (MQ=255)
                                                                                                                             ccccTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                         >  1:550126/1‑133 (MQ=255)
                                                                                                                                cTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAt                      >  1:1281527/1‑149 (MQ=255)
                                                                                                                                 tCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGga                                                               <  3:300425/107‑1 (MQ=255)
                                                                                                                                   tacatatcccTCTCCCCCGGGGAGAGGGGCCCGTCCGGCGCTGT‑TGGCCTCCCCCCCCCGTGGACCGACCCCGCGC‑CTTTTTTTCGCCCCCGCCCCCCCCGGGGAGGGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtttt                   <  2:6426477/142‑1 (MQ=255)
                                                                                                                                   aCCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGGagag                                                            >  1:3860669/1‑108 (MQ=255)
                                                                                                                                       tAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGc               >  2:5997934/1‑149 (MQ=255)
                                                                                                                                          cccTCTCCCTGAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTAt            >  1:2191177/1‑149 (MQ=255)
                                                                                                                                           ccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGgaga                                                             >  4:117498/1‑99 (MQ=255)
                                                                                                                                           ccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAgg                                       <  4:2802568/121‑1 (MQ=255)
                                                                                                                                              ctcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGga                                                               <  4:716022/94‑1 (MQ=255)
                                                                                                                                                  cTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcccc    >  4:6562101/1‑149 (MQ=255)
                                                                                                                                                   tCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGGAGAGGGCCgggg                                                    <  2:2136596/100‑1 (MQ=255)
                                                                                                                                                   tCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccccc   >  3:1465400/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCG‑CCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCG  >  NC_000913/3780903‑3781199

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