Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,451,416 C→A 6.6% G90V (GGC→GTC)  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,4160CA 6.6% 274.1 / 8.0 91G90V (GGC→GTC) rplW50S ribosomal subunit protein L23
Reads supporting (aligned to +/- strand):  ref base C (39/46);  new base A (4/2);  total (43/48)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-01

ACGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGCACAGCAGCTTTGAT  >  NC_000913/3451268‑3451558
                                                                                                                                                    |                                                                                                                                              
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                                                                                                                                 <  1:669494/149‑1 (MQ=255)
               tGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGa                                                                                                                                 <  2:865066/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:499328/1‑149 (MQ=255)
                          aGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCtttt                                                                                                                      >  1:1247130/1‑149 (MQ=255)
                          aGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCtttt                                                                                                                      <  2:904515/149‑1 (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:1332847/149‑1 (MQ=255)
                                      aCCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTa                                                                                                          <  2:1335668/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                                                                                 <  2:439048/149‑1 (MQ=255)
                                                               ccGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGAc                                                                                 >  1:1218447/1‑149 (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:38959/1‑149 (MQ=255)
                                                                           gTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTcc                                                                     <  1:929589/149‑1 (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:737759/1‑149 (MQ=255)
                                                                                 aTTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTaa                                                               >  1:749553/1‑149 (MQ=255)
                                                                                 aTTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTaa                                                               >  2:1242741/1‑149 (MQ=255)
                                                                                 aTTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTaa                                                               <  1:493305/149‑1 (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:352051/1‑149 (MQ=255)
                                                                                                                      tCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCaaa                          <  2:498931/149‑1 (MQ=255)
                                                                                                                       cAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAAc                         <  1:606303/149‑1 (MQ=255)
                                                                                                                              ctgaaagAAGGCCAGAATCTGGACTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCt                  >  1:1403678/7‑149 (MQ=255)
                                                                                                                               ccAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTg                 >  2:422585/1‑149 (MQ=255)
                                                                                                                               ccAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTg                 >  1:256327/1‑149 (MQ=255)
                                                                                                                               ccAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTg                 <  2:749553/149‑1 (MQ=255)
                                                                                                                                cAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGc                >  2:700629/1‑149 (MQ=255)
                                                                                                                                 aaCGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGCCGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGca               >  1:877704/1‑149 (MQ=255)
                                                                                                                                  aagAAGGCCAGAATCTGGACTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGcac              >  2:343361/3‑149 (MQ=255)
                                                                                                                                  aagAAGGCCAGAATCTGGACTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGcac              >  1:1584940/3‑149 (MQ=255)
                                                                                                                                   agAAGGCCAGAATCTGGACTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGcaca             >  1:1414324/2‑149 (MQ=255)
                                                                                                                                    gAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGCAcag            <  2:1517761/149‑1 (MQ=255)
                                                                                                                                          cAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGCACAGCAGCtt      <  1:828238/149‑1 (MQ=255)
                                                                                                                                           aGATTCTGGCCTTCTTTCAGGGTGACGTAATCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGCACTGCAGCttt     <  1:929016/149‑1 (MQ=255)
                                                                                                                                              ttCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGCACAGCAGCTTTGAt  <  2:159943/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
ACGGAGCAAAAGGTTTGCCCTTGTGCAGCTCAGGGTTAACCACTTTAACTACGTGGCGACGACCCGGAGATGTCGGTTTACATTTAACAACTGCCATTGTATTACTCCTCCGACTTACTCAGCGCCGCCAACGAAGTCCAGATTCTGGCCTTCTTTCAGGGTGACGTAAGCTTTTTTCCAGTCGCTACGACGACCGATACGCTGTCCGTGACGTTTAACTTTCCCTTTAACTACCAGGGTGTTAACGACTTCGACTTCGACTTCAAACAGTTTCTGCACAGCAGCTTTGAT  >  NC_000913/3451268‑3451558

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: