Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,451,404 T→G 6.1% D94A (GAC→GCC)  rplW ← 50S ribosomal subunit protein L23

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,451,4040TG 6.1% 313.6 / 7.9 109D94A (GAC→GCC) rplW50S ribosomal subunit protein L23
Reads supporting (aligned to +/- strand):  ref base T (49/53);  new base G (5/2);  total (54/55)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 111 total aligned reads displayed.

TTTTTTTCCAGCAACGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGCACAGCAGCT  >  NC_000913/3451255‑3451553
                                                                                                                                                     |                                                                                                                                                     
tttttttCCAGCAACGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGt                                                                                                                                                       >  2:309685/1‑150 (MQ=255)
 ttttttCCAGCAACGGAGCAAAAGGGTTGCCCTTGTGCAGCTCATGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACTACTGCCATTGTATTACTCCTCCGTCGTACTCCGCGCCGCCAACGAAGTc                                                                                                                                                      >  2:305176/1‑150 (MQ=255)
    tttCCAGCAACGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAg                                                                                                                                                   >  2:1467613/1‑150 (MQ=255)
    tttCCAGCAACGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAg                                                                                                                                                   <  2:1566763/150‑1 (MQ=255)
    tttCCAGCAACGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAg                                                                                                                                                   <  2:1607934/150‑1 (MQ=255)
      tCCAGCAACGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGAt                                                                                                                                                 >  1:1693417/1‑150 (MQ=255)
       ccAGCAACGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCATCCGACTTACTCAGCGCCGCCAACGAAGTCCAGAtt                                                                                                                                                >  2:929995/1‑150 (MQ=255)
        cAGCAACGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTc                                                                                                                                               >  1:1051438/1‑150 (MQ=255)
            aaCGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGc                                                                                                                                           <  1:919398/150‑1 (MQ=255)
             aCGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGcc                                                                                                                                          >  1:1221902/1‑150 (MQ=255)
              cGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCct                                                                                                                                         <  2:557033/150‑1 (MQ=255)
                   cAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCttt                                                                                                                                    <  1:754231/150‑1 (MQ=255)
                    aaaaGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTc                                                                                                                                   <  1:470378/150‑1 (MQ=255)
                    aaaaGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGGCCAGAATCTGGACTTCgttg                                                                                                                                   <  1:1733395/150‑2 (MQ=255)
                       cggTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCTACGACCGGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAgg                                                                                                                                <  2:851694/149‑1 (MQ=255)
                        ggTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAggg                                                                                                                               <  2:182914/150‑1 (MQ=255)
                        ggTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAggg                                                                                                                               <  2:795719/150‑1 (MQ=255)
                        ggTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAggg                                                                                                                               >  1:665946/1‑150 (MQ=255)
                          tttGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTg                                                                                                                             >  1:59823/1‑150 (MQ=255)
                             gCCCTTGTGCAGCTCAGGGTTAACGACTTTAACTACGGGGCGACGACCCGGAGATGTCGGTGTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAAGGAAGGCCAGATGCTGGCCTTCTTTCAGGGTGACg                                                                                                                          <  2:1220240/150‑1 (MQ=255)
                             gCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACg                                                                                                                          <  1:1714479/150‑1 (MQ=255)
                             gCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCAGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACg                                                                                                                          >  1:1052684/1‑150 (MQ=255)
                                    tgCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCtt                                                                                                                   >  1:122558/1‑150 (MQ=255)
                                        gCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCtttttt                                                                                                               >  1:939145/1‑150 (MQ=255)
                                        gCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCAGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCtttttt                                                                                                               >  1:1121072/1‑150 (MQ=255)
                                         cTCACGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGATTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGTCGCCAACGAACTCCAGATTCTGCACTTCTTTCAGGGTGACGTAAGCTTATTTc                                                                                                              >  2:503961/1‑150 (MQ=255)
                                          taatggTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTcc                                                                                                             <  2:131578/146‑1 (MQ=255)
                                          tCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTcc                                                                                                             <  1:788043/150‑1 (MQ=255)
                                           cAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCa                                                                                                            >  1:1611785/1‑150 (MQ=255)
                                           cAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCa                                                                                                            <  2:1628364/150‑1 (MQ=255)
                                            aGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAg                                                                                                           <  2:1614873/150‑1 (MQ=255)
                                            aGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCAAg                                                                                                           >  2:1535660/1‑150 (MQ=255)
                                             gggTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGt                                                                                                          <  1:461932/150‑1 (MQ=255)
                                               gTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCGGGGTGACGTAAGCTTTTTTCCAGTCg                                                                                                        >  1:1712714/1‑150 (MQ=255)
                                                 tAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCt                                                                                                      >  2:1602034/1‑150 (MQ=255)
                                                  aaCCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTa                                                                                                     >  1:520683/1‑150 (MQ=255)
                                                     cACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTacgc                                                                                                  >  1:1346408/1‑149 (MQ=255)
                                                               aCGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATAc                                                                                        >  1:1437261/1‑150 (MQ=255)
                                                                cGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACg                                                                                       >  2:1175512/1‑150 (MQ=255)
                                                                  tGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCt                                                                                     <  2:960528/150‑1 (MQ=255)
                                                                   ggCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGTCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTg                                                                                    >  2:816778/1‑150 (MQ=255)
                                                                     cgacgaCCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTc                                                                                  >  1:422772/1‑150 (MQ=255)
                                                                       acgacCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCg                                                                                <  2:59823/150‑1 (MQ=255)
                                                                         gacCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTg                                                                              <  1:592514/150‑1 (MQ=255)
                                                                          acCCGGAGATGTCGGTTTTCATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGa                                                                             <  1:495645/150‑1 (MQ=255)
                                                                            ccGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCAGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACg                                                                           <  2:3140/150‑1 (MQ=255)
                                                                             cGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGt                                                                          <  1:882082/150‑1 (MQ=255)
                                                                              gcagaTGTCGGATTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGtt                                                                         <  2:689764/148‑1 (MQ=255)
                                                                                        gTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATACTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTccc                                                               >  2:547654/1‑150 (MQ=255)
                                                                                         tttACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGGCGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCt                                                              <  1:429237/150‑1 (MQ=255)
                                                                                         tttACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCt                                                              >  1:1698683/1‑150 (MQ=255)
                                                                                           tACATTTAACAACTTCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCttt                                                            <  1:1229553/150‑1 (MQ=255)
                                                                                           tACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCttt                                                            <  2:1069071/150‑1 (MQ=255)
                                                                                           tACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCttt                                                            >  1:680768/1‑150 (MQ=255)
                                                                                           tACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCttt                                                            <  1:888276/150‑1 (MQ=255)
                                                                                              aTTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAAc                                                         >  2:446974/1‑150 (MQ=255)
                                                                                               tttAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACt                                                        >  1:1301185/1‑150 (MQ=255)
                                                                                                  aacaacTGCCATTGTATTACTCCTCCGCCTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCATCTTTCAGGGTGACGTAAGCTTTTATCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTccc                                                     >  1:77029/1‑150 (MQ=255)
                                                                                                  aacaacTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTAcc                                                     >  2:753219/1‑150 (MQ=255)
                                                                                                   acaacTGCCTTTGTATTACTCCTCCGACTTACTCTGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCa                                                    <  2:172780/150‑1 (MQ=255)
                                                                                                         gCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGgtgt                                              >  2:745073/1‑150 (MQ=255)
                                                                                                           cATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTa                                            >  1:658124/1‑150 (MQ=255)
                                                                                                           cATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTa                                            <  2:192317/150‑1 (MQ=255)
                                                                                                           cATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTa                                            >  1:1366982/1‑150 (MQ=255)
                                                                                                               gTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcga                                        >  2:632489/1‑150 (MQ=255)
                                                                                                                tATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgac                                       <  1:1124867/150‑1 (MQ=255)
                                                                                                                   tACTCCTCCGACTTACTCACCGCCGCCAACGAAGTCCAGATTCTGGCATTCTGTCAGGGTGACGAAAGCTTTATTCCAGTCGCTACGACGACCGATACGCTCTCAGTGACGTTTAACTTTCCAATTAACTACCAGGGTGTTAAcgacttc                                    >  1:671630/1‑150 (MQ=255)
                                                                                                                    aCTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcg                                   <  2:1475740/150‑1 (MQ=255)
                                                                                                                    aCTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcg                                   <  2:1249307/150‑1 (MQ=255)
                                                                                                                      tcctccGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgac                                 <  2:651149/150‑1 (MQ=255)
                                                                                                                        ctccGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgactt                               >  2:734636/1‑150 (MQ=255)
                                                                                                                          ccGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGTGTGTTAAcgacttcgacttcg                             <  2:1260007/150‑1 (MQ=255)
                                                                                                                            gACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgacttcgac                           <  2:1052684/150‑1 (MQ=255)
                                                                                                                              cTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgacttcgactt                         <  2:1303696/150‑1 (MQ=255)
                                                                                                                               ttACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTACCTATCCCTTTAACTTCCAGGGTGTTAACgacttcgacttcgacttc                        <  2:1049096/150‑1 (MQ=255)
                                                                                                                               ttACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACgacttcgacttcgacttc                        <  2:298960/150‑1 (MQ=255)
                                                                                                                                tACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCa                       >  1:521517/1‑150 (MQ=255)
                                                                                                                                 aCTCCGCGCCTCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTACGCTTTTTTCCCGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCaa                      <  1:131735/150‑1 (MQ=255)
                                                                                                                                    cAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGCCGTTTAACTTTCCCTTTAACTACCAGGGTGTTACCGCCTTCGACTTCGACTTCAAACa                   >  1:1047546/1‑150 (MQ=255)
                                                                                                                                    cAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACa                   <  2:122558/150‑1 (MQ=255)
                                                                                                                                     aGCGCCGCCAACGAATTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAg                  <  1:1210099/150‑1 (MQ=255)
                                                                                                                                      gcgcCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGt                 <  1:1602034/150‑1 (MQ=255)
                                                                                                                                       ggccgccAACGAAGTCGAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAAGTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGtt                <  2:360838/149‑1 (MQ=255)
                                                                                                                                        gccgccAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTGACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGttt               >  1:1411971/1‑150 (MQ=255)
                                                                                                                                         ccgccAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTc              <  1:753219/150‑1 (MQ=255)
                                                                                                                                         ccgccAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTc              <  1:1467613/150‑1 (MQ=255)
                                                                                                                                          cgccAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCt             >  1:748077/1‑150 (MQ=255)
                                                                                                                                          cctgaaagAAGGCCAGAATCTGGACTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCt             >  1:1528758/8‑150 (MQ=255)
                                                                                                                                           gccAACGAAGTCCAGATTCTGGCCTTCTTTCCGGGGGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTg            <  1:1643957/150‑1 (MQ=255)
                                                                                                                                           gccAACGAAGTCAAGATTCTGGACTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGATACGACGACCGATACGATGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCCACTTCGACTTAAAACAGTTTCTg            >  2:1203698/1‑150 (MQ=255)
                                                                                                                                            ccAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGc           <  2:1693417/150‑1 (MQ=255)
                                                                                                                                             gaaagAAGGCCAGAATCTGGACTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGca          >  1:979014/5‑150 (MQ=255)
                                                                                                                                              aaagAAGGCCAGAATCTGGACTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGcac         >  2:889412/4‑150 (MQ=255)
                                                                                                                                               aCGAAGTCCAGATTCTGGCCGTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGCCGCCCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGgaca        >  2:1210421/1‑150 (MQ=255)
                                                                                                                                                cGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTATAACTACCAGTGTGTTAACGACTTCGACATCGACTTCAAACAGTTTCTGCAcag       >  1:1104361/1‑150 (MQ=255)
                                                                                                                                                agAAGGCCAGAATCTGGACTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGCAcag       >  2:212257/2‑150 (MQ=255)
                                                                                                                                                 gAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGGCCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTAGACTTCGACTTCAAACAGTTTCTGCAcagc      >  1:914575/1‑150 (MQ=255)
                                                                                                                                                  aaGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGCCGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGCAcagca     >  1:1303519/1‑150 (MQ=255)
                                                                                                                                                  aaGGCCAGAATCTGGACTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGCAcagca     >  2:665648/1‑150 (MQ=255)
                                                                                                                                                     tCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGAGTTAACGACTTCGACTTCTACTTCAAACAGTTTCTGCACAGCAGCt  >  1:1726943/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTTTTTCCAGCAACGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGCACAGCAGCT  >  NC_000913/3451255‑3451553

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: