Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 165,917 C→T A396A (GCC→GCT mrcB → fused glycosyl transferase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913165,9170CT98.8% 342.7 / ‑4.4 86A396A (GCC→GCTmrcBfused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (38/47);  total (38/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

AGGCGCTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCCCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGTAAAAGATCTCTCCGGCGTGAAGATCTTCACTACCTTTGACTCG  >  NC_000913/165769‑166061
                                                                                                                                                    |                                                                                                                                                
agGCGCTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCTAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCg                                                                                                                                                <  1:90073/151‑1 (MQ=255)
agGCGCTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCg                                                                                                                                                <  1:694680/151‑1 (MQ=255)
   cgcTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTcc                                                                                                                                             <  2:469074/151‑1 (MQ=255)
   cgcTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTcc                                                                                                                                             <  1:978664/151‑1 (MQ=255)
    gcTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCg                                                                                                                                            <  2:552630/151‑1 (MQ=255)
            tCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGt                                                                                                                                    >  2:240324/1‑151 (MQ=255)
            tCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGt                                                                                                                                    >  2:750798/1‑151 (MQ=255)
             cGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGtt                                                                                                                                   <  2:947991/151‑1 (MQ=255)
               gTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCa                                                                                                                                 <  2:1462069/151‑1 (MQ=255)
               gTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCa                                                                                                                                 <  2:374806/151‑1 (MQ=255)
                 aTGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAg                                                                                                                                <  2:44863/150‑1 (MQ=255)
                 aTGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGc                                                                                                                               <  2:100973/151‑1 (MQ=255)
                 aTGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGc                                                                                                                               >  2:36224/1‑151 (MQ=255)
                   ggTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGcc                                                                                                                              <  2:755094/150‑1 (MQ=255)
                     tGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCcgcg                                                                                                                           >  1:532358/1‑151 (MQ=255)
                        aaGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCgg                                                                                                                        >  1:731850/1‑151 (MQ=255)
                        aaGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCgg                                                                                                                        >  2:172757/1‑151 (MQ=255)
                         aGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCgg                                                                                                                        >  2:788037/1‑150 (MQ=255)
                         aGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGt                                                                                                                       >  1:786229/1‑151 (MQ=255)
                          ggggCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTg                                                                                                                      <  2:1030389/151‑1 (MQ=255)
                           gggCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTgg                                                                                                                     <  1:1094566/151‑1 (MQ=255)
                              cGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGt                                                                                                                  <  1:1714443/151‑1 (MQ=255)
                               gTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTg                                                                                                                 >  2:1353006/1‑151 (MQ=255)
                                tCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTg                                                                                                                 >  1:1159102/1‑150 (MQ=255)
                                  cATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGAtc                                                                                                              <  1:478641/151‑1 (MQ=255)
                                      tACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGAtctctc                                                                                                          >  1:686548/1‑151 (MQ=255)
                                       aCAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGAtctctc                                                                                                          <  1:1507092/150‑1 (MQ=255)
                                         aaCCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTctcctc                                                                                                       >  1:106149/1‑151 (MQ=255)
                                           cccGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTcag                                                                                                     <  2:462372/151‑1 (MQ=255)
                                              gTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTcagcca                                                                                                  <  2:1003528/151‑1 (MQ=255)
                                                  cGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCt                                                                                              <  2:7554/151‑1 (MQ=255)
                                                     aaCCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCttt                                                                                            >  2:818448/1‑150 (MQ=255)
                                                     aaCCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCttt                                                                                            <  1:1090864/150‑1 (MQ=255)
                                                            aaCTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACt                                                                                    >  2:1474353/1‑151 (MQ=255)
                                                                ggCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGt                                                                                 <  1:1093395/150‑1 (MQ=255)
                                                                ggCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCCCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTgg                                                                                  <  2:990261/149‑1 (MQ=255)
                                                                 gcgcTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGc                                                                               <  1:1137125/151‑1 (MQ=255)
                                                                      ggAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTc                                                                            >  1:1612662/1‑149 (MQ=255)
                                                                      ggAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCa                                                                           <  2:271075/150‑1 (MQ=255)
                                                                       gAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGGGCGTCAgg                                                                         >  1:1528993/1‑151 (MQ=255)
                                                                        aGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGa                                                                        <  1:1056533/151‑1 (MQ=255)
                                                                          cGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGc                                                                      >  1:1016408/1‑151 (MQ=255)
                                                                           gACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCt                                                                     >  2:778831/1‑151 (MQ=255)
                                                                           gACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCt                                                                     <  2:1060370/151‑1 (MQ=255)
                                                                           gACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCt                                                                     >  1:1071721/1‑151 (MQ=255)
                                                                           gACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCt                                                                     >  1:1372879/1‑151 (MQ=255)
                                                                            aCGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTg                                                                    >  1:1098896/1‑151 (MQ=255)
                                                                              gTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGc                                                                   >  1:650583/1‑150 (MQ=255)
                                                                              gTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCa                                                                  <  2:432546/151‑1 (MQ=255)
                                                                                  tCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGc                                                               >  1:394233/1‑150 (MQ=255)
                                                                                     ggTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCaaaa                                                           <  2:1345886/151‑1 (MQ=255)
                                                                                         ctgcGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTgg                                                       >  2:1248481/1‑151 (MQ=255)
                                                                                           gcGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTggg                                                      <  2:279777/150‑1 (MQ=255)
                                                                                           gcGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGc                                                     <  2:685563/151‑1 (MQ=255)
                                                                                           gcGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGc                                                     <  2:500586/151‑1 (MQ=255)
                                                                                           gcGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGc                                                     <  1:1308040/151‑1 (MQ=255)
                                                                                             gTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCg                                                    <  1:709034/150‑1 (MQ=255)
                                                                                                 gctgcAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATaa                                                <  2:475282/150‑1 (MQ=255)
                                                                                                 gctgcAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATa                                                 <  2:862807/149‑1 (MQ=255)
                                                                                                 gctgcAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAg                                               <  2:167042/151‑1 (MQ=255)
                                                                                                 gctgcAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAg                                               <  2:1301761/151‑1 (MQ=255)
                                                                                                    gcaacagcaacaGATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGt                                             <  1:306887/150‑1 (MQ=255)
                                                                                                            aacagATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGTAAAAGAtct                                    >  2:16757/1‑151 (MQ=255)
                                                                                                             acagATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGTAAAAGAtctc                                   <  1:1561948/151‑1 (MQ=255)
                                                                                                             acagATTATTGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGTAAAAGAtctc                                   <  1:1103885/151‑1 (MQ=255)
                                                                                                                    attGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGTAAAAGATCTCTCCGGCg                            >  2:362512/1‑151 (MQ=255)
                                                                                                                      tGATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGTAAAAGATCTCTCCGGCGt                           <  1:432659/150‑1 (MQ=255)
                                                                                                                       gATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGTAAAAGATCTCTCCGGCGTg                          <  1:1232834/150‑1 (MQ=255)
                                                                                                                       gATCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGTAAAAGATCTCTCCGGCGTGa                         <  2:1313732/151‑1 (MQ=255)
                                                                                                                        aTCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGTAAAAGATCTCTCCGGCGTGaa                        <  2:154012/151‑1 (MQ=255)
                                                                                                                        aTCAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGTAAAAGATCTCTCCGGCGTGa                         >  2:331615/1‑150 (MQ=255)
                                                                                                                          cAAGAACTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGTAAAAGATCTCTCCGGCGTGa                         >  2:549147/1‑148 (MQ=255)
                                                                                                                              aaCTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGTAAAAGATCTCTCCGGCGTGAAGATCtt                  >  2:174898/1‑151 (MQ=255)
                                                                                                                              aaCTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGTAAAAGATCTCTCCGGCGTGAAGATCtt                  <  1:818383/151‑1 (MQ=255)
                                                                                                                               aCTCTATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGTAAAAGATCTCTCCGGCGTGAAGATCTTc                 >  2:535471/1‑151 (MQ=255)
                                                                                                                                ctctATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGTAAAAGATCTCTCCGGCGTGAAGATCtt                  >  1:547558/1‑149 (MQ=255)
                                                                                                                                ctctATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGTAAAAGATCTCTCCGGCGTGAAGATCTTc                 >  1:1160095/1‑150 (MQ=255)
                                                                                                                                  ctATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGTAAAAGATCTCTCCGGCGTGAAGATCTTCACt              >  1:1055463/1‑151 (MQ=255)
                                                                                                                                  ctATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGTAAAAGATCTCTCCGGCGTGAAGATCTTCACt              >  1:1618957/1‑151 (MQ=255)
                                                                                                                                   tATGACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGTAAAAGATCTCTCCGGCGTGAAGATCTTCACt              >  2:619839/1‑150 (MQ=255)
                                                                                                                                      gACATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGTAAAAGATCTCTCCGGCGTGAAGATCTTCACTAcc           <  2:1305362/150‑1 (MQ=255)
                                                                                                                                        cATGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGTAAAAGATCTCTCCGGCGTGAAGATCTTCACTACCttt        >  1:1327310/1‑151 (MQ=255)
                                                                                                                                         aTGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGTAAAAGATCTCTCCGGCGTGAAGATCTTCACTACCTTTg       >  1:120017/1‑151 (MQ=255)
                                                                                                                                         aTGTTGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGTAAAAGATCTCTCCGGCGTGAAGATCTTCACTACCTTTg       >  2:1078784/1‑151 (MQ=255)
                                                                                                                                             tGAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGTAAAAGATCTCTCCGGCGTGAAGATCTTCACTACCTTTGACt    <  1:685152/150‑1 (MQ=255)
                                                                                                                                              gAGTGCTCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGTAAAAGATCTCTCCGGCGTGAAGATCTTCACTACCTTTGACTCg  <  2:1612760/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
AGGCGCTGTTAGTCGGTATGGTGAAAGGGGCGTCCATCTACAACCCGTGGCGTAACCCAAAACTGGCGCTGGAGCGACGTAATCTGGTGCTGCGTCTGCTGCAACAGCAACAGATTATTGATCAAGAACTCTATGACATGTTGAGTGCCCGTCCGCTGGGGGTTCAGCCGCGCGGTGGGGTGATCTCTCCTCAGCCAGCCTTTATGCAACTGGTGCGTCAGGAGCTGCAGGCAAAACTGGGCGATAAGGTAAAAGATCTCTCCGGCGTGAAGATCTTCACTACCTTTGACTCG  >  NC_000913/165769‑166061

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