Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,282,832 C→G 5.6% intergenic (+31/+108) rplY → / ← yejK 50S ribosomal subunit protein L25/spermidine nucleoid‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,282,8320CG 5.6% 366.7 / 10.8 107intergenic (+31/+108)rplY/yejK50S ribosomal subunit protein L25/spermidine nucleoid‑associated protein
Reads supporting (aligned to +/- strand):  ref base C (62/39);  new base G (2/4);  total (64/43)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

Only 100 of 112 total aligned reads displayed.

TACAGCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGT  >  NC_000913/2282685‑2282974
                                                                                                                                                   |                                                                                                                                               
tACAGCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCcg                                                                                                                                                >  4:372801/1‑149 (MQ=255)
    gCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCcg‑gcgg                                                                                                                                           >  3:538119/1‑149 (MQ=255)
    gCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCcg‑gcg                                                                                                                                            >  1:749326/1‑148 (MQ=255)
      gAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCgggg                                                                                                                                         >  4:418997/1‑149 (MQ=255)
          ttCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCCGGA‑GCGGGGtttt                                                                                                                                     >  4:329409/1‑149 (MQ=255)
          ttCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCCGGA‑GCGGGGttt                                                                                                                                      <  2:363014/148‑1 (MQ=255)
           tCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCCGGA‑GCGGGGttttt                                                                                                                                    <  4:664943/149‑1 (MQ=255)
           tCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGC‑CGGAGCGGGGttttt                                                                                                                                    <  3:329409/149‑1 (MQ=255)
              gACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGt                                                                                                                                 <  4:682676/149‑1 (MQ=255)
               aCCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTa                                                                                                                                >  4:846358/1‑149 (MQ=255)
                ccATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTAt                                                                                                                               <  3:418997/149‑1 (MQ=255)
                ccATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTAt                                                                                                                               >  3:63199/1‑149 (MQ=255)
                   tCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTg                                                                                                                            >  4:56558/1‑149 (MQ=255)
                    cGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTg                                                                                                                            >  1:526209/1‑148 (MQ=255)
                          tgACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGAtt                                                                                                                      >  1:679476/1‑148 (MQ=255)
                           gACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGGTTTTTGTATCTGCAGATTAt                                                                                                                    >  3:236296/1‑149 (MQ=255)
                           gACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGC‑CGGAGCGGGGtttt                                                                                                                                     <  3:824024/132‑1 (MQ=255)
                             cGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGc                                                                                                                  <  4:191305/149‑1 (MQ=255)
                                 aaaGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTg                                                                                                               >  2:759738/1‑148 (MQ=255)
                                 aaaGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGa                                                                                                              >  4:728555/1‑149 (MQ=255)
                                  aaGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCCCCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGa                                                                                                              <  2:844952/148‑1 (MQ=255)
                                    gAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTa                                                                                                           >  3:335276/1‑149 (MQ=255)
                                       aTCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACgg                                                                                                        >  3:128101/1‑149 (MQ=255)
                                        tCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACgg                                                                                                        >  1:130169/1‑148 (MQ=255)
                                         cAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTa                                                                                                      >  3:948533/1‑149 (MQ=255)
                                          aaaGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTAt                                                                                                     >  4:831049/1‑149 (MQ=255)
                                           aaGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTAt                                                                                                     <  1:50153/148‑1 (MQ=255)
                                           aaGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTAt                                                                                                     <  2:949779/148‑1 (MQ=255)
                                            aGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTAtt                                                                                                    <  2:749326/148‑1 (MQ=255)
                                              ttAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCt                                                                                                 >  4:504452/1‑149 (MQ=255)
                                               tAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCt                                                                                                 >  2:644289/1‑148 (MQ=255)
                                               tAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCt                                                                                                 >  2:783946/1‑148 (MQ=255)
                                                aaaGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTAt                                                                                               <  4:895029/149‑1 (MQ=255)
                                                 aaGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGttttt                                                                                                                                    <  1:128494/111‑1 (MQ=255)
                                                 aaGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGttttt                                                                                                                                    >  2:128494/1‑111 (MQ=255)
                                                  aGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTAttt                                                                                             >  3:390250/1‑149 (MQ=255)
                                                     tCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTAtttttt                                                                                          <  3:609435/149‑1 (MQ=255)
                                                     tCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTAttttt                                                                                           <  1:183793/148‑1 (MQ=255)
                                                      cAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTAtttttt                                                                                          >  2:938731/1‑148 (MQ=255)
                                                       aGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCcg‑gcg                                                                                                                                            >  4:725901/1‑97 (MQ=255)
                                                       aGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCcg‑gcg                                                                                                                                            <  3:725901/97‑1 (MQ=255)
                                                       aGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGc                                                                                        <  4:128101/149‑1 (MQ=255)
                                                        ggACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGc                                                                                        >  1:763564/1‑148 (MQ=255)
                                                           cGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAgg                                                                                     >  1:594443/1‑148 (MQ=255)
                                                                  cGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATaa                                                                             <  3:504452/149‑1 (MQ=255)
                                                                    tCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAgg                                                                           >  4:371091/1‑149 (MQ=255)
                                                                     cACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAgg                                                                           >  1:818825/1‑148 (MQ=255)
                                                                      aCCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCg                                                                         >  3:486734/1‑149 (MQ=255)
                                                                       cccGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGt                                                                        >  3:517026/1‑149 (MQ=255)
                                                                        ccGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGtt                                                                       >  3:779968/1‑149 (MQ=255)
                                                                            aCAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAcg                                                                   >  4:886753/1‑149 (MQ=255)
                                                                              aaaCCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCAcgcc                                                                 >  3:722727/1‑149 (MQ=255)
                                                                                      gCTGCAGCACATCGACTTCGTTCGCGCTTAATAGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTgg                                                         <  4:323955/149‑1 (MQ=255)
                                                                                        tGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAGAACCCCGCTCCG‑TCGGGGTTTTTTGTATCTGCAGATTATGCCTGACTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCa                                                       >  4:908569/1‑149 (MQ=255)
                                                                                            gcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAATCCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTc                                                    <  2:161251/148‑1 (MQ=255)
                                                                                            gcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCa                                                   >  4:459203/1‑149 (MQ=255)
                                                                                            gcaCATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCa                                                   <  3:372801/149‑1 (MQ=255)
                                                                                             cacaGCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAg                                                  >  4:850906/1‑149 (MQ=255)
                                                                                               caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGaa                                                <  4:486734/149‑1 (MQ=255)
                                                                                               caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGaa                                                <  3:846358/149‑1 (MQ=255)
                                                                                               caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGa                                                 <  2:379857/148‑1 (MQ=255)
                                                                                               caTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACACCGCATCTGGCATTCAGaa                                                <  3:789505/149‑1 (MQ=255)
                                                                                                aTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGaag                                               >  3:85748/1‑148 (MQ=255)
                                                                                                aTCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTAACGCCGCATCTGGCATTCAGaa                                                >  1:247758/1‑148 (MQ=255)
                                                                                                 tCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACCCCCCATCTGGCATTCAGAac                                               >  2:876324/1‑148 (MQ=255)
                                                                                                   gACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAacac                                             >  1:688152/1‑148 (MQ=255)
                                                                                                           tcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGa                                     <  1:783946/148‑1 (MQ=255)
                                                                                                           tcgCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGAt                                    >  3:495258/1‑149 (MQ=255)
                                                                                                               gcTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTaa                                 <  1:938731/148‑1 (MQ=255)
                                                                                                                  tAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTg                              >  2:591828/1‑148 (MQ=255)
                                                                                                                   aaTTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTgcc                            <  3:403931/149‑1 (MQ=255)
                                                                                                                   aaTTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTgc                             <  1:310494/148‑1 (MQ=255)
                                                                                                                      tGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTgccgc                          >  1:434124/1‑148 (MQ=255)
                                                                                                                       gCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCa                        >  3:754982/1‑149 (MQ=255)
                                                                                                                          gAATAAGTTGTAAAAAACCCCGCTCCG‑GCTGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGATCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCTGCCAGa                      <  2:50551/148‑1 (MQ=255)
                                                                                                                          gAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGa                      <  2:594443/148‑1 (MQ=255)
                                                                                                                           aaTAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGAc                     >  2:819976/1‑148 (MQ=255)
                                                                                                                           aaTAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACg                    >  4:151289/1‑149 (MQ=255)
                                                                                                                            aTAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAg                                                                                      <  4:315157/82‑1 (MQ=255)
                                                                                                                            aTAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAg                                                                                      >  3:315157/1‑82 (MQ=255)
                                                                                                                            aTAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGa                                                                                >  2:386207/1‑88 (MQ=255)
                                                                                                                            aTAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGa                                                                                <  1:386207/88‑1 (MQ=255)
                                                                                                                            aTAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACCCAAAACCGATTAATTGCCGCCAGACg                    <  2:818825/148‑1 (MQ=255)
                                                                                                                            aTAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACg                    <  2:141435/148‑1 (MQ=255)
                                                                                                                            aTAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACg                    >  1:170446/1‑148 (MQ=255)
                                                                                                                            aTAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGt                   >  4:79431/1‑149 (MQ=255)
                                                                                                                             tAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGt                   <  1:287378/148‑1 (MQ=255)
                                                                                                                             tAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTc                  >  3:844528/1‑149 (MQ=255)
                                                                                                                             tAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTc                  >  4:634311/1‑149 (MQ=255)
                                                                                                                                 ttGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATGCGGCGTGAACGCCTTATCTGGcc                                                       >  3:163327/1‑107 (MQ=38)
                                                                                                                                 ttGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATGCGGCGTGAACGCCTTATCTGGcc                                                       <  4:163327/108‑2 (MQ=38)
                                                                                                                                  tGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATGCGGCGTGAACGCCTTATCTGGc                                                        >  1:926866/1‑106 (MQ=38)
                                                                                                                                  tGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGgc              >  2:663899/1‑148 (MQ=255)
                                                                                                                                    tAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCt           >  4:917664/1‑149 (MQ=255)
                                                                                                                                     aaaaaaCCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTg          >  4:409824/1‑149 (MQ=255)
                                                                                                                                        aaaCCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATaa                                                                             <  2:857450/79‑1 (MQ=255)
                                                                                                                                        aaaCCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATaa                                                                             >  1:857450/1‑79 (MQ=255)
                                                                                                                                        aaaCCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCaa       >  4:763246/1‑149 (MQ=255)
                                                                                                                                             ccGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTgg   >  1:671594/1‑148 (MQ=255)
                                                                                                                                              cGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGt  >  1:694450/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
TACAGCGAAGTTCTGACCATCGTTGTTGACGGTAAAGAAATCAAAGTTAAAGCTCAGGACGTACAGCGTCACCCGTACAAACCGAAGCTGCAGCACATCGACTTCGTTCGCGCTTAATTGCTGAATAAGTTGTAAAAAACCCCGCTCCG‑GCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGT  >  NC_000913/2282685‑2282974

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