Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,781,038 T→A 17.4% intergenic (+21/‑177) 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,0380TA17.4% 990.0 / 49.6 410intergenic (+21/‑177)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (45/250);  new base A (57/2);  total (156/254)
Fisher's exact test for biased strand distribution p-value = 5.97e-35
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 414 total aligned reads displayed.

ATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGC  >  NC_000913/3780889‑3781186
                                                                                                                                                     |                                                                                                                                                    
  gAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCGCACCCaaa                                                                                                                                                    >  1:940994/1‑150 (MQ=255)
  gAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCCCACCCGaa                                                                                                                                                    >  1:8442462/1‑150 (MQ=255)
   aaaGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAAc                                                                                                                                                   <  2:3741979/150‑1 (MQ=255)
      gTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAAccca                                                                                                                                                >  1:2119732/1‑149 (MQ=255)
        ggCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCctct                                                                                                                                              <  2:3790081/150‑1 (MQ=255)
         gCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCtctc                                                                                                                                             <  2:4568108/150‑1 (MQ=255)
         gCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGGCCCCGCCCCCGAACCCtctc                                                                                                                                             >  1:3274569/1‑150 (MQ=255)
         gCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAAGGACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCtctc                                                                                                                                             <  2:5610479/150‑1 (MQ=255)
          cGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTcc                                                                                                                                            <  2:8011529/150‑1 (MQ=255)
           gATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCGCCCCCGAAACCGCTacc                                                                                                                                           >  2:1593017/1‑150 (MQ=255)
             tGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTc                                                                                                                                         <  1:7464687/150‑1 (MQ=255)
              gACGCTGACTGGCGCGAAATCTATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCa                                                                                                                                        <  2:1036718/150‑1 (MQ=255)
              gACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCa                                                                                                                                        <  1:2959878/150‑1 (MQ=255)
              gACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCGAACCCACACCCTCa                                                                                                                                        >  1:4907831/1‑150 (MQ=255)
                cGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAgg                                                                                                                                      <  1:4650533/150‑1 (MQ=255)
                 gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCCCCCGAACCCACTCCCTCAggg                                                                                                                                     >  1:3419747/1‑150 (MQ=255)
                 gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCCCACCCCAACCCCCAcccccccgg                                                                                                                                     >  2:6686768/1‑144 (MQ=255)
                   tGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCACACCCCAAACCACACCCACACGgag                                                                                                                                   >  1:7838746/1‑150 (MQ=255)
                      cTGGCGCGAAATCGATCAGCGACAGTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCAGGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAGGGAGAgg                                                                                                                                <  2:7093333/150‑1 (MQ=255)
                       tGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAggg                                                                                                                               <  1:5600495/150‑1 (MQ=255)
                       tGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAggg                                                                                                                               <  1:2457294/150‑1 (MQ=255)
                       tGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCAAACCCACAACCGCAGGGGGAggg                                                                                                                               >  1:6632151/1‑150 (MQ=255)
                         gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGa                                                                                                                             <  1:6756773/150‑1 (MQ=255)
                         gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGa                                                                                                                             <  2:5759352/150‑1 (MQ=255)
                         gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCACACCCGAACCCACTACCACAGGGCGAggggc                                                                                                                             >  1:7705132/1‑149 (MQ=255)
                          cgcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAGGGAGAgggtta                                                                                                                            <  2:1115741/150‑4 (MQ=255)
                          cgcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAgggtta                                                                                                                            <  2:5277881/150‑4 (MQ=255)
                          cgcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGAc                                                                                                                            <  1:7247670/150‑1 (MQ=255)
                            cgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAGGGAGAgggttagg                                                                                                                          <  2:3132981/150‑6 (MQ=255)
                            cgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAACCCGCACCCACAGGGAGAGGGGGccc                                                                                                                          >  1:688505/1‑149 (MQ=255)
                              aaaTCGATCAGCGAAATAACGCAAGATTCGCTGGTTCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAATGGAATTCGGCATAGTCGTTTTCGACCCCACACACTCACCCACTCCCGCCGGGTGAGGGGAccctg                                                                                                                        >  2:3015301/1‑147 (MQ=255)
                                 tCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCgg                                                                                                                     <  2:6631835/150‑1 (MQ=255)
                                 tCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCgg                                                                                                                     <  2:3928897/150‑1 (MQ=255)
                                  cGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGTGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgc                                                                                                                    <  1:6449404/150‑1 (MQ=255)
                                  cGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgc                                                                                                                    <  1:2219957/150‑1 (MQ=255)
                                   gATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgcg                                                                                                                   <  1:5017830/150‑1 (MQ=255)
                                     tCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCt                                                                                                                 <  2:2496821/150‑1 (MQ=255)
                                       aGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCACACCCGAACCCACACCCACAGGGAGAGGGGACCGTTCGGCGCtgt                                                                                                               >  1:8045119/1‑150 (MQ=255)
                                        gCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgta                                                                                                              <  2:3532413/150‑1 (MQ=255)
                                        gCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCCCCCCAAACCACAACCACCCGGGGCGGGGACCCTTCGgcgcggga                                                                                                              >  1:4982323/1‑146 (MQ=255)
                                         cGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtat                                                                                                             <  1:8278519/150‑1 (MQ=255)
                                         cGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCCACCCCAAACCCGCACCCGCCGGGAGGGGGGGCCGGTAGgcgcgggtt                                                                                                             >  1:368843/1‑145 (MQ=255)
                                             aTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTAc                                                                                                         <  1:3488511/150‑1 (MQ=255)
                                             aTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTAc                                                                                                         <  1:5775960/150‑1 (MQ=255)
                                              ttACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACt                                                                                                        <  2:4369157/150‑1 (MQ=255)
                                               tACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTc                                                                                                       <  1:596765/150‑1 (MQ=255)
                                                  gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTTCGCCCCCTCACCCTAACCCACTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCACt                                                                                                    >  2:4822239/1‑150 (MQ=255)
                                                  gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCt                                                                                                    <  2:448434/150‑1 (MQ=255)
                                                    aaGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGACTGCGGCATAGTCGTGTGCGCCCCCACACCCAAAACCCAGACCGCAGGGCGATGGGACCCGTAGGCGCGGGATGTGCTCCCTcc                                                                                                  >  2:524376/1‑149 (MQ=255)
                                                      gATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCAct                                                                                                <  2:4727255/150‑1 (MQ=255)
                                                        ttCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCAAACCCACTCCCACAGGGAGAGGGGACCGTTAGGCGCGGTATGTACTCCCGCAcacg                                                                                              >  2:5537565/1‑147 (MQ=255)
                                                            cTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaaa                                                                                          <  2:8259162/150‑1 (MQ=255)
                                                            cTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaaa                                                                                          <  2:6294815/150‑1 (MQ=255)
                                                               gTGCTGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGa                                                                                       <  1:5325478/150‑1 (MQ=255)
                                                               gTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCACCCAAACCCACACCCACAGGGGGCGGGGGCCGTTCGGCGCGGTATGTTCTCCCTCACACCGaaaaca                                                                                       >  2:6153980/1‑147 (MQ=255)
                                                                 gCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGTGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGaca                                                                                     <  2:6531708/150‑1 (MQ=255)
                                                                 gCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGaca                                                                                     <  2:2821454/150‑1 (MQ=255)
                                                                 gCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCCCCCCAACCCACACCCACCCGGCCCGGGGCCCCGTCCGCGCCGGATGGACACCCGCACACAGAAACGacc                                                                                     >  2:3398181/1‑149 (MQ=255)
                                                                           ggTAAAGAGGTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCtt                                                                           <  2:1532419/150‑1 (MQ=255)
                                                                             tAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTGGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCtttt                                                                         <  1:1023524/150‑1 (MQ=255)
                                                                             tAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCCCCCAAACCCCCACCCTCAGGGGGTGGGGGCCCTTCGGCGCAGTATGTACTCCCACAAACTGAAACGACACCGCAccgagtt                                                                         >  2:7724049/1‑144 (MQ=255)
                                                                             tAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCGAACCCACTCCCACAGGGAGAGGGGGCCGGTCGGCGCGGTATGTACTCCCGCACTCAGAAACGACACCCCACACattt                                                                         >  2:6366260/1‑150 (MQ=255)
                                                                              aaaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttttt                                                                        <  1:5909587/150‑1 (MQ=255)
                                                                                agagTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAAACCACAACCGCAGGGGGAGGGGGCCCGTGGGCGCGGGATTTAAAACCACACTCTGTAAAGACACCCCCCTCatttttt                                                                      >  1:5457707/1‑150 (MQ=255)
                                                                                       ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCAAACCCGCTCCCACAGGGCGAGGGGCCCCGTCGGCGCTGTTTGTTATACCGCACTCTGAAAAGACACCGCACGCTTTTTTTCTcccgc                                                               >  2:1285935/1‑148 (MQ=255)
                                                                                        cTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCg                                                              <  2:4095132/150‑1 (MQ=255)
                                                                                        cTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCg                                                              <  1:1433160/150‑1 (MQ=255)
                                                                                         tGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCACAACCAATCACACTCCCACCGGGGGCGGGGCCCCGTAGGCGCGGTTTGTTCGCCCACACACTGGAAAGACACCCCACACTTTTTTTATACCgcgc                                                             >  2:2661099/1‑150 (MQ=255)
                                                                                              cACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGTGCCCCCTCACCCTAACCCTCTACCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACACTGCAAAGACACCGCCCTCTTTTTTTCTCCCTCGccccac                                                        >  2:7687820/1‑148 (MQ=255)
                                                                                                  ggCTCCCCTGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg                                                    <  2:5000103/150‑1 (MQ=255)
                                                                                                  ggCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg                                                    <  2:2589085/150‑1 (MQ=255)
                                                                                                  ggCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg                                                    <  2:8355534/150‑1 (MQ=255)
                                                                                                   gcgcccAGGGCGAAAGGGAATGCGGCATAGGCGTTTGCGCCCCCTCACCCTACCCCTCTCCCTCCGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgggg                                                   <  2:5744737/147‑1 (MQ=255)
                                                                                                   gCTCCCATGGCGAAAGGGAATGCGACATAGTCGTGTGCGCCCCCACCCCCAAAACCACTCCCGCAGGGAGAGGGGACCCGTCCGCGCGGGATGTACGACCGCACTCTGAAACGACCCCGCACACTTTTTTTCTCCCACGCCCCAccccgg                                                   >  2:7005718/1‑146 (MQ=255)
                                                                                                    cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga                                                  <  2:14059/150‑1 (MQ=255)
                                                                                                    cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga                                                  <  1:3907334/150‑1 (MQ=255)
                                                                                                      cccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgaga                                                <  1:270109/150‑1 (MQ=255)
                                                                                                      cccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgaga                                                <  1:8333480/150‑1 (MQ=255)
                                                                                                       ccATGGCGAATGGGAATGGGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGATTGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTTAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGagag                                               <  1:1825566/150‑1 (MQ=255)
                                                                                                       ccATGGCGAAAGGGAATGCGGCATAGTCGTTTGGGCCCCCACACCCGAACCCACACCCAAAGGGAGAGGGGACCGGTCGGCGCGGTATGTACTACCACACTCTGAAACGAAACCGCACACATTTTTTATACCGCGCCCCACCGGGGggag                                               >  2:5187648/1‑150 (MQ=255)
                                                                                                         aTGGCGAAAGGGAATGCGGCATAGTCGTTGGCGCCCCCTCACCCTACCCCTCTCCCTCCGGGAGAGGGGACCGTGCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAggg                                             <  1:6683846/150‑1 (MQ=255)
                                                                                                            gCGAAAGGGAATGCGGCATAGTCGTTTTCGCCCCCTCAACCGAAACCACTCCCACACGGAGCTGGGGCCCGTGGGCGCGGGATGTTCAACCCCACCCTGAAACGACCCCCCCCACGTTTTTTCTCCCCCGCCCCACCGCGGCGGGGGGCg                                          >  2:1864557/1‑150 (MQ=255)
                                                                                                             cGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgg                                         <  1:1258372/150‑1 (MQ=255)
                                                                                                              gAACGGGAATGCGGCATAGTCGTGTGCGCCCCCTCACCCTAACCCCCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGGACTCCCTCACTCTGAAACGACACCGCACTCTTGTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                        <  2:5359149/150‑1 (MQ=255)
                                                                                                                  gggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                                    <  1:6233396/150‑1 (MQ=255)
                                                                                                                    gAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                  <  2:5075972/150‑1 (MQ=255)
                                                                                                                    gAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCGCCCTCGCCCCTCCGGGGAGATGGCCGGGGTGAgg                                  <  1:4229927/150‑1 (MQ=255)
                                                                                                                     caTGCGGCCGAGTCGTGTGCGCCCCCGCCCCCGACCCCTCGCCCTCCGGGAGAGGGGACCGTGCGGCGCTGTATGTACTCCCGCCCTCTGAAACGACCCCGCACTCTTTTTTTCGCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAggg                                 <  2:6819514/149‑1 (MQ=255)
                                                                                                                          ggCAGAGTCGTTTGCGCCCCCTCCCCCTAACCCTCTCCCTCAGGGAGAGGTGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                            <  2:7815787/150‑1 (MQ=255)
                                                                                                                                 tCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGGATGTACGCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCa                     >  2:4268956/1‑150 (MQ=255)
                                                                                                                                 tCGTTTGCGCCCCCTCACCCAAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCGGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGGGAGGGGAAAAGGCCGcc                     >  2:4646553/1‑149 (MQ=255)
                                                                                                                                    tttGCGCCCCCTCACCCGAACCCTCACCCTCAGGGAGAGGGGACCCTTCGGCGCTGTATGTACTCCCACACACAGAAACCACACCGCACTCATTTTTTCTCCCTCGCCCCACCGGGGAGAGGGCCGGGGGGAGGGGAAAAGGCCGCACTg                  >  2:5947419/1‑150 (MQ=255)
                                                                                                                                         gcCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCAc             >  2:7475368/1‑150 (MQ=255)
                                                                                                                                               tCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAttt       >  2:7229314/1‑150 (MQ=255)
                                                                                                                                                   ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg   >  1:4527583/1‑150 (MQ=255)
                                                                                                                                                   ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg   >  1:5657313/1‑150 (MQ=255)
                                                                                                                                                   ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg   >  1:523807/1‑150 (MQ=255)
                                                                                                                                                    cTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGc  >  2:5421382/1‑150 (MQ=255)
                                                                                                                                                    cTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGc  >  2:2894383/1‑150 (MQ=255)
                                                                                                                                                    cTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGAGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGc  <  2:2076856/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGC  >  NC_000913/3780889‑3781186

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: