Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,428,954 G→A intergenic (‑192/‑282) rrsD ← / → yrdA 16S ribosomal RNA of rrnD operon/bacterial transferase hexapeptide domain protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,428,9540GA85.7% 177.1 / 24.0 77intergenic (‑192/‑282)rrsD/yrdA16S ribosomal RNA of rrnD operon/bacterial transferase hexapeptide domain protein
Reads supporting (aligned to +/- strand):  ref base G (5/6);  new base A (37/29);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 5.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

CCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAATTTTTTGCAC  >  NC_000913/3428857‑3429081
                                                                                                 |                                                                                                                               
ccGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGAtttttt                                                                                        >  1:1068344/1‑139 (MQ=11)
ccGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGACAACCCTGAATTTCAGGAtttttt                                                                                        >  1:306260/1‑139 (MQ=2)
  gTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTtc                                                                                      >  1:346298/1‑139 (MQ=17)
  gTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTtc                                                                                      >  2:1463679/1‑139 (MQ=17)
   tATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTtct                                                                                     >  1:1605612/1‑139 (MQ=17)
    aTCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTtctc                                                                                    >  2:62197/1‑139 (MQ=18)
       ttCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTc                                                                                 >  1:991568/1‑139 (MQ=25)
        tCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCg                                                                                >  1:1640218/1‑139 (MQ=37)
        tCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCg                                                                                <  1:62197/139‑1 (MQ=37)
        tCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCg                                                                                <  2:1605612/139‑1 (MQ=37)
               cccACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTcc                                                                         >  1:1235341/1‑139 (MQ=255)
               cccACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTcc                                                                         >  1:1341190/1‑139 (MQ=255)
                 cacacaGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGa                                                                       <  2:295422/139‑1 (MQ=255)
                 cacacaGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGa                                                                       <  2:566921/139‑1 (MQ=255)
                 cacacaGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGa                                                                       >  1:1293837/1‑139 (MQ=255)
                   cacaGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGAcc                                                                     >  1:775490/1‑139 (MQ=255)
                     caGATTGTCTGATAAATAGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCAt                                                                   >  2:1032596/1‑139 (MQ=255)
                        aTTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCt                                                                >  1:1002703/1‑139 (MQ=255)
                        aTTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCt                                                                <  2:1293837/139‑1 (MQ=255)
                        aTTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCt                                                                >  2:1200722/1‑139 (MQ=255)
                        aTTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCt                                                                >  1:24457/1‑139 (MQ=255)
                            tCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGa                                                            >  2:1592475/1‑139 (MQ=255)
                              tGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAg                                                          >  1:1492583/1‑139 (MQ=255)
                              tGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAg                                                          >  1:876336/1‑139 (MQ=255)
                              tGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAg                                                          >  1:1268614/1‑139 (MQ=255)
                                aTAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCgttgt                                                                            >  2:546762/1‑119 (MQ=255)
                                aTAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCgttgt                                                                            <  1:546762/119‑1 (MQ=255)
                                    aTTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTc                                                    <  1:1150250/139‑1 (MQ=255)
                                    aTTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTc                                                    <  2:1478976/139‑1 (MQ=255)
                                      tGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGAcc                                                                     <  1:782016/120‑1 (MQ=255)
                                      tGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGAcc                                                                     >  2:782016/1‑120 (MQ=255)
                                      tGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCAc                                                  >  2:422801/1‑139 (MQ=255)
                                        ttAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGt                                                >  2:440033/1‑139 (MQ=255)
                                        ttAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGt                                                <  1:1032596/139‑1 (MQ=255)
                                          aaaGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACTATCCTGTGTAGTGTTTCACGTTg                                              >  2:445424/1‑139 (MQ=255)
                                           aaGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTgt                                             <  2:1419775/139‑1 (MQ=255)
                                           aaGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTgt                                             <  1:422801/139‑1 (MQ=255)
                                            agagCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTgtc                                            <  1:1902/139‑1 (MQ=255)
                                                cAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATc                                                                  >  2:1107584/1‑113 (MQ=255)
                                                cAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATc                                                                  <  1:1107584/113‑1 (MQ=255)
                                                cAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGtct                                        >  2:1378830/1‑139 (MQ=255)
                                                  gTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCa                                      >  1:433082/1‑139 (MQ=255)
                                                  gTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGa                                                                       >  1:1032567/1‑106 (MQ=255)
                                                  gTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGa                                                                       <  2:1032567/106‑1 (MQ=255)
                                                  gTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGAcc                                                                     >  1:780348/1‑108 (MQ=255)
                                                  gTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGAcc                                                                     <  2:780348/108‑1 (MQ=255)
                                                  gTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCa                                      >  2:1206394/1‑139 (MQ=255)
                                                     gCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACg                                   >  2:842033/1‑139 (MQ=255)
                                                     gCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACg                                   >  1:593781/1‑139 (MQ=255)
                                                      cGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACgg                                  <  2:1260399/139‑1 (MQ=255)
                                                      cGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCTGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACgg                                  <  2:991568/139‑1 (MQ=255)
                                                        aCGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGAttttttt                                                                                       >  2:1115852/1‑84 (MQ=32)
                                                        aCGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGAttttttt                                                                                       <  1:1115852/84‑1 (MQ=32)
                                                        aCGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAAGCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTGCCGACCATc                                                                  >  1:1069755/1‑105 (MQ=255)
                                                        aCGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATc                                                                  <  2:1069755/105‑1 (MQ=255)
                                                        aCGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACg                                                 <  2:1083469/122‑1 (MQ=255)
                                                        aCGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACg                                                 >  1:1083469/1‑122 (MQ=255)
                                                        aCGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAg                                <  1:1378830/139‑1 (MQ=255)
                                                        aCGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAg                                <  2:346298/139‑1 (MQ=255)
                                                        aCGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAg                                <  1:445424/139‑1 (MQ=255)
                                                         cgcgGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAgg                               <  2:1341190/139‑1 (MQ=255)
                                                         cgcgGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAgg                               <  1:1206394/139‑1 (MQ=255)
                                                           cgGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGgcg                             <  2:1492583/139‑1 (MQ=255)
                                                           cgGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGgcg                             <  2:1640218/139‑1 (MQ=255)
                                                           cgGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGgcg                             <  2:801611/139‑1 (MQ=255)
                                                           cgGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGgcg                             <  1:1200722/139‑1 (MQ=255)
                                                            ggCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGgcgc                            <  2:1268614/139‑1 (MQ=255)
                                                             gCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCa                           >  1:1147410/1‑139 (MQ=255)
                                                                    gCTCACTGTCGCGAGGTGGCGTATATTACACTTTCCTCTCTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAgg                    <  2:1002703/139‑1 (MQ=255)
                                                                           gTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCa             >  1:1370940/1‑139 (MQ=255)
                                                                           gTCGCGAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCa             >  2:8491/1‑139 (MQ=255)
                                                                             cgcgAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGaa                                                           <  2:274263/91‑1 (MQ=255)
                                                                             cgcgAGGTGGCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGaa                                                           >  1:274263/1‑91 (MQ=255)
                                                                                      gCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAATTTTTTGcac  <  1:440033/139‑1 (MQ=255)
                                                                                      gCGTATATTACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAATTTTTTGcac  >  1:178754/1‑139 (MQ=255)
                                                                                             ttACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAgg                    >  1:385992/1‑114 (MQ=255)
                                                                                             ttACACTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAgg                    <  2:385992/114‑1 (MQ=255)
                                                                                                 |                                                                                                                               
CCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAATTTTTTGCAC  >  NC_000913/3428857‑3429081

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: