Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,439,588 G→C 5.9% intergenic (‑79/+28) yhdN ← / ← rplQ DUF1992 family protein/50S ribosomal subunit protein L17

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,439,5880GC 5.9% 1598.2 / 73.0 441intergenic (‑79/+28)yhdN/rplQDUF1992 family protein/50S ribosomal subunit protein L17
Reads supporting (aligned to +/- strand):  ref base G (210/205);  new base C (5/21);  total (215/226)
Fisher's exact test for biased strand distribution p-value = 2.01e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

Only 100 of 459 total aligned reads displayed.

AGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTC  >  NC_000913/3439439‑3439736
                                                                                                                                                     |                                                                                                                                                    
      tCAAACTCACCGTTCGCTTGCGCTTCTGCTATATGGCGCTCGGCCCACTGGTCAAGTAACCACATAACTACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCgg                                                                                                                                                <  2:6357947/150‑1 (MQ=255)
        aaaCTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGt                                                                                                                                              >  1:6780892/1‑150 (MQ=255)
         aaCTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTTGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGtt                                                                                                                                             >  1:502578/1‑150 (MQ=255)
          aCTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCGGCCCACTGGTCAAGTAACCACATAAATACGCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATACGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGttt                                                                                                                                            <  1:2294489/150‑1 (MQ=255)
          aCTCACCTGTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGCAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGttt                                                                                                                                            <  2:4746234/150‑1 (MQ=255)
           cTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGtttt                                                                                                                                           <  2:6783714/150‑1 (MQ=255)
            tCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAACTACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGttttt                                                                                                                                          <  2:1940626/150‑1 (MQ=255)
            tCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTATGCCCAAAGGTAAAGTAACCACAAAAAAAAAACAGGGGATGAACAAAAAGAGTACAGCTATTTTAGATAAGTGGGGATACTACGGGTATAAAAAAAAACGCCGGGGCGGGttttt                                                                                                                                          >  1:2801973/1‑150 (MQ=255)
            tCACCTTTCGCTTGCGCTTCTGCTATAAGGCGCTCTGCCCACTGGTCAAGTAACCACATAAAAACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACAACGGGTATAAAAAAAACCGCCGGGGCGGGttttt                                                                                                                                          >  1:1420499/1‑150 (MQ=255)
             cACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGtttttt                                                                                                                                         <  2:4246569/150‑1 (MQ=255)
             cACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGtttttt                                                                                                                                         <  2:2908434/150‑1 (MQ=255)
             cACCTTTCGCTTGCGCTCCTGCTATATGCCGCTCTGCCCACTGGACAAGTAACCACATAACTACTCCAGGGTATTAACAAAAATAGTACAGCTTTTGTAGATACGTGGGTATACTACGGTTATAAAAAAACCCGCCGGGGCGGGtttttt                                                                                                                                         <  1:382534/150‑1 (MQ=255)
             cACCTGTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCCGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGtttttt                                                                                                                                         <  2:5072060/150‑1 (MQ=255)
              aCCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGttttttt                                                                                                                                        <  1:764107/150‑1 (MQ=255)
                ctgtCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTAc                                                                                                                                      >  1:3325926/4‑150 (MQ=255)
                cTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTAc                                                                                                                                      >  2:3232124/1‑150 (MQ=255)
                cTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACAGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCCCCGGGGCGGGTTTTTTTAc                                                                                                                                      >  1:6053212/1‑150 (MQ=255)
                  ttCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGt                                                                                                                                    >  2:3971165/1‑150 (MQ=255)
                       ttGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCtt                                                                                                                               >  2:3611024/1‑150 (MQ=255)
                        tGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTc                                                                                                                              >  2:1367789/1‑150 (MQ=255)
                         gcgcTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCa                                                                                                                             >  2:5318518/1‑150 (MQ=255)
                           gcTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGa                                                                                                                           >  1:6891174/1‑150 (MQ=255)
                           gcTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGa                                                                                                                           >  1:2468010/1‑150 (MQ=255)
                             ttCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGAta                                                                                                                         >  1:4238873/1‑149 (MQ=255)
                              tCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTa                                                                                                                        <  2:1040764/150‑1 (MQ=255)
                              tCTGCTATATGGCGCACTGCCCACTGGTCAAGTAACCACCTACATACTCCAGGGGATGAACAAAAAGAGCACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTa                                                                                                                        <  2:1017110/150‑1 (MQ=255)
                                tGCGATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTAct                                                                                                                      <  2:3678829/150‑1 (MQ=255)
                                  cTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTActct                                                                                                                    >  1:2940366/1‑150 (MQ=255)
                                    atatGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCCGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgc                                                                                                                  <  2:5946560/150‑1 (MQ=255)
                                       tGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTCCGGGTATAAAAAAAACCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTCCTCTgcagc                                                                                                               >  1:2113096/1‑150 (MQ=255)
                                          cgcTCGGCCCACTGGTCAAGTAACCACATCAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAACAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCagcagc                                                                                                            <  2:3354651/150‑1 (MQ=255)
                                               tGCCCACTGGTCAAGTAACCACATAAACACTCCAGGGGATGAACAAAAAGAGTACAGCTAGTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCGGCTTCTg                                                                                                       <  2:5464460/150‑1 (MQ=255)
                                                gCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGc                                                                                                      <  1:6491716/150‑1 (MQ=255)
                                                 cccACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCt                                                                                                     >  1:4541599/1‑150 (MQ=255)
                                                   cACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCttt                                                                                                   <  2:2265936/150‑1 (MQ=255)
                                                    aCTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTtc                                                                                                  >  1:5391935/1‑150 (MQ=255)
                                                    aCTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTtc                                                                                                  >  2:7150499/1‑150 (MQ=255)
                                                       ggTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGAGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTctct                                                                                               <  2:7113437/150‑1 (MQ=255)
                                                        gTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTCCAGCTCTTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTg                                                                                              >  2:3840631/1‑150 (MQ=255)
                                                         tCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGa                                                                                             >  1:4794172/1‑150 (MQ=255)
                                                          cAAGTAACCACATAAATACTCCAGGGGCTGAACCAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTGTTTACGTTGCTTCAGAGTACTCTGCCGCAGCTTCTGCTTTCTCTGac                                                                                            >  1:6156293/1‑149 (MQ=255)
                                                          cAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGaa                                                                                            <  2:4183368/150‑1 (MQ=255)
                                                           aaGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAAc                                                                                           >  2:1698021/1‑150 (MQ=255)
                                                           aaGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAAc                                                                                           <  1:4207735/150‑1 (MQ=255)
                                                           aaGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAAc                                                                                           >  1:6558792/1‑150 (MQ=255)
                                                           aaGTAACCACATAAAAACTCCAGGGGATGAACAAAAAGAGTACACCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAAc                                                                                           <  1:6606146/150‑1 (MQ=255)
                                                             gTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGAAACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGa                                                                                         <  2:6703822/150‑1 (MQ=255)
                                                                aCCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCa                                                                                      <  1:1008176/150‑1 (MQ=255)
                                                                 ccACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCaa                                                                                     >  1:6412205/1‑150 (MQ=255)
                                                                  cacaTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAc                                                                                    <  2:894009/150‑1 (MQ=255)
                                                                   acaTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAcc                                                                                   >  2:3158053/1‑150 (MQ=255)
                                                                   acaTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAcc                                                                                   <  1:3536144/150‑1 (MQ=255)
                                                                   acaTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAcc                                                                                   >  1:1999354/1‑150 (MQ=255)
                                                                     aTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAg                                                                                 >  1:1041660/1‑150 (MQ=255)
                                                                      caaaTACTCCAGGGGATGAACACAATGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGAATACTCCGCAGCAGCTTATGCTTTCTCTGAACGATCACCCAGc                                                                                <  2:2671279/149‑1 (MQ=255)
                                                                       aaaTACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCt                                                                               <  1:4949375/150‑1 (MQ=255)
                                                                        aagaCTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTc                                                                              <  1:3929641/147‑1 (MQ=255)
                                                                               cAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTaa                                                                       >  2:6492320/1‑150 (MQ=255)
                                                                                aGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAg                                                                      >  2:2158112/1‑150 (MQ=255)
                                                                                aGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAg                                                                      >  2:4983704/1‑150 (MQ=255)
                                                                                 ggggATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGc                                                                     >  1:69631/1‑150 (MQ=255)
                                                                                          cacaaaGAGTACAGCGATTGTAGATAAGTGGGGATACTACGGGTAGAAACAAACCCGCCGGGGCGGGATTATTTACGGTGCTTCAGATTACTCGGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgcgc                                                            <  2:5148160/147‑1 (MQ=255)
                                                                                          cAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgcgc                                                            <  1:3399460/150‑1 (MQ=255)
                                                                                           aaaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGcgcg                                                           >  1:6721076/1‑150 (MQ=255)
                                                                                           aaaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGcgcg                                                           >  2:4452686/1‑150 (MQ=255)
                                                                                           aaaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAAAGATCAACCAGCTCGATGTAAGCCATCGGcgcg                                                           >  1:5232932/1‑150 (MQ=255)
                                                                                             aaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGtt                                                         >  2:2743219/1‑150 (MQ=255)
                                                                                             aaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGtt                                                         >  2:4857893/1‑150 (MQ=255)
                                                                                                gagTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTc                                                      <  2:2380906/150‑1 (MQ=255)
                                                                                                gagTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTc                                                      <  2:2952107/150‑1 (MQ=255)
                                                                                                  gTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTATTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGc                                                    <  1:4614248/150‑1 (MQ=255)
                                                                                                   tACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGcc                                                   <  1:5539297/150‑1 (MQ=255)
                                                                                                    aCAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCt                                                  <  1:1534492/150‑1 (MQ=255)
                                                                                                      aGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGc                                                >  2:5705791/1‑150 (MQ=255)
                                                                                                       gCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCa                                               >  1:1215923/1‑150 (MQ=255)
                                                                                                         tATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACg                                             <  1:835930/150‑1 (MQ=255)
                                                                                                         tATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACg                                             >  2:1485194/1‑150 (MQ=255)
                                                                                                          aTTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACgg                                            <  1:4869465/150‑1 (MQ=255)
                                                                                                            tGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGaa                                          <  1:5133392/150‑1 (MQ=255)
                                                                                                             gTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAg                                         >  1:5278716/1‑150 (MQ=255)
                                                                                                              tAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGc                                        <  1:5623161/150‑1 (MQ=255)
                                                                                                                 aTAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTAAGATTACTCTGCAGCAGCTTCTGCTTTCTATGAACGATCAACCAGATCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCcac                                     >  2:3999402/1‑150 (MQ=255)
                                                                                                                     gTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACtt                                 <  1:359615/150‑1 (MQ=255)
                                                                                                                      tGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTc                                >  1:5812881/1‑150 (MQ=255)
                                                                                                                               tacGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGa                       <  2:4612716/150‑1 (MQ=255)
                                                                                                                                  gggTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCCGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAgtg                    >  1:1307333/1‑150 (MQ=255)
                                                                                                                                    gTATAAAAAATCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTa                  <  1:4683391/150‑1 (MQ=255)
                                                                                                                                    gTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGGTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTa                  >  1:5933088/1‑150 (MQ=255)
                                                                                                                                      ataAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAac                <  1:7111664/150‑1 (MQ=255)
                                                                                                                                      ataAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAac                <  1:1597476/150‑1 (MQ=255)
                                                                                                                                         aaaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGGTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAaccac             >  2:1121927/1‑150 (MQ=255)
                                                                                                                                         aaaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCAACGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAaccac             >  1:5776611/1‑150 (MQ=255)
                                                                                                                                          aaaaaCCCGCCGGGGCGGGTTTGTTTACGTTGCTTCAGATTACTCTGCAGCAGCGTCTGCTTTCTCTGAACGATCAACCAGCTCGCTGTAAGCCATCGGCGCGTGGTCGCCTGCACGGAAGCCACACGTCAGAATACGAGTGTAaccacc            >  2:7135248/1‑150 (MQ=255)
                                                                                                                                          aaaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAaccacc            >  2:6898916/1‑150 (MQ=255)
                                                                                                                                           aaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCg           >  2:6936804/1‑150 (MQ=255)
                                                                                                                                               cccGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCAc       <  2:282081/150‑1 (MQ=255)
                                                                                                                                                 cgccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACgg     <  1:1747198/150‑1 (MQ=255)
                                                                                                                                                 cgccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACgg     <  1:3201410/150‑1 (MQ=255)
                                                                                                                                                   ccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCt   <  2:5573971/150‑1 (MQ=255)
                                                                                                                                                    cgGGGCGGGTTTTTTTACGGTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCATTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTc  <  2:6198920/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGGTTATCAAACTCACCTTTCGCTTGCGCTTCTGCTATATGGCGCTCTGCCCACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCGCCTGCACGGAAGCCACACTTCAGAATACGAGTGTAACCACCGGCACGGCTC  >  NC_000913/3439439‑3439736

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: