Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,451,404 T→G 6.4% 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.4% 282.6 / 9.8 94D94A (GAC→GCC) rplW50S ribosomal subunit protein L23
Reads supporting (aligned to +/- strand):  ref base T (43/45);  new base G (4/2);  total (47/47)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

TTTTTTCCAGCAACGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGCACAG  >  NC_000913/3451256‑3451548
                                                                                                                                                    |                                                                                                                                                
ttttttCCAGCAACGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGt                                                                                                                                                  >  1:498931/1‑149 (MQ=255)
   tttCCAGCAACGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCa                                                                                                                                               >  1:1080350/1‑149 (MQ=255)
        aGCAACGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTc                                                                                                                                          >  1:233383/1‑149 (MQ=255)
        aGCAACGGAGCAAAAGGTTTGCCCTTGGGCAGCTCAGGGTTAACCAATTTAACTACGTGGCGACGCCCCGGAGATGTCGGGTTACATTTAACAACTGCCATTGTATTTCTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTc                                                                                                                                          >  2:828238/1‑149 (MQ=255)
         gCAACGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGAAAACCACTTTAACTACCTGGCGACGACCCGGAGAAGTCGGTTTACATTCAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGACGCCAACGAAGTCCAGATTCt                                                                                                                                         >  2:1135496/1‑149 (MQ=255)
           aaCGGAGCAAATGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTgg                                                                                                                                       <  2:943449/149‑1 (MQ=255)
            aCGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGc                                                                                                                                      <  1:925299/149‑1 (MQ=255)
               gAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTACTCCGACTTACACAGCGCCGCCAACGAAGTCCAGATTCTGGCctt                                                                                                                                   >  2:929016/1‑149 (MQ=255)
                  cAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGGCCAGAATCTGGActtcgt                                                                                                                                <  2:1428237/149‑3 (MQ=255)
                   aaaaGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCTCCGCCAACGAAGTCCAGAATCTGGCCTTCttt                                                                                                                               >  1:1001858/1‑149 (MQ=255)
                   aaaaGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCttt                                                                                                                               >  1:499191/1‑149 (MQ=255)
                    aaaGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGCTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGGCCAGAATCTGGACTTCgttg                                                                                                                              <  1:1290511/149‑2 (MQ=255)
                      aGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAg                                                                                                                            <  1:1497484/149‑1 (MQ=255)
                         tttGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGt                                                                                                                         >  1:439048/1‑149 (MQ=255)
                          ttGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTg                                                                                                                        <  2:1368359/149‑1 (MQ=255)
                           tGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGa                                                                                                                       <  2:865066/149‑1 (MQ=255)
                           tGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGa                                                                                                                       <  1:669494/149‑1 (MQ=255)
                                ttGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTaa                                                                                                                  <  1:348621/149‑1 (MQ=255)
                                 tgtgCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAg                                                                                                                 <  2:630696/149‑1 (MQ=255)
                                 tgtgCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGGCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAg                                                                                                                 <  1:1131756/149‑1 (MQ=255)
                                      aGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCtttt                                                                                                            <  2:904515/149‑1 (MQ=255)
                                      aGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCtttt                                                                                                            >  1:1247130/1‑149 (MQ=255)
                                      aGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCtttt                                                                                                            >  2:499328/1‑149 (MQ=255)
                                       gCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCttttt                                                                                                           <  1:1032415/149‑1 (MQ=255)
                                       gCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCttttt                                                                                                           <  1:306559/149‑1 (MQ=255)
                                        cTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCtttttt                                                                                                          >  2:1040083/1‑149 (MQ=255)
                                        cTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATGCTGGCCTTCGTTCAGGGGGACGTAAGCtttttt                                                                                                          >  1:1071000/1‑149 (MQ=255)
                                           aGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCa                                                                                                       <  2:1345146/149‑1 (MQ=255)
                                             ggTTTAGCAATTTTACTACGTGGCGTCGACCCGGTGTTGTCGGTTTACATTTATCAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGt                                                                                                     <  2:113533/149‑1 (MQ=255)
                                                  aCCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTa                                                                                                >  1:1517761/1‑149 (MQ=255)
                                                  aCCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTa                                                                                                <  2:1335668/149‑1 (MQ=255)
                                                  aCCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTa                                                                                                <  2:1332847/149‑1 (MQ=255)
                                                   ccACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACACCTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCAGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTac                                                                                               >  2:1083693/1‑149 (MQ=255)
                                                   ccACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTac                                                                                               >  2:365838/1‑149 (MQ=255)
                                                     aCTTTAACTACGTGGCGACGACCCGGAGACGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTacga                                                                                             >  2:606303/1‑149 (MQ=255)
                                                           aCTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCg                                                                                       >  2:936398/1‑149 (MQ=255)
                                                            cTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGa                                                                                      >  2:731324/1‑149 (MQ=255)
                                                             tACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGccgaccgaa                                                                                     >  2:333702/1‑148 (MQ=255)
                                                                 tGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGc                                                                                 >  2:577805/1‑149 (MQ=255)
                                                                  ggCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTTGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCt                                                                                <  2:800952/149‑1 (MQ=255)
                                                                   gCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTg                                                                               >  2:1182592/1‑149 (MQ=255)
                                                                    cgacgaCCCGGAGATGTCGGTTTACATTTAACAACTGCAATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTTAGCTTTTTTCCAGTCGATACGACGACCGATACGCTGt                                                                              <  2:701406/149‑1 (MQ=255)
                                                                        gacCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGt                                                                          <  2:1041507/149‑1 (MQ=255)
                                                                           ccGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGCCGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGAc                                                                       >  1:1433588/1‑149 (MQ=255)
                                                                           ccGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGAc                                                                       >  1:1218447/1‑149 (MQ=255)
                                                                           ccGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGAc                                                                       <  2:439048/149‑1 (MQ=255)
                                                                             ggAGATGTGGGTTTACATTTATAATCTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGACGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGt                                                                     <  1:1004362/149‑1 (MQ=255)
                                                                              gagaTGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGtt                                                                    >  1:668392/1‑149 (MQ=255)
                                                                                  tGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGCCTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAAc                                                                <  1:306453/149‑1 (MQ=255)
                                                                                  tGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAAc                                                                >  1:1592755/1‑149 (MQ=255)
                                                                                    tCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACtt                                                              >  1:1045373/1‑149 (MQ=255)
                                                                                     cGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGGGACGTTTAACttt                                                             <  2:469904/149‑1 (MQ=255)
                                                                                     cGGTTTACATTCAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCTGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACttt                                                             <  1:1135496/149‑1 (MQ=255)
                                                                                       gTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTcc                                                           <  1:929589/149‑1 (MQ=255)
                                                                                       gTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTcc                                                           >  1:38959/1‑149 (MQ=255)
                                                                                        tttACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTccc                                                          <  1:692677/149‑1 (MQ=255)
                                                                                        tttACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTccc                                                          <  1:1083693/149‑1 (MQ=255)
                                                                                        tttACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCACCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTccc                                                          <  2:1327436/149‑1 (MQ=255)
                                                                                           aCATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCttt                                                       <  2:1189391/149‑1 (MQ=255)
                                                                                             aTTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTaa                                                     >  2:1242741/1‑149 (MQ=255)
                                                                                             aTTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTaa                                                     >  2:737759/1‑149 (MQ=255)
                                                                                             aTTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTaa                                                     <  1:493305/149‑1 (MQ=255)
                                                                                             aTTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTaa                                                     >  1:749553/1‑149 (MQ=255)
                                                                                             aTTTAACAACTGCCATTGTATTACTCCTCAGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTaa                                                     <  1:1185310/149‑1 (MQ=255)
                                                                                                    aacTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAg                                              <  2:959219/149‑1 (MQ=255)
                                                                                                     acTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAgg                                             <  2:1080350/149‑1 (MQ=255)
                                                                                                      cTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAggg                                            >  2:617422/1‑149 (MQ=255)
                                                                                                       tGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGgt                                           <  1:1137283/149‑1 (MQ=255)
                                                                                                         ccATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTGTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGgtgt                                         >  2:526455/1‑149 (MQ=255)
                                                                                                          cATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGtt                                        >  1:121824/1‑149 (MQ=255)
                                                                                                          cATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGtt                                        >  1:1056364/1‑149 (MQ=255)
                                                                                                             tGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAc                                     <  2:233383/149‑1 (MQ=255)
                                                                                                               tATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcga                                   <  2:1247130/149‑1 (MQ=255)
                                                                                                                aTTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTATTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTTCCAGGGTGTTAAcgac                                  >  2:1148988/1‑149 (MQ=255)
                                                                                                                       ctccGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgact                           >  1:235029/1‑149 (MQ=255)
                                                                                                                         ccGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgacttc                         <  2:499191/149‑1 (MQ=255)
                                                                                                                             cTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgacttcgact                     <  1:417843/149‑1 (MQ=255)
                                                                                                                              ttccTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCGTTTATCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgacttcgactt                    <  1:333702/146‑1 (MQ=255)
                                                                                                                              ttACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAAcgacttcgacttcgactt                    <  1:406196/149‑1 (MQ=255)
                                                                                                                               tACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACgacttcgacttcgacttc                   <  2:915100/149‑1 (MQ=255)
                                                                                                                                aCTCAGCGCCGCCAACGAAGTCCAGCTTCTGGCCTTCTTTCAGGGTGAAGTAAGCTTTTTTCCCGTCGCTACGACGACCGATAAGATGTCCGTGACGTTTAACTTTCCCTTTACCTACCAGGGTGTTAAcgacttcgacttcgccttaa                  >  1:866316/1‑147 (MQ=255)
                                                                                                                                  tCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGCCCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCaaa                >  2:1456733/1‑149 (MQ=255)
                                                                                                                                  tCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCaaa                <  2:498931/149‑1 (MQ=255)
                                                                                                                                  tCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCaaa                >  2:352051/1‑149 (MQ=255)
                                                                                                                                   cAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAAc               <  1:606303/149‑1 (MQ=255)
                                                                                                                                          ctgaaagAAGGCCAGAATCTGGACTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCt        >  1:1403678/7‑149 (MQ=255)
                                                                                                                                           ccAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTg       >  1:256327/1‑149 (MQ=255)
                                                                                                                                           ccAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTg       <  2:749553/149‑1 (MQ=255)
                                                                                                                                           ccAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTg       >  2:422585/1‑149 (MQ=255)
                                                                                                                                            cAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGc      >  2:700629/1‑149 (MQ=255)
                                                                                                                                             aaCGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGCCGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGca     >  1:877704/1‑149 (MQ=255)
                                                                                                                                              aagAAGGCCAGAATCTGGACTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGcac    >  1:1584940/3‑149 (MQ=255)
                                                                                                                                              aagAAGGCCAGAATCTGGACTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGcac    >  2:343361/3‑149 (MQ=255)
                                                                                                                                               agAAGGCCAGAATCTGGACTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGcaca   >  1:1414324/2‑149 (MQ=255)
                                                                                                                                                gAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGCAcag  <  2:1517761/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TTTTTTCCAGCAACGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGCACAG  >  NC_000913/3451256‑3451548

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