Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 428,067 A→G T141A (ACC→GCC)  secD → SecYEG protein translocase auxillary subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913428,0670AG99.2% 417.8 / ‑4.5 117T141A (ACC→GCC) secDSecYEG protein translocase auxillary subunit
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (64/52);  total (64/53)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

Only 100 of 117 total aligned reads displayed.

ATGGGCGTTATGGGTGACAAATACGTCGTGGCGCTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACACCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGTGATGAAGCC  >  NC_000913/427920‑428213
                                                                                                                                                   |                                                                                                                                                  
aTGGGCGTTATGGGTGACAAATACGTCGTGGCGCTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCcg                                                                                                                                                 >  4:939307/1‑151 (MQ=255)
   ggCGTTATGGGTGACAAATACGTCGTGGCGCTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCgcgc                                                                                                                                              >  3:131930/1‑151 (MQ=255)
    gCGTTATGGGTGACAAATACGTCGTGGCGCTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCt                                                                                                                                             >  1:34557/1‑151 (MQ=255)
     cGTTATGGGTGACAAATACGTCGTGGCGCTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCtt                                                                                                                                            >  4:1659877/1‑151 (MQ=255)
         aTGGGTGACAAATACGTCGTGGCGCTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGc                                                                                                                                         >  4:1478883/1‑150 (MQ=255)
           gggTGACAAATACGTCGTGGCGCTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCaaa                                                                                                                                      <  4:491701/151‑1 (MQ=255)
            ggTGACAAATACGTCGTGGCGCTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAAc                                                                                                                                     >  4:764349/1‑151 (MQ=255)
             gTGACAAATACGTCGTGGCGCTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACt                                                                                                                                    <  3:262173/151‑1 (MQ=255)
             gTGACAAATACGTCGTGGCGCTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACt                                                                                                                                    >  2:109775/1‑151 (MQ=255)
               gACAAATACGTCGTGGCGCTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTcc                                                                                                                                  >  2:601310/1‑151 (MQ=255)
               gACAAATACGTCGTGGCGCTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTc                                                                                                                                   <  2:365510/150‑1 (MQ=255)
                aCAAATACGTCGTGGCGCTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTcc                                                                                                                                  >  3:148882/1‑150 (MQ=255)
                aCAAATACGTCGTGGCGCTTAACCTTGCCCCAGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGCCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCTAGCTTGGCCAACt                                                                                                                                    >  1:241070/1‑148 (MQ=255)
                   aaTACGTCGTGGCGCTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGa                                                                                                                              <  1:113609/151‑1 (MQ=255)
                    aTACGTCGTGGCGCTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGaa                                                                                                                             <  2:1047440/151‑1 (MQ=255)
                    aTACGTCGTGGCGCTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGa                                                                                                                              >  4:567910/1‑150 (MQ=255)
                     tACGTCGTGGCGCTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAAc                                                                                                                            >  1:437164/1‑151 (MQ=255)
                      aCGTCGTGGCGCTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACa                                                                                                                           <  2:1573907/151‑1 (MQ=255)
                       cgtcgtGGCGCTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACa                                                                                                                           >  2:544754/1‑150 (MQ=255)
                        gtcgtGGCGCTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACaaa                                                                                                                         <  4:22723/151‑1 (MQ=255)
                        gtcgtGGCGCTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACaaa                                                                                                                         <  4:798132/151‑1 (MQ=255)
                          cgtGGCGCTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAat                                                                                                                       <  4:1188462/151‑1 (MQ=255)
                           gtGGCGCTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAat                                                                                                                       <  4:1289926/150‑1 (MQ=255)
                            tGGCGCTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAatat                                                                                                                     >  3:265187/1‑151 (MQ=255)
                             ggCGCTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATc                                                                                                                    <  2:539160/151‑1 (MQ=255)
                               cgcTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGa                                                                                                                  <  1:471792/151‑1 (MQ=255)
                                gcTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTTACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCACACTCCAGGAACAAAATATCGAt                                                                                                                 <  1:118285/151‑1 (MQ=255)
                                gcTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGAt                                                                                                                 <  4:686253/151‑1 (MQ=255)
                                    aaCCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGcc                                                                                                             >  4:341784/1‑151 (MQ=255)
                                     aCCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCt                                                                                                            <  4:1483554/151‑1 (MQ=255)
                                      ccTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCt                                                                                                            >  2:1081386/1‑150 (MQ=255)
                                         tGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTAc                                                                                                          <  3:61975/149‑1 (MQ=255)
                                         tGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTAc                                                                                                          <  2:15104/149‑1 (MQ=255)
                                         tGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGc                                                                                                        >  2:1645182/1‑151 (MQ=255)
                                          gCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGc                                                                                                        >  4:1310040/1‑150 (MQ=255)
                                          gCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGc                                                                                                        >  4:1422313/1‑150 (MQ=255)
                                                gCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGAcc                                                                                                 <  4:1366301/151‑1 (MQ=255)
                                                 cAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCt                                                                                                >  4:627856/1‑151 (MQ=255)
                                                 cAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAAAATCGATAGCCTACGCAGTGAcc                                                                                                 >  2:1042506/1‑150 (MQ=255)
                                                  aaCGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTg                                                                                               >  4:490675/1‑151 (MQ=255)
                                                   aCGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTgc                                                                                              >  1:36620/1‑151 (MQ=255)
                                                   aCGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTgc                                                                                              >  2:612454/1‑151 (MQ=255)
                                                    cgccgcGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTgcg                                                                                             <  4:791341/151‑1 (MQ=255)
                                                     gccgcGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTgcgc                                                                                            >  4:384651/1‑151 (MQ=255)
                                                      ccgcGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGcgcg                                                                                           >  2:675054/1‑151 (MQ=255)
                                                       cgcgCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGcgcg                                                                                           >  4:752176/1‑150 (MQ=255)
                                                       cgcgCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACACCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGa                                                                                          <  4:1561064/151‑1 (MQ=255)
                                                          gctggctgGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGaaa                                                                                        >  4:806607/1‑150 (MQ=255)
                                                           ctggctggCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGaaa                                                                                        >  2:1350897/1‑149 (MQ=255)
                                                                  gcAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATca                                                                               <  2:244903/151‑2 (MQ=255)
                                                                    aGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCg                                                                             >  1:119055/1‑151 (MQ=255)
                                                                        aTTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGtata                                                                         <  2:337329/151‑1 (MQ=255)
                                                                              gCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTg                                                                   >  1:224228/1‑151 (MQ=255)
                                                                                tGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGtt                                                                 >  2:1269812/1‑151 (MQ=255)
                                                                                tGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGtt                                                                 >  4:644871/1‑151 (MQ=255)
                                                                                 gAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGtt                                                                 <  2:763329/150‑1 (MQ=255)
                                                                                 gAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGtt                                                                 <  2:938090/150‑1 (MQ=255)
                                                                                     cGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGACAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTa                                                             >  4:337490/1‑150 (MQ=255)
                                                                                         gAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGaa                                                        <  3:339283/151‑1 (MQ=255)
                                                                                         gAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGaa                                                        <  1:36063/151‑1 (MQ=255)
                                                                                         gAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGa                                                         <  4:920907/150‑1 (MQ=255)
                                                                                          aaGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGaaa                                                       <  3:382765/151‑1 (MQ=255)
                                                                                           aGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGaaaa                                                      <  2:667381/151‑1 (MQ=255)
                                                                                           aGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGaaa                                                       >  4:45637/1‑150 (MQ=255)
                                                                                               cGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAGGGGCATCCCGTATACCACTGTTCGTAAAGAAaaca                                                    <  3:144184/149‑1 (MQ=255)
                                                                                               cGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAaacaac                                                  >  4:102345/1‑151 (MQ=255)
                                                                                               cGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAaacaac                                                  <  2:627012/151‑1 (MQ=255)
                                                                                               cGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAaacaac                                                  >  4:1392130/1‑151 (MQ=255)
                                                                                                ggCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACt                                                 <  4:613289/151‑1 (MQ=255)
                                                                                                ggCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACt                                                 >  3:153247/1‑151 (MQ=255)
                                                                                                 gCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTa                                                <  2:1520970/151‑1 (MQ=255)
                                                                                                   cTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACg                                              >  4:724895/1‑151 (MQ=255)
                                                                                                   cTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACg                                              >  4:51180/1‑151 (MQ=255)
                                                                                                    ttGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACgg                                             >  4:157345/1‑151 (MQ=255)
                                                                                                    ttGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACgg                                             >  2:292243/1‑151 (MQ=255)
                                                                                                     tGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACgg                                             <  4:217364/150‑1 (MQ=255)
                                                                                                     tGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACg                                              <  4:888776/149‑1 (MQ=255)
                                                                                                     tGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACg                                              <  4:1106622/149‑1 (MQ=255)
                                                                                                     tGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGc                                            >  2:1474704/1‑151 (MQ=255)
                                                                                                      gACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGcc                                           <  4:1669364/151‑1 (MQ=255)
                                                                                                       aCCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCt                                          <  2:1224405/151‑1 (MQ=255)
                                                                                                           gCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAg                                       <  2:567130/150‑1 (MQ=255)
                                                                                                           gCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGc                                      >  3:36038/1‑151 (MQ=255)
                                                                                                            cGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGc                                      <  2:1669738/150‑1 (MQ=255)
                                                                                                                  ggcgTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACtt                               <  4:713295/151‑1 (MQ=255)
                                                                                                                     gTTCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTcc                            <  4:536067/151‑1 (MQ=255)
                                                                                                                      ttCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCcg                           >  2:491266/1‑151 (MQ=255)
                                                                                                                      ttCACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCcg                           >  4:392797/1‑151 (MQ=255)
                                                                                                                        cACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCcgcg                         >  1:218100/1‑151 (MQ=255)
                                                                                                                        cACTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCcgcg                         >  4:1487823/1‑151 (MQ=255)
                                                                                                                          cTTCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGAt                       >  4:789905/1‑151 (MQ=255)
                                                                                                                           ttCCTGATGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATg                      >  4:1680522/1‑151 (MQ=255)
                                                                                                                                  tGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTCCGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCCCTGTTCGTAAAGAAAACAACTAAGGCCTGAGCATCACTTTCCGCGATGCTAAAGc               >  2:953606/1‑151 (MQ=255)
                                                                                                                                  tGGAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGc               >  4:1507479/1‑151 (MQ=255)
                                                                                                                                    gAAGTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTc             <  4:1616963/151‑1 (MQ=255)
                                                                                                                                       gTGGATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGtg          <  4:932409/151‑1 (MQ=255)
                                                                                                                                         ggATATGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGtgat        <  2:1587698/151‑1 (MQ=255)
                                                                                                                                           atatGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGtgat        >  1:452097/1‑149 (MQ=255)
                                                                                                                                              tGGACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACTCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGTGATGAAGc   >  4:1656400/1‑151 (MQ=255)
                                                                                                                                               ggACGCCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGTGATGAAGcc  >  4:907509/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATGGGCGTTATGGGTGACAAATACGTCGTGGCGCTTAACCTTGCCCCGGCAACGCCGCGCTGGCTGGCAGCTATTCACGCTGAGCCGATGAAGCTCGGCCTTGACCTGCGTGGCGGCGTTCACTTCCTGATGGAAGTGGATATGGACACCGCGCTTGGCAAACTCCAGGAACAAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGTGATGAAGCC  >  NC_000913/427920‑428213

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