Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,912,482 Δ2 bp coding (94‑95/882 nt) htpX ← putative endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,912,4820A.98.3% 280.1 / ‑2.1 60coding (95/882 nt)htpXputative endopeptidase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (39/20);  total (39/21)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
*NC_0009131,912,4830C.98.3% 283.9 / ‑3.1 60coding (94/882 nt)htpXputative endopeptidase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (39/20);  total (39/21)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGAACGCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGCATCATAATTTTCTTTTTACCTCAGTTTAACAAAACGCAATATGCGAT  >  NC_000913/1912342‑1912618
                                                                                                                                            ||                                                                                                                                       
ccAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACt                                                                                                                              <  1:713152/151‑1 (MQ=255)
 cAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTg                                                                                                                             <  2:1195036/151‑1 (MQ=255)
  aGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGt                                                                                                                            >  2:102199/1‑151 (MQ=255)
  aGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGt                                                                                                                            >  2:756244/1‑151 (MQ=255)
     cAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATc                                                                                                                         >  1:1145201/1‑151 (MQ=255)
      aaCGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATcc                                                                                                                        >  2:275861/1‑151 (MQ=255)
      aaCGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATcc                                                                                                                        <  2:1253647/151‑1 (MQ=255)
      aaCGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCACTGA‑‑GCTGCTCGACTGTATcc                                                                                                                        >  2:876117/1‑151 (MQ=255)
            cccTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTc                                                                                                                   >  2:432395/1‑150 (MQ=255)
            cccTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTc                                                                                                                   >  1:925365/1‑150 (MQ=255)
                ttCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCt                                                                                                              >  2:743012/1‑151 (MQ=255)
                  cGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTc                                                                                                             >  1:417750/1‑150 (MQ=255)
                  cGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTc                                                                                                             >  2:833058/1‑150 (MQ=255)
                    ttACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGt                                                                                                          >  1:368800/1‑151 (MQ=255)
                              gCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCg                                                                                                 >  2:27245/1‑150 (MQ=255)
                                 cGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGaaa                                                                                              >  2:10573/1‑150 (MQ=255)
                                   aTCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACg                                                                                           >  2:1090071/1‑151 (MQ=255)
                                    tCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACg                                                                                           >  2:1017514/1‑150 (MQ=255)
                                         tCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCAtt                                                                                     >  1:951635/1‑151 (MQ=255)
                                         tCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCAt                                                                                      <  2:1121366/150‑1 (MQ=255)
                                           cccGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTAc                                                                                   >  2:33713/1‑151 (MQ=255)
                                               ccAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGcc                                                                               >  2:1079754/1‑151 (MQ=255)
                                               ccAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGcc                                                                               >  2:846330/1‑151 (MQ=255)
                                                    aGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGtt                                                                          <  2:304516/151‑1 (MQ=255)
                                                     gATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTc                                                                         >  2:886674/1‑151 (MQ=255)
                                                      aTCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCg                                                                        >  2:359328/1‑151 (MQ=255)
                                                          tAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGtt                                                                      >  2:284912/1‑149 (MQ=255)
                                                           aaTGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGc                                                                   >  2:503439/1‑151 (MQ=255)
                                                              gCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAgg                                                                <  2:639302/151‑1 (MQ=255)
                                                               ccatccatTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAgg                                                                >  1:657886/1‑150 (MQ=255)
                                                                        tGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGcgcg                                                       <  2:335545/150‑1 (MQ=255)
                                                                           aCATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGAtt                                                    >  1:1164887/1‑150 (MQ=255)
                                                                             aTCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGc                                                 >  1:1068550/1‑151 (MQ=255)
                                                                                aGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGcatc                                              <  1:794051/151‑1 (MQ=255)
                                                                                aGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGcatc                                              <  1:210988/151‑1 (MQ=255)
                                                                                aGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGcatc                                              >  2:796063/1‑151 (MQ=255)
                                                                                  aaGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGcatcat                                            <  1:1039739/151‑1 (MQ=255)
                                                                                  aaGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGcatcat                                            <  2:478138/151‑1 (MQ=255)
                                                                                       aaaCGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGCATCATAAttt                                       >  2:793825/1‑151 (MQ=255)
                                                                                        aaCGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGCATCATAAtttt                                      >  2:424132/1‑151 (MQ=255)
                                                                                         aCGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGCATCATAAtttt                                      >  2:1243595/1‑150 (MQ=255)
                                                                                           gAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGCATCATAATTTTCtt                                   <  2:177976/151‑1 (MQ=255)
                                                                                            aaGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGCATCATAATTTTCttt                                  >  1:344281/1‑151 (MQ=255)
                                                                                            aaGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGCATCATAATTTTCtt                                   <  2:509379/150‑1 (MQ=255)
                                                                                             aGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGCATCATAATTTTCtttt                                 <  1:783464/151‑1 (MQ=255)
                                                                                              ggAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGCATCATAATTTTCttttt                                >  2:814251/1‑151 (MQ=255)
                                                                                                 accaccAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGCATCATAATTTTCTTTTTAcc                             <  1:549265/151‑1 (MQ=255)
                                                                                                  ccaccaAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGCATCATAATTTTCTTTTTACCt                            >  1:1035761/1‑151 (MQ=255)
                                                                                                    accaAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGCATCATAATTTTCTTTTTACCTc                           >  1:898882/1‑150 (MQ=255)
                                                                                                      caAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGCATCATAATTTTCTTTTTACCTCAGt                        >  1:28739/1‑151 (MQ=255)
                                                                                                        aaaCCGAACAGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGCATCATAATTTTCTTTTTACCTCAGttt                      >  2:876563/1‑151 (MQ=255)
                                                                                                              aaCAGCAAGGCCATGATCATCAGCCCCGGAACGGTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCACGAAATCGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGCATCATAATTTTCTTTTTACCTCAGTTTAACa                  <  1:1136738/151‑1 (MQ=255)
                                                                                                                 aGCAAGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGCATCATAATTTTCTTTTTACCTCAGTTTAACAAAACg             <  1:703466/151‑1 (MQ=255)
                                                                                                                     aGGCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGCATCATAATTTTCTTTTTACCTCAGTTTAACAAAACGCAat         >  1:699290/1‑151 (MQ=255)
                                                                                                                      ggCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGCATCATAATTTTCTTTTTACCTCAGTTTAACAAAACGCAata        <  2:786209/151‑1 (MQ=255)
                                                                                                                       gCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGCATCATAATTTTCTTTTTACCTCAGTTTAACAAAACGCAatat       <  2:1019268/151‑1 (MQ=255)
                                                                                                                       gCCATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGCATCATAATTTTCTTTTTACCTCAGTTTAACAAAACGCAata        >  2:844780/1‑150 (MQ=255)
                                                                                                                        ccATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGCATCATAATTTTCTTTTTACCTCAGTTTAACAAAACGCAATATg      <  1:404419/151‑1 (MQ=255)
                                                                                                                        ccATGATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGCATCATAATTTTCTTTTTACCTCAGTTTAACAAAACGCAATATg      >  2:164149/1‑151 (MQ=255)
                                                                                                                            gATCATCAGCCCCTGA‑‑GCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGCATCATAATTTTCTTTTTACCTCAGTTTAACAAAACGCAATATGCGAt  <  1:868170/151‑1 (MQ=255)
                                                                                                                                            ||                                                                                                                                       
CCAGCCAACGTTCCCTTTCGTTACGCGGTTGCTCGATCACTTCCCCGCCAACAGATCGTAATGCCATCCATTTGGACATCAGAAGCGAAACGAAGGAACCACCAAAACCGAACAGCAAGGCCATGATCATCAGCCCCTGAACGCTGCTCGACTGTATCCCTGTCAGGCTCAGTACCAGCCCGAAAACGACCATTACGGCCAGGTTCGTTAGCAGGAAGAGCGCGATTCGCATCATAATTTTCTTTTTACCTCAGTTTAACAAAACGCAATATGCGAT  >  NC_000913/1912342‑1912618

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