Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,781,022 T→G 10.6% intergenic (+5/‑193) 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,0220TG10.6% 1506.8 / 46.4 545intergenic (+5/‑193)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (103/372);  new base G (67/3);  total (170/375)
Fisher's exact test for biased strand distribution p-value = 6.19e-35
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 555 total aligned reads displayed.

TCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCA  >  NC_000913/3780873‑3781167
                                                                                                                                                     |                                                                                                                                                  
 cTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGttt                                                                                                                                                   <  1:5698567/150‑1 (MQ=255)
    aTCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTgc‑g                                                                                                                                               >  1:4296036/1‑150 (MQ=255)
    aTCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTgc‑g                                                                                                                                               >  1:4607779/1‑150 (MQ=255)
     tCGAAAAAGAGATGACAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTgc‑gc                                                                                                                                              <  1:1440315/150‑1 (MQ=255)
       gAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑Gccc                                                                                                                                            >  2:463231/1‑150 (MQ=255)
        aaaaaGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑Gcccc                                                                                                                                           <  1:4918205/150‑1 (MQ=255)
        aaaaaGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑Gcccc                                                                                                                                           >  2:1454164/1‑150 (MQ=255)
         aaaaGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCAGGGCGAAAGGGAATGCGGCATAGTCGTGTGC‑Gccccc                                                                                                                                          >  1:1180473/1‑150 (MQ=255)
           aaGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTTC‑Gcccccac                                                                                                                                        >  2:689113/1‑148 (MQ=255)
           aaGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTc                                                                                                                                        <  2:3564182/150‑1 (MQ=255)
           aaGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTc                                                                                                                                        <  2:2947645/150‑1 (MQ=255)
            agagATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCa                                                                                                                                       <  2:2720848/150‑1 (MQ=255)
             gagaTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCAc                                                                                                                                      <  2:5801052/150‑1 (MQ=255)
              agaTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCAcc                                                                                                                                     <  2:2576468/150‑1 (MQ=255)
                    taGTGGCGTTGACGCTGACTGGCGCGAACTCGATGATCGAAATTCCGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTTCCTGCGGCACTGGCTCCCATGGCGTCAGGGTATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAAcc                                                                                                                               <  1:2769601/149‑1 (MQ=255)
                    aaGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGGCCCCCCCAccaaacc                                                                                                                                >  1:2032094/1‑145 (MQ=255)
                         gCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCtctc                                                                                                                          <  2:2153810/150‑1 (MQ=255)
                         gCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCtctc                                                                                                                          <  2:5037461/150‑1 (MQ=255)
                         gCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCGCCCCCAAAcccacac                                                                                                                          >  2:41369/1‑146 (MQ=255)
                          cGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTcc                                                                                                                         <  2:2813861/150‑1 (MQ=255)
                          cGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGAGTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTcc                                                                                                                         <  2:556556/150‑1 (MQ=255)
                           gATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCCCCCCCAAACCCCCTccc                                                                                                                        >  2:4686695/1‑150 (MQ=255)
                           gATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGC‑GCCCCCCCACCCCAAaccccaacc                                                                                                                        >  2:1320234/1‑144 (MQ=255)
                              gACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCGCCCCCGAACCCACTcccgcc                                                                                                                     >  2:837596/1‑147 (MQ=255)
                               aCGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAg                                                                                                                    <  1:4197820/150‑1 (MQ=255)
                               aCGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAg                                                                                                                    <  1:2558778/150‑1 (MQ=255)
                                cGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAgg                                                                                                                   <  2:4224421/150‑1 (MQ=255)
                                 gCNGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGACTGCGGCATAGTCGTGTGC‑GCCCCCCCCCCCCCACCCCCGcccccccgg                                                                                                                  >  1:1279285/1‑144 (MQ=255)
                                     aCTGGCGCGAATTCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGagag                                                                                                              <  2:434348/150‑1 (MQ=255)
                                     aCTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCGCACCCAAACCCACACCCACAGGGagag                                                                                                              >  1:641635/1‑150 (MQ=255)
                                      cTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTGAGGGAGAgg                                                                                                             <  1:5774436/150‑1 (MQ=255)
                                      cTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTGAGGGAGAgg                                                                                                             <  1:4760798/150‑1 (MQ=255)
                                      cTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAgg                                                                                                             <  1:1593865/150‑1 (MQ=255)
                                       tGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAggg                                                                                                            <  2:2100526/150‑1 (MQ=255)
                                         gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGa                                                                                                          <  1:5338354/150‑1 (MQ=255)
                                         gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGa                                                                                                          <  1:3019099/150‑1 (MQ=255)
                                          cgcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCCTGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTGAGGGAGAgggtta                                                                                                         <  2:5699367/150‑4 (MQ=255)
                                          cgcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTGAGGGAGAgggtta                                                                                                         <  2:230048/150‑4 (MQ=255)
                                           gcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGAcc                                                                                                        <  1:2981260/150‑1 (MQ=255)
                                              aaaTCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCACTCCCACAGGGAGAGGGGGCCggt                                                                                                     >  1:5604553/1‑148 (MQ=255)
                                               aaTCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTc                                                                                                    <  1:4695609/150‑1 (MQ=255)
                                                 tCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCgg                                                                                                  <  1:9864/150‑1 (MQ=255)
                                                   gATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcg                                                                                                <  1:3941500/150‑1 (MQ=255)
                                                   gATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCAAACCCACTCCCCCAGGGCGCGGGGACCCGTCGgcg                                                                                                >  1:3203378/1‑150 (MQ=255)
                                                   gATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcg                                                                                                <  2:3037357/150‑1 (MQ=255)
                                                     tCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCt                                                                                              <  2:3310363/150‑1 (MQ=255)
                                                       aGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCGCACCCGAACCCACACCCACAGGGAGAGGGGACCCTTCGGCGCtgt                                                                                            >  1:2411000/1‑150 (MQ=255)
                                                        gCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGACTGCGGCATAGTCGTTTGC‑GCCCCCGCACCCAAACCCACTCCCGCAGGGGGGGGGGGCCCGTCGGCGCagta                                                                                           >  2:5163366/1‑150 (MQ=255)
                                                         cGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGTCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtat                                                                                          <  1:5411008/150‑1 (MQ=255)
                                                         cGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtat                                                                                          <  2:1681164/150‑1 (MQ=255)
                                                         cGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCCCACCCCAAACCACACCCGCAGGGGGAGGGGGCCCGTGGgcgcggggt                                                                                          >  1:170817/1‑145 (MQ=255)
                                                          gAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtatg                                                                                         <  1:1769816/150‑1 (MQ=255)
                                                           aaaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtatgt                                                                                        <  2:3931530/150‑1 (MQ=255)
                                                            aaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtatgta                                                                                       <  1:2833658/150‑1 (MQ=255)
                                                            aaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtatgta                                                                                       <  2:1141349/150‑1 (MQ=255)
                                                                 cGCAAGATTCGGAGGTGGAGGGACTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGTAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTccc                                                                                  <  2:4438410/150‑1 (MQ=255)
                                                                  gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCCCCCCCAAAACCACACCCAACGGGCGAGGGGGACGTTAGGCGCggtatgtaaaacca                                                                                 >  2:320956/1‑144 (MQ=255)
                                                                  gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCCCCCCCAAAAACCCCCCCCCACGGAGAGGGGCCCGGGCGGCGCGGGATTTACAcccc                                                                                 >  1:2311496/1‑149 (MQ=255)
                                                                  gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGC‑GCCCCCGCACCCGAACCCACTCCCGCCCGGGGAGGGGGCCCGGCGGCGCGGTATGTACAcccg                                                                                 >  1:702352/1‑149 (MQ=255)
                                                                   cAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGG‑GCCCCCACACACTAAACCAAAAACGCAGGGGGTGGGGACCGTTAGGCGCTGTATGTTATcccac                                                                                >  1:3601950/1‑148 (MQ=255)
                                                                   cAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCCCCCCCAAAACCACACCCGCAGGGAGAGGGGGCCCGTCGGCGCGGTGTGTAAACCCTc                                                                                >  1:2605943/1‑150 (MQ=255)
                                                                   cAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCACACCCAAACCCTCAACCACAGGGAGGGGGGCCCGGTCGGCGCGGGATGTTCTaccac                                                                                >  1:3408393/1‑148 (MQ=255)
                                                                       aTTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCActc                                                                            <  1:911521/150‑1 (MQ=255)
                                                                        ttCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCActct                                                                           <  1:4508314/150‑1 (MQ=255)
                                                                            cTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaaa                                                                       <  2:3631560/150‑1 (MQ=255)
                                                                                 gctggggCTGGGTAATGTGTTGCCTGCGGCAGTGGCTCCCTTGGCGAATGGTATTGCGGCATTTTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCTGGTAGTGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGaca                                                                  <  1:4446216/147‑1 (MQ=255)
                                                                                         tGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCActc                                                          <  2:5401303/150‑1 (MQ=255)
                                                                                         tGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCACACCCAAAACCACAACCACACGGAGAGGGGACCGGTAGGCGCAGTATGTTATACCACACAATGAAACCACACCCCAcac                                                          >  1:3505013/1‑148 (MQ=255)
                                                                                         ggggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCCCTCTGAAACGACCCCGCActc                                                          <  2:1523856/149‑1 (MQ=255)
                                                                                           ggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCtt                                                        <  1:5209309/150‑1 (MQ=255)
                                                                                             tAAAGAGTTGCCTGCGGCACTGGCTCCCTTGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCtttt                                                      <  1:5026932/150‑1 (MQ=255)
                                                                                              aaaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttttt                                                     <  2:1180473/150‑1 (MQ=255)
                                                                                               aaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCtttttt                                                    <  2:2255271/150‑1 (MQ=255)
                                                                                                 gagTTGCCGGCGGCACTGGCTCCCAGGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtc                                                  <  2:5281031/150‑1 (MQ=255)
                                                                                                    ttGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCCCACCCTAAACCACTCCCTCAGGGCGAGGGGCCCGGTCGGCGCTGTATGTACTCCCCCACTCTGAAACGACACCGCACGCGTTTTTTATcc                                               >  2:4287951/1‑150 (MQ=255)
                                                                                                     tGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTccc                                              <  1:3113722/150‑1 (MQ=255)
                                                                                                      cccTGCGGCACGGGCGCCCCGGGCGCACGGGACTGCGCCATCGGCGTTTGC‑GCCCCCGCCCCCGACCCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCCCTCTGAAACGACAGCGCACTCTTTTTTTCTcccc                                             <  2:5213397/149‑2 (MQ=255)
                                                                                                       ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTc                                            <  1:100272/150‑1 (MQ=255)
                                                                                                       ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCGCACCCGAACCCTCTACCACAGGGTGATGGGGCCGGTCGGCGCTGTATGTTCTCCCGCACGCCGAAACGCCACCCCACTCTTTTTTTCTACCTc                                            >  2:4782610/1‑150 (MQ=255)
                                                                                                       ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGC‑GCCCCCGCACCCCAACCCACAACCACCCGGAGAGGGGGCCCGTCGGCGCGGGGTGTACTCCCGCCCACAGAAAAGACCCCCCCCACgttttttttaccgc                                            >  1:4237815/1‑143 (MQ=255)
                                                                                                         tGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGc                                          <  2:4015845/150‑1 (MQ=255)
                                                                                                          gCGGCACGGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCCGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGcc                                         <  2:4528510/150‑1 (MQ=255)
                                                                                                            ggCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGcccc                                       <  1:5516497/150‑1 (MQ=255)
                                                                                                            ggCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTCCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGcccc                                       <  1:2148978/150‑1 (MQ=255)
                                                                                                                 tGGCTCCCATGGCGAAAGGGAATGCGGCATCGTCGTGTGC‑GCCCCCTCACCCTAACCCTCGCCCTCCGGGAGAGGGGACCGGTCGGCGCTGTATGTACTCCCTCCCTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgg                                  <  2:1923136/150‑1 (MQ=255)
                                                                                                                 tGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgg                                  <  1:4052951/150‑1 (MQ=255)
                                                                                                                 tGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgg                                  <  2:3368354/150‑1 (MQ=255)
                                                                                                                 tGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCGAACCCAAACCCACTCCCGCCGGGTGGGGGGCCCCGTCGGCGCAGTATGTACTCCCTCACACTGAAACGACACCGCCCTCTTTTTTTCTACCAAGccccgaccg                                  >  2:3884408/1‑145 (MQ=255)
                                                                                                                 tGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGC‑GCCCCCCCACCCCAACCCACGCCCCCACGGGGGCGGGACCCGTAGGCGCTGTCTTGACTGCCTCACACTGAAACGACACCGCACACTTTTTTTTTCCCTAGCCCCGCCgg                                  >  2:4264298/1‑150 (MQ=255)
                                                                                                                  ggCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg                                 <  2:5420105/150‑1 (MQ=255)
                                                                                                                    cgcccATGGCGACAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCCCCCTAACCCGCTCCCACAGGGAGAGGGGACCGTTCGGCGCTGTATGTACGCCCGCCCTCGGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga                               <  2:1116125/148‑1 (MQ=255)
                                                                                                                      cccATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgaga                             <  2:4854139/150‑1 (MQ=255)
                                                                                                                        cATGGCGAAAGGGAATGCGGCATAGTCGTGTGC‑GCCCCCCCCCCCCAACCCACCCCCCCCCGGGGGGGGGCCCCGTCGGCGCGGGGTGTACTCCCCCACTCTGAAAACCCACCCCACACATTTTTTCTCCCACGCCCCACCCCGGCGAgg                           >  2:1797033/1‑150 (MQ=255)
                                                                                                                         aTGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCTCCCCGGAGGGGCGAggg                          <  1:3707608/150‑1 (MQ=255)
                                                                                                                             cGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTACCCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgg                      <  1:5143754/150‑1 (MQ=255)
                                                                                                                                  gggAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCCTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                 <  1:264664/150‑1 (MQ=255)
                                                                                                                                  gggAATGCGGCATAGTCGTTTGC‑GCCCCCACACCCAAAACCACACCCACAGGGAGCGGGGACCCGTCGGCGCGGTATTTACTCCCACACTCGGAAACCACACCCCACTCTTTTTTTCTCCCGCGCCCCACCCGGGCGCGGGCCGGGGGGa                 >  2:3731098/1‑150 (MQ=255)
                                                                                                                                  gggAATGCGGCATAGTCGTTAGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                 <  2:3470667/150‑1 (MQ=255)
                                                                                                                                   ggAATGCGGCATAGTCGTTGGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAg                <  1:1360920/150‑1 (MQ=255)
                                                                                                                                    gattGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTTTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg               <  1:3293533/147‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGC‑GCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCA  >  NC_000913/3780873‑3781167

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