Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,781,022 T→G 10.4% 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.4% 1209.4 / 51.2 451intergenic (+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 (112/283);  new base G (52/4);  total (164/287)
Fisher's exact test for biased strand distribution p-value = 3.43e-21
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 463 total aligned reads displayed.

TCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTG  >  NC_000913/3780873‑3781170
                                                                                                                                                     |                                                                                                                                                     
tCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGtt                                                                                                                                                       >  1:7099243/1‑150 (MQ=255)
 cTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGttt                                                                                                                                                      <  1:8205887/150‑1 (MQ=255)
  tGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTg                                                                                                                                                     >  2:4595129/1‑150 (MQ=255)
  tGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTg                                                                                                                                                     >  2:5040332/1‑150 (MQ=255)
   gATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTgc                                                                                                                                                    <  1:6143454/150‑1 (MQ=255)
       gAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTCAAGAGTTGCCTGCGGCACTGGCTCCCAGGGCGAAAGGGAATGCGGCATAGTCGTTTGCGccc                                                                                                                                                >  1:1698185/1‑150 (MQ=255)
       gAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGccc                                                                                                                                                >  1:7166471/1‑150 (MQ=255)
        aaaaaGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGcccc                                                                                                                                               <  2:4245209/150‑1 (MQ=255)
        aaaaaGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGcccc                                                                                                                                               <  2:4797690/150‑1 (MQ=255)
         aaaaGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGccccc                                                                                                                                              <  2:7195568/150‑1 (MQ=255)
           aaGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGcccccac                                                                                                                                            >  2:2125226/1‑148 (MQ=255)
            agagATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCa                                                                                                                                           <  2:292757/150‑1 (MQ=255)
             gagaTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTctc                                                                                                                                          <  1:3160073/150‑3 (MQ=255)
             gagaTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCAc                                                                                                                                          <  1:1976197/150‑1 (MQ=255)
             gagaTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCAc                                                                                                                                          <  1:2777433/150‑1 (MQ=255)
             gagaTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCAc                                                                                                                                          <  1:5139582/150‑1 (MQ=255)
              agaTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCAcc                                                                                                                                         <  2:6982506/150‑1 (MQ=255)
              agaTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCAcc                                                                                                                                         <  2:4001897/150‑1 (MQ=255)
              agaTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCAcc                                                                                                                                         <  2:1313656/150‑1 (MQ=255)
               gaTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTAGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGcccccccccac                                                                                                                                        >  2:733771/1‑144 (MQ=255)
                   aaaGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCGAAc                                                                                                                                    >  1:258838/1‑150 (MQ=255)
                      gTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAAccca                                                                                                                                 >  1:2119732/1‑149 (MQ=255)
                       tGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCACACCCTAAaccac                                                                                                                                >  2:4498177/1‑148 (MQ=255)
                        ggCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCctct                                                                                                                               <  1:4488573/150‑1 (MQ=255)
                        ggCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCctct                                                                                                                               <  1:3687188/150‑1 (MQ=255)
                         gCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCtctc                                                                                                                              <  2:5610438/150‑1 (MQ=255)
                         gCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCtctc                                                                                                                              <  2:1014414/150‑1 (MQ=255)
                          cGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTcc                                                                                                                             <  2:6797737/150‑1 (MQ=255)
                          cGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTcc                                                                                                                             <  2:8011529/150‑1 (MQ=255)
                           gATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAACCCACAccc                                                                                                                            >  1:4976102/1‑150 (MQ=255)
                           gATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCGCCCCCGAAACCGCTacc                                                                                                                            >  2:1593017/1‑150 (MQ=255)
                           gATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCCCCCCCCAAaccccgacc                                                                                                                            >  2:3850993/1‑144 (MQ=255)
                            aTGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCt                                                                                                                           <  2:4682869/150‑1 (MQ=255)
                              gACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCa                                                                                                                         <  2:124793/150‑1 (MQ=255)
                              gACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCGAACCCACACCCTCa                                                                                                                         >  1:4907831/1‑150 (MQ=255)
                                cGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAgg                                                                                                                       <  1:4650533/150‑1 (MQ=255)
                                   tGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCACACCCCAAACCACACCCACACGgag                                                                                                                    >  1:7838746/1‑150 (MQ=255)
                                    gACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGgaga                                                                                                                   <  1:7586516/150‑1 (MQ=255)
                                    gACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGgaga                                                                                                                   <  2:3030255/150‑1 (MQ=255)
                                    gACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGgaga                                                                                                                   <  2:7045332/150‑1 (MQ=255)
                                      cTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAGGGAGAgg                                                                                                                 <  2:3211917/150‑1 (MQ=255)
                                       tGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAGGGAGAggg                                                                                                                <  2:76427/150‑1 (MQ=255)
                                       tGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAggg                                                                                                                <  2:4940847/150‑1 (MQ=255)
                                        ggCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAGGGAGAgggt                                                                                                               <  1:7308595/150‑2 (MQ=255)
                                         gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTTCGCCCCCACACCCAAACCCACACCCACACGGAGGGGGGa                                                                                                              >  1:5585246/1‑150 (MQ=255)
                                         gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGa                                                                                                              <  2:5759352/150‑1 (MQ=255)
                                           gcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAgggttag                                                                                                            <  2:3744845/150‑5 (MQ=255)
                                            cgAACTCGAGCAGCGAAAGTCCGCACGCTTCGCTGGTGCCGGGGCTGGGTAACGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAGGGAGAgggttagg                                                                                                           <  2:845111/150‑6 (MQ=255)
                                            cgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTTCGCCCCCTCACCCAAACCCACTCCCACAGGGAGAGGGGAccc                                                                                                           >  1:3442875/1‑149 (MQ=255)
                                            cgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAGGGAGAgggttagg                                                                                                           <  2:5711479/150‑6 (MQ=255)
                                                  cGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGTGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgc                                                                                                     <  1:6449404/150‑1 (MQ=255)
                                                  cGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgc                                                                                                     <  1:2219957/150‑1 (MQ=255)
                                                   gATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcg                                                                                                    <  1:5017830/150‑1 (MQ=255)
                                                     tCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCTAACCCACACCCCCAGGGAGAGGGGCCCGTTGGgcgcg                                                                                                  >  1:5751358/1‑149 (MQ=255)
                                                       aGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCCCCCTCAGGGAGAGGGGACCGTTCGGCGCtgt                                                                                                <  2:6360724/150‑1 (MQ=255)
                                                        gCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgta                                                                                               <  2:3532413/150‑1 (MQ=255)
                                                         cGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGGTCGGCGCtgtat                                                                                              <  2:5369491/150‑1 (MQ=255)
                                                          gAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCACACCCAAAAACACACCCAAAGGGATAGGGGGCCGATCGGCGCggtttt                                                                                             >  1:6651387/1‑147 (MQ=255)
                                                            aaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtatgta                                                                                           <  1:3280338/150‑1 (MQ=255)
                                                            aaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAAACCACTCCCACACGGAGAGGGGACCCGTAGGCGCggtatgta                                                                                           >  1:8413134/1‑150 (MQ=255)
                                                              ttACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACt                                                                                         <  2:4369157/150‑1 (MQ=255)
                                                                        ttCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCAAACCCACTCCCACAGGGAGAGGGGACCGTTAGGCGCGGTATGTACTCCCGCAcacg                                                                               >  2:5537565/1‑147 (MQ=255)
                                                                         tCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTg                                                                              <  1:6238733/150‑1 (MQ=255)
                                                                         tCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTg                                                                              <  1:5416908/150‑1 (MQ=255)
                                                                            cTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaaa                                                                          <  2:2531169/150‑1 (MQ=255)
                                                                            cTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaaa                                                                          <  1:2346293/150‑1 (MQ=255)
                                                                            cTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaaa                                                                          <  2:6294290/150‑1 (MQ=255)
                                                                             tGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAAc                                                                         <  1:5823291/150‑1 (MQ=255)
                                                                               gTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGa                                                                       <  1:5325376/150‑1 (MQ=255)
                                                                                           ggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCTAAACCACTCCCTCAGGGTGTGGGGACCGGTAGGGGGGGTTTGTACTGCCTCACTCTGAAACCACACCGcccaatt                                                           >  2:5087825/1‑144 (MQ=255)
                                                                                             tAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTGGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCtttt                                                         <  1:1023524/150‑1 (MQ=255)
                                                                                             tAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCCCCCAAACCCCCACCCTCAGGGGGTGGGGGCCCTTCGGCGCAGTATGTACTCCCACAAACTGAAACGACACCGCAccgagtt                                                         >  2:7724049/1‑144 (MQ=255)
                                                                                                  agTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtct                                                    <  1:561856/150‑1 (MQ=255)
                                                                                                  agTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtct                                                    <  2:6576874/150‑1 (MQ=255)
                                                                                                  agTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCCCCCCAAACCACAACCTCAGGGGGGGGGGTCCGGTCGGCGCTGTATTTACGCCCCCACAAAGGAAACACACCCCACGCtttttttat                                                    >  2:8327812/1‑148 (MQ=255)
                                                                                                   gTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtctc                                                   <  1:7995054/150‑1 (MQ=255)
                                                                                                    ttGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTcc                                                  <  2:5332119/150‑1 (MQ=255)
                                                                                                        cTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCg                                              <  1:1433160/150‑1 (MQ=255)
                                                                                                         tGCGGCACTGGCTCCCATGGCGAATGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACGCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGc                                             <  2:6656758/150‑1 (MQ=255)
                                                                                                         tGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTTCGCCCCCCCAACCAAAACCACAACCGCAGGGGGGGGGGGCCCGTCGGGGCGGGTTTTTATCCCACCCGCTGTAAACGCACCGCCCGCTTTTTTTTTccccccc                                             >  2:259623/1‑146 (MQ=255)
                                                                                                                 tGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgg                                     <  2:1453454/150‑1 (MQ=255)
                                                                                                                 tGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgg                                     <  1:931563/150‑1 (MQ=255)
                                                                                                                   gCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTTCCCCCCCACACCCAAACCCAAAACCACAGGGAGGGGGGGCCCTTCCGCGCTGTATTTACTCCCACACACTGAAAACACAACCCACACTTTTTTTAACCCACGCCCCACCcggg                                   >  2:2794773/1‑150 (MQ=255)
                                                                                                                    cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga                                  <  2:6442072/150‑1 (MQ=255)
                                                                                                                    cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga                                  <  1:830057/150‑1 (MQ=255)
                                                                                                                       ccATGGCGAAAGGGAATGCGGCATAGTCGTTTTCGCCCCCACACCCAAAACCGCACCCACAGGGGGTGGGGGCCCTTCGGCGCTGTATTTACTCCCACACACAGAAACGACACCGCACACTTTTTTTATCCCGCGCCCCACCCGGGcgag                               >  2:1962885/1‑150 (MQ=255)
                                                                                                                       ccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGagag                               <  1:6774600/150‑1 (MQ=255)
                                                                                                                         ctGGCGAAAGGGAATGCGGCATAGTCGTTGGCGCCCCCTCACCCTCCCCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAggg                             <  1:6684003/149‑1 (MQ=255)
                                                                                                                         aTGGCGAAAGGGAATGCGGCATAGTCGTTTTCGCCCCCACACCCGAACCCAAAACCCAACGGTGTGGGGAACCGTCGGCGCAGGTTGTAACACCGCACGCAGAAACCACACCGCCAACATTTTTTATCCCGCCCCCCGCCCGGGCGgggg                             >  2:1083766/1‑150 (MQ=255)
                                                                                                                             cGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgg                         <  1:3889986/150‑1 (MQ=255)
                                                                                                                             cGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgg                         <  1:1258372/150‑1 (MQ=255)
                                                                                                                              gAAAGGGAATGCGGCATAGTCGTTTGTGCCCCCGCACCCGAAACCACTACCTCAGGGAGAGGGGACCGTTCGGCGCAGTATTTACTCCCACAAACAGAAACGACAACGCACACATTTTTTATACCTCGCCCCTCCGGGGAGAGGGCCggg                        >  2:8255165/1‑150 (MQ=255)
                                                                                                                                aaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                      <  1:5959875/150‑1 (MQ=255)
                                                                                                                                    gAATGCGGCATAGTCGTTTTCGCCCCCCCACCCAAACCCTATCCCTCCGGGAGAGGGGACCCTTCGGCGCTGTATGTAATACCTCACTATGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCAACGGGGGGAGGGGCGGGGGGAgg                  >  2:5255296/1‑150 (MQ=255)
                                                                                                                                    gAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCCGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                  <  2:3419747/150‑1 (MQ=255)
                                                                                                                                    gAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCGCCCTCGCCCCTCCGGGGAGATGGCCGGGGTGAgg                  <  1:4229927/150‑1 (MQ=255)
                                                                                                                                    gAATGCGGCATAGTCGTTTGCGCCCCCACCCCCCAACCCACACCCCCCCGGCGCGGGGACCCGGCCGCGCCGTAAGTACACCCACACACAGAAACGACACCCCACACATTTTTTCTCCCGCCCCCCACCCGGGAGAGGGCCGGGGGGAgg                  >  2:522659/1‑150 (MQ=255)
                                                                                                                                          ggCAGAGTCGTTTGCGCCCCCTCCCCCTAACCCTCTCCCTCAGGGAGAGGTGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa            <  2:7815787/150‑1 (MQ=255)
                                                                                                                                                 tCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGGATGTACGCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCa     >  2:4268956/1‑150 (MQ=255)
                                                                                                                                                    tttGCGCCCCCTCACCCGAACCCTCACCCTCAGGGAGAGGGGACCCTTCGGCGCTGTATGTACTCCCACACACAGAAACCACACCGCACTCATTTTTTCTCCCTCGCCCCACCGGGGAGAGGGCCGGGGGGAGGGGAAAAGGCCGCACTg  >  2:5947419/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTG  >  NC_000913/3780873‑3781170

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: