Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,142,815 T→C 5.7% F302L (TTC→CTC)  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,8150TC 5.7% 260.5 / 4.0 87F302L (TTC→CTC) PP_3646aldehyde dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base T (38/44);  new base C (2/3);  total (40/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01

GCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGCCGAAGGCGCCAAACTG  >  NC_002947/4142675‑4142964
                                                                                                                                            |                                                                                                                                                     
gcAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAACTCGTCGaag                                                                                                                                              <  1:2267502/150‑3 (MQ=255)
  aaGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAgcgc                                                                                                                                            >  1:1818835/1‑150 (MQ=255)
   aGTCGCCGAACATCATCTTCGCCGACGCCCACCCGGCAAGCACAAACTATCATGCCGTTGCCGGCATCTATGCCGCAGCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGcc                                                                                                                                           <  1:1116489/150‑1 (MQ=255)
     tCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTg                                                                                                                                         >  2:2123825/1‑150 (MQ=255)
      cGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGa                                                                                                                                        <  2:474346/150‑1 (MQ=255)
       gCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGAt                                                                                                                                       <  2:2020014/150‑1 (MQ=255)
        ccGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGAtt                                                                                                                                      <  2:766904/150‑1 (MQ=255)
          gAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGc                                                                                                                                    >  2:1450759/1‑150 (MQ=255)
            aCATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTACTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGccc                                                                                                                                  <  2:1492389/150‑1 (MQ=255)
               tcatcTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCcgcg                                                                                                                               >  2:1607755/1‑150 (MQ=255)
                  tcTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCa                                                                                                                            >  2:103660/1‑150 (MQ=255)
                    ttCGCCGACGCCGGTCTGGAGAGCGCCATCAATGGCGCCGTCGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGTATGAAATCTTCGACTAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAg                                                                                                                          <  2:552748/150‑1 (MQ=255)
                    ttCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCTCGTCCGGGCAGAGCGGCGTGCCCGGACGTCGCCGGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAg                                                                                                                          >  1:168689/1‑150 (MQ=255)
                    ttCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAg                                                                                                                          <  1:1021975/150‑1 (MQ=255)
                     tCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAgc                                                                                                                         <  2:2149518/150‑1 (MQ=255)
                     tCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAgc                                                                                                                         <  1:314121/150‑1 (MQ=255)
                     tCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGCCAGCTGACCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAgc                                                                                                                         <  1:147072/150‑1 (MQ=255)
                         cgacgccgacCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGcgcat                                                                                                                     >  1:417428/1‑150 (MQ=255)
                             gccgacCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGcgcatccgc                                                                                                                 >  1:2028693/1‑150 (MQ=255)
                                gacCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCt                                                                                                                                          <  2:2224603/122‑1 (MQ=255)
                                gacCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCt                                                                                                                                          >  1:2224603/1‑122 (MQ=255)
                                gacCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGcgcatccgcatc                                                                                                              <  1:2357813/150‑1 (MQ=255)
                                gacCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAACTCGTCGaagat                                                                                                                                            >  2:1913170/1‑116 (MQ=255)
                                gacCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAACTCGTCGaagat                                                                                                                                            <  1:1913170/120‑5 (MQ=255)
                                 acCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCg                                                                                                             >  2:345700/1‑150 (MQ=255)
                                  ccTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCgg                                                                                                            >  2:2335939/1‑150 (MQ=255)
                                   cTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGcgcat                                                                                                                     <  2:653581/140‑1 (MQ=255)
                                   cTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGcgcat                                                                                                                     >  1:653581/1‑140 (MQ=255)
                                      gACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAAc                                                                                                        <  2:1171570/150‑1 (MQ=255)
                                           cgcCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCa                                                                                                   <  1:953395/150‑1 (MQ=255)
                                            gcCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAg                                                                                                  >  1:130224/1‑150 (MQ=255)
                                            gcCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAg                                                                                                  <  1:1127999/150‑1 (MQ=255)
                                            gcCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAg                                                                                                  >  2:842489/1‑150 (MQ=255)
                                            gcCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAg                                                                                                  >  2:2341121/1‑150 (MQ=255)
                                              cATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTCGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATGCGCCACCCGCCGga                                                                                                >  2:170337/1‑150 (MQ=255)
                                               aTCACTGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGAGCATCCGCATCGGCAACCCGCAGgac                                                                                               <  1:824516/150‑1 (MQ=255)
                                               aTCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGgac                                                                                               <  2:1777332/150‑1 (MQ=255)
                                                 cAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGgacga                                                                                             <  1:2041964/150‑1 (MQ=255)
                                                  aaTGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCa                                                                                                   >  1:2307029/1‑143 (MQ=255)
                                                  aaTGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCa                                                                                                   <  2:2307029/143‑1 (MQ=255)
                                                  aaTGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGgacgac                                                                                            <  2:104386/150‑1 (MQ=255)
                                                   aTGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGacgacg                                                                                           >  2:175784/1‑150 (MQ=255)
                                                           gTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGaa                                                                                   >  2:105961/1‑150 (MQ=255)
                                                             ggCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAAt                                                                                 >  2:2055853/1‑150 (MQ=255)
                                                              gCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATg                                                                                <  2:1818835/150‑1 (MQ=255)
                                                               ccGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATgg                                                                               <  2:1820029/150‑1 (MQ=255)
                                                               ccGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATgg                                                                               >  1:191799/1‑150 (MQ=255)
                                                                cGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTAGAGCGCCTGATTGCCCgcgc                                                                                                                              >  2:498439/1‑102 (MQ=255)
                                                                cGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTAGAGCGCCTGATTGCCCgcgc                                                                                                                              <  1:498439/102‑1 (MQ=255)
                                                                cGGCATCTATGCCACGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATggg                                                                              >  2:690264/1‑150 (MQ=255)
                                                                 ggCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGAAGCCAGCGAAATGGGc                                                                             <  1:2299695/150‑1 (MQ=255)
                                                                  gCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGcc                                                                            <  2:468852/150‑1 (MQ=255)
                                                                    aTCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGcgcatc                                                                                                                    >  2:212623/1‑108 (MQ=255)
                                                                    aTCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGcgcatc                                                                                                                    <  1:212623/108‑1 (MQ=255)
                                                                     tCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGa                                                                         <  1:345700/150‑1 (MQ=255)
                                                                       tATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATg                                                                       <  1:1450759/150‑1 (MQ=255)
                                                                           ccGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATGGGCAACCCGCAGGACGACGCCAGCGAAATGGGc                                                                             <  2:2112090/140‑1 (MQ=255)
                                                                           ccGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATGGGCAACCCGCAGGACGACGCCAGCGAAATGGGc                                                                             >  1:2112090/1‑140 (MQ=255)
                                                                                 ccGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCgcgc                                                             <  2:2028693/150‑1 (MQ=255)
                                                                                    ggCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCACGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAACTGGGCCCGATGGCCACCGCGCAAc                                                          >  2:311818/1‑150 (MQ=255)
                                                                                      cAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGgacga                                                                                             <  2:1825081/113‑1 (MQ=255)
                                                                                      cAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGgacga                                                                                             >  1:1825081/1‑113 (MQ=255)
                                                                                         aGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTg                                                     >  1:1168357/1‑150 (MQ=255)
                                                                                          gCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGtt                                                                                                                                                      >  2:256352/1‑52 (MQ=255)
                                                                                          gCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGtt                                                                                                                                                      <  1:256352/52‑1 (MQ=255)
                                                                                            tGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGcgcatc                                                                                                                    <  2:735789/84‑1 (MQ=255)
                                                                                            tGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGcgcatc                                                                                                                    >  1:735789/1‑84 (MQ=255)
                                                                                            tGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCa                                                  <  1:243121/150‑1 (MQ=255)
                                                                                                     ggCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGaa                                         >  2:463427/1‑150 (MQ=255)
                                                                                                        tCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGc                                      >  2:37754/1‑150 (MQ=255)
                                                                                                            gcCGGCGGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCtgg                                  <  2:853139/150‑1 (MQ=255)
                                                                                                                   ggTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATGGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAAGGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTggcggc                           <  1:842489/150‑1 (MQ=255)
                                                                                                                     tGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGc                                                                             >  2:1159897/1‑98 (MQ=255)
                                                                                                                     tGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGc                                                                             <  1:1159897/98‑1 (MQ=255)
                                                                                                                         ggATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAAt                                                                                 <  1:1150297/90‑1 (MQ=255)
                                                                                                                         ggATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAAt                                                                                 >  2:1150297/1‑90 (MQ=255)
                                                                                                                            tGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGc                  >  2:1237570/1‑150 (MQ=255)
                                                                                                                             gAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGcc                 >  1:759005/1‑150 (MQ=255)
                                                                                                                               aaTCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGgacgac                                                                                            <  2:543621/73‑1 (MQ=255)
                                                                                                                               aaTCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCAGCCGCATCGGCAACCCGCAGgacgac                                                                                            >  1:543621/1‑73 (MQ=255)
                                                                                                                                 tCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGCCGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGCCGAAg             <  1:1973439/150‑1 (MQ=255)
                                                                                                                                 tCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGCCGAAg             <  1:817046/150‑1 (MQ=255)
                                                                                                                                  gctCGACGAACTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAAt                                                                                 >  1:1275141/3‑81 (MQ=255)
                                                                                                                                  gctCGACGAACTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAAt                                                                                 <  2:1275141/79‑1 (MQ=255)
                                                                                                                                   ttCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGCAGCAAGGCTACCGGGCCCGACGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGCCGAAGgc           <  1:2123825/150‑1 (MQ=255)
                                                                                                                                      gacgaGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCGCCTACATGGGCACGCAATCCACTGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGCCGAAGGCGcc        >  1:542798/1‑150 (MQ=255)
                                                                                                                                       acgaGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGCCGAAGgccttg       <  2:10554/150‑5 (MQ=255)
                                                                                                                                       acgaGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGCCGAAGgccttg       <  2:13263/150‑5 (MQ=255)
                                                                                                                                           gTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGCCGAAGGCGCCAAACt   >  2:1090190/1‑150 (MQ=255)
                                                                                                                                            ttCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGCCGAAGGCGCCAAACTg  >  1:929036/1‑150 (MQ=255)
                                                                                                                                            ttCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGCCGAAGGCGCCAAACTg  >  1:1167558/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                     
GCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGCCGAAGGCGCCAAACTG  >  NC_002947/4142675‑4142964

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