Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,142,724 C→T 7.0% I271I (ATC→ATT PP_3646 → aldehyde dehydrogenase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,142,7240CT 7.0% 404.8 / 12.7 138I271I (ATC→ATTPP_3646aldehyde dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base C (60/68);  new base T (5/5);  total (65/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 138 total aligned reads displayed.

GGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACG  >  NC_002947/4142577‑4142872
                                                                                                                                                   |                                                                                                                                                    
ggTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCaa                                                                                                                                                    <  1:2154066/150‑1 (MQ=255)
 gggcgcAAGATCGCCTTCACCGGCGGTGCAGCCACGGCCCGCCACGTGGTACGCATCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAAt                                                                                                                                                   <  2:727539/148‑1 (MQ=255)
 gTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAAt                                                                                                                                                   >  1:1931150/1‑150 (MQ=255)
     gcAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGgcg                                                                                                                                               >  2:813919/1‑150 (MQ=255)
     gcAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGgcg                                                                                                                                               >  2:202433/1‑150 (MQ=255)
     gcAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATAATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGgcg                                                                                                                                               >  2:2864516/1‑150 (MQ=255)
         gATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGt                                                                                                                                           >  1:1461970/1‑150 (MQ=255)
         gATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGt                                                                                                                                           <  2:2864835/150‑1 (MQ=255)
         gATCGCCGTCCCCGGCGGTGCCGCCCCCGCCCGACACGTGGTACGCCTCAGCGCCGCAAACTGCGCCAAGCTGTCACTGGAACTGGGCGGCAAGTCGCTGAACATCATCTTCGCCGACGCCGACATGGACAGCGCCATCAATGGCGccct                                                                                                                                           <  2:542442/150‑3 (MQ=255)
          aTCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTCCGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACGGCGCCATTGATGGCGCTGtt                                                                                                                                          <  2:2102506/150‑2 (MQ=255)
           tCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTAGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTgg                                                                                                                                         <  1:290088/150‑1 (MQ=255)
              ccTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCg                                                                                                                                      >  2:1538235/1‑150 (MQ=255)
              ccTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCg                                                                                                                                      >  2:438245/1‑150 (MQ=255)
              ccTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCg                                                                                                                                      <  2:183850/150‑1 (MQ=255)
               cTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCgg                                                                                                                                     >  1:970893/1‑150 (MQ=255)
                ttCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATTGATGgcgctgtcc                                                                                                                                         >  1:2663022/1‑140 (MQ=255)
                ttCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATTGATGgcgctatcc                                                                                                                                         <  2:2663022/145‑6 (MQ=255)
                ttCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGc                                                                                                                                    >  2:835490/1‑150 (MQ=255)
                     cggcggTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTa                                                                                                                               >  1:2378071/1‑150 (MQ=255)
                      ggcggTGGCGGCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTAt                                                                                                                              >  1:942366/1‑150 (MQ=255)
                      ggcggTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGTGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTAt                                                                                                                              <  1:29522/150‑1 (MQ=255)
                             ccgccACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGt                                                                                                                       >  1:2958610/1‑150 (MQ=255)
                             ccgccACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGt                                                                                                                       >  2:2678350/1‑150 (MQ=255)
                               gccACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTct                                                                                                                     >  2:394061/1‑149 (MQ=255)
                                 cACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCgg                                                                                                                   <  1:698210/150‑1 (MQ=255)
                                  aCCGCCCGCCACGTGGTACGCAGCAGCCCCGAAAACTTCGCTACACTGGCGCTGGAACTGGGCGGCAAGTCGCCGCCTATCATCCTCGCTCCTGTTCGTTGGTACGACGCCGGCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCggg                                                                                                                  <  1:835490/150‑1 (MQ=255)
                                   ccgcccgcCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGc                                                                                                                 >  2:2127431/1‑150 (MQ=255)
                                       ccgccACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAAt                                                                                                                                                   <  2:2515539/112‑1 (MQ=255)
                                       ccgccACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAAt                                                                                                                                                   >  1:2515539/1‑112 (MQ=255)
                                       ccgccACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAg                                                                                                             >  1:844021/1‑150 (MQ=255)
                                       ccgccACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAg                                                                                                             >  2:2318086/1‑150 (MQ=255)
                                           cACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGc                                                                                                         <  2:1931150/150‑1 (MQ=255)
                                           cACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGc                                                                                                         <  2:1165969/150‑1 (MQ=255)
                                           cACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGc                                                                                                         <  1:783633/150‑1 (MQ=255)
                                              gTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATg                                                                                                                                                  >  1:2777391/1‑106 (MQ=255)
                                              gTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTg                                                                                                      <  1:1320684/150‑1 (MQ=255)
                                               tGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTAGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTgg                                                                                                     >  1:591268/1‑150 (MQ=255)
                                                ggTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGc                                                                                                    <  2:1591877/150‑1 (MQ=255)
                                                  tACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCg                                                                                                  <  1:1949674/150‑1 (MQ=255)
                                                        gcagcGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTcgc                                                                                            <  1:616424/150‑1 (MQ=255)
                                                           gcgcCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGcc                                                                                         >  2:450490/1‑150 (MQ=255)
                                                            cgcCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCct                                                                                        <  2:942366/150‑1 (MQ=255)
                                                              ccGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCctgc                                                                                      <  2:796927/150‑1 (MQ=255)
                                                               cGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGGGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCctgct                                                                                     >  1:2847104/1‑150 (MQ=255)
                                                                gAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCctgctg                                                                                    >  2:2588161/1‑150 (MQ=255)
                                                                gAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCctgctg                                                                                    >  1:2557522/1‑150 (MQ=255)
                                                                gAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCctgctg                                                                                    <  1:3014786/150‑1 (MQ=255)
                                                                gAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGAGCCATCAATGGCGTCGGTTCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCctgctg                                                                                    <  1:2678350/150‑1 (MQ=255)
                                                                  aaaCTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGt                                                                                  <  2:2844287/150‑1 (MQ=255)
                                                                    aCTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGc                                                                                <  2:3004731/150‑1 (MQ=255)
                                                                     cTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCa                                                                               >  2:1480630/1‑150 (MQ=255)
                                                                      ttCGCCAAGCTGTCGCTGGAGCTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAg                                                                              >  2:1855089/1‑150 (MQ=255)
                                                                      ttCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCGTGCCGGCGCAg                                                                              >  2:838/1‑150 (MQ=255)
                                                                      ttCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCAGGGCAAAGCTGCGTGGCCGGGTCGCCGCGGCTTGTGCAg                                                                              <  2:1978151/150‑1 (MQ=255)
                                                                        cGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGa                                                                            <  1:1538235/150‑1 (MQ=255)
                                                                          ccAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATg                                                                          <  2:1132380/150‑1 (MQ=255)
                                                                           cAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTg                                                                                                      <  1:93888/121‑1 (MQ=255)
                                                                           cAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTg                                                                                                      >  2:93888/1‑121 (MQ=255)
                                                                           cAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTg                                                                                 >  1:1958965/1‑142 (MQ=255)
                                                                            aaGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGaa                                                                        >  2:2000205/1‑150 (MQ=255)
                                                                             aGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGaaa                                                                       <  1:153733/150‑1 (MQ=255)
                                                                              gCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATg                                                                          >  1:1070698/1‑146 (MQ=255)
                                                                                  tCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTc                                                                  <  1:260119/150‑1 (MQ=255)
                                                                                    gCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGCGAAAATCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTcga                                                                <  2:2733002/150‑1 (MQ=255)
                                                                                         aaCTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGt                                                           <  2:2557522/150‑1 (MQ=255)
                                                                                              ggcggcAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTtcgtc                                                      <  2:1888912/150‑1 (MQ=255)
                                                                                              ggcggcAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTtcgtc                                                      <  2:844021/150‑1 (MQ=255)
                                                                                              ggcggcAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTtcgtc                                                      <  1:438245/150‑1 (MQ=255)
                                                                                                 ggcAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAACTCGTCGaa                                                   <  1:813919/150‑2 (MQ=255)
                                                                                                    aaGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAgcgc                                                <  2:1549858/150‑1 (MQ=255)
                                                                                                     aGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGcc                                               <  2:2847104/150‑1 (MQ=255)
                                                                                                      gTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCt                                              <  1:2270422/150‑1 (MQ=255)
                                                                                                              aCATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGccc                                      <  2:1890183/150‑1 (MQ=255)
                                                                                                                atcatcTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCcgc                                    >  1:1610131/1‑150 (MQ=255)
                                                                                                                 tcatcTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCcgcg                                   >  2:122347/1‑150 (MQ=255)
                                                                                                                    tcTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCa                                >  2:2445381/1‑150 (MQ=255)
                                                                                                                      ttCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTtcgtc                                                      <  2:147433/126‑1 (MQ=255)
                                                                                                                      ttCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAg                              >  2:2265990/1‑150 (MQ=255)
                                                                                                                             acgccgacCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGcgcatcc                       >  1:2711268/1‑150 (MQ=255)
                                                                                                                             acgccgacCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGcgcatcc                       <  1:1480630/150‑1 (MQ=255)
                                                                                                                               gccgacCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGcgcatccgc                     >  2:22416/1‑150 (MQ=255)
                                                                                                                                 cgacCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTtcgt                                                       <  2:498836/114‑1 (MQ=255)
                                                                                                                                 cgacCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGcgcatccgcat                   <  1:2588161/150‑1 (MQ=255)
                                                                                                                                 cgacCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGcgcatccgcat                   >  2:2413307/1‑150 (MQ=255)
                                                                                                                                 cgacCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGcgcatccgcat                   <  1:1623712/150‑1 (MQ=255)
                                                                                                                                  gacCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGcgcatccgcatc                  <  1:2499745/150‑1 (MQ=255)
                                                                                                                                     cTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGc               >  2:644526/1‑150 (MQ=255)
                                                                                                                                       ggACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCATGATGAAATCTTCGACGGGTTCGTCTGCGGCCTGTTTGGGCGCTCGAAGGGCATCCGCATCGGGaa             >  2:1412637/1‑150 (MQ=255)
                                                                                                                                        gACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACAAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAAc            <  1:2364319/150‑1 (MQ=255)
                                                                                                                                         acggcgcCATTGATGGCGCTGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGaa                                                                        <  2:2955911/86‑1 (MQ=255)
                                                                                                                                         acggcgcCATTGATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGAtt                                          >  2:2802251/4‑119 (MQ=255)
                                                                                                                                         acggcgcCATTGATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGAtt                                          <  1:2802251/116‑1 (MQ=255)
                                                                                                                                         acggcgcCATTGATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGcc                                       <  1:2986717/119‑1 (MQ=255)
                                                                                                                                          cAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAAccc          <  1:352845/150‑1 (MQ=255)
                                                                                                                                          cAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAAccc          <  2:2378071/150‑1 (MQ=255)
                                                                                                                                              gcCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAg      <  1:2096985/150‑1 (MQ=255)
                                                                                                                                              gcCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAg      <  2:1461970/150‑1 (MQ=255)
                                                                                                                                               ccATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCTGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAgg     <  1:394061/150‑1 (MQ=255)
                                                                                                                                                 aTCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGgac   >  2:2654166/1‑150 (MQ=255)
                                                                                                                                                  tCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAAGCGTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGgacg  >  2:1493084/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGTGCGCAAGATCGCCTTCACCGGCGGTGCCGCCACCGCCCGCCACGTGGTACGCAGCAGCGCCGAAAACTTCGCCAAGCTGTCGCTGGAACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACG  >  NC_002947/4142577‑4142872

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