Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,003,786 G→A 7.1% R182Q (CGG→CAG)  PP_2625 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,003,7860GA 7.1% 553.2 / 22.1 182R182Q (CGG→CAG) PP_2625hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (84/85);  new base A (5/8);  total (89/93)
Fisher's exact test for biased strand distribution p-value = 5.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 191 total aligned reads displayed.

TGCATACAGGCCATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACATGTCGGTG  >  NC_002947/3003638‑3003933
                                                                                                                                                    |                                                                                                                                                   
 gCATACAGGCCATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCcgggcggg                                                                                                                                                   >  2:2706741/1‑150 (MQ=255)
     aCAGGCCATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCt                                                                                                                                               >  1:2318123/1‑150 (MQ=255)
     aCAGGCCATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCt                                                                                                                                               <  1:4303322/150‑1 (MQ=255)
      cAGGCCATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGTGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTg                                                                                                                                              <  2:842095/150‑1 (MQ=255)
       aGGCCATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGc                                                                                                                                             >  2:1772024/1‑150 (MQ=255)
        ggCCATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGcc                                                                                                                                            >  2:3014367/1‑150 (MQ=255)
         gCCATCTCCCGGCTACAGCTGAGTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGccc                                                                                                                                           <  2:244782/150‑1 (MQ=255)
           cATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGc                                                                                                                                         >  1:4427622/1‑150 (MQ=255)
           cATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGc                                                                                                                                         >  2:2709675/1‑150 (MQ=255)
            aTCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGcc                                                                                                                                        >  1:3760223/1‑150 (MQ=255)
              ctcCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACTCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCAt                                                                                                                                      <  1:3251469/150‑1 (MQ=255)
               tcCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATc                                                                                                                                     <  2:4553024/150‑1 (MQ=255)
               tcCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATc                                                                                                                                     >  1:591018/1‑150 (MQ=255)
               tcCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATc                                                                                                                                     >  1:5198134/1‑150 (MQ=255)
                 ccGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGGCGGGCAGGCCCCGCCCGGCATCGa                                                                                                                                   <  1:1835191/150‑1 (MQ=255)
                 ccGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGGCGGGCAGGCCCCGCCCGGCATCGa                                                                                                                                   <  2:3704194/150‑1 (MQ=255)
                 ccGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGGCGGGCAGGCCCCGCCCGGCATCGa                                                                                                                                   <  2:1598205/150‑1 (MQ=255)
                 ccGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGa                                                                                                                                   >  1:3461616/1‑150 (MQ=255)
                 ccGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGa                                                                                                                                   <  2:914499/150‑1 (MQ=255)
                 ccGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGa                                                                                                                                   <  1:3108205/150‑1 (MQ=255)
                 ccGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGa                                                                                                                                   <  1:4860291/150‑1 (MQ=255)
                 ccGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGa                                                                                                                                   <  2:1855500/150‑1 (MQ=255)
                      tACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGa                                                                                                                              <  2:780976/150‑1 (MQ=255)
                         aGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCCAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCCGCGAGCAGATCa                                                                                                                           >  1:4113480/1‑150 (MQ=255)
                         aGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCa                                                                                                                           <  1:2048974/150‑1 (MQ=255)
                            tGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGa                                                                                                                        >  1:1112426/1‑150 (MQ=255)
                            tGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGAGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGa                                                                                                                        <  2:3865806/150‑1 (MQ=255)
                             gATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGaa                                                                                                                       <  1:2706741/150‑1 (MQ=255)
                             gATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGaa                                                                                                                       <  1:2343875/150‑1 (MQ=255)
                             gATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGaa                                                                                                                       >  1:2150826/1‑150 (MQ=255)
                              aTTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCcgggcgg                                                                                                                                                    >  1:82136/1‑120 (MQ=255)
                              aTTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCcgggcgg                                                                                                                                                    <  2:82136/120‑1 (MQ=255)
                              aTTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAg                                                                                                                      >  2:4849013/1‑150 (MQ=255)
                                  ccAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGa                                                                                                                  >  2:3403682/1‑150 (MQ=255)
                                  ccAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGa                                                                                                                  >  2:3403886/1‑150 (MQ=255)
                                   cAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGAt                                                                                                                 >  1:4361177/1‑150 (MQ=255)
                                   cAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGAt                                                                                                                 >  1:4502020/1‑150 (MQ=255)
                                   cAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGAt                                                                                                                 >  1:2126646/1‑150 (MQ=255)
                                    aGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGTCCGCCATCGAGCAGATCAGGAAGCTGATc                                                                                                                <  1:3278784/150‑1 (MQ=255)
                                         tCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGc                                                                                                           >  1:5310272/1‑150 (MQ=255)
                                            gATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCa                                                                                                        >  2:4633998/1‑150 (MQ=255)
                                              tCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGc                                                                                                      <  1:4374954/150‑1 (MQ=255)
                                              tCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGc                                                                                                      <  2:4427622/150‑1 (MQ=255)
                                                 aGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGc                                                                                                                                                 <  2:1334537/104‑1 (MQ=255)
                                                         tGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGtt                                                                                           <  1:5141981/150‑1 (MQ=255)
                                                         tGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGtt                                                                                           <  1:338235/150‑1 (MQ=255)
                                                           ggTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTtgt                                                                                         <  2:2936733/150‑1 (MQ=255)
                                                            gTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTtgtg                                                                                        >  2:516960/1‑150 (MQ=255)
                                                             tCATGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTgtga                                                                                       <  2:3142575/150‑1 (MQ=255)
                                                             tCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTgtga                                                                                       <  2:3604343/150‑1 (MQ=255)
                                                                 ggTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTgtgagtga                                                                                   <  2:2128300/150‑1 (MQ=255)
                                                                   tACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATc                                                                                 <  2:387710/150‑1 (MQ=255)
                                                                      gCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACg                                                                              >  2:4726572/1‑150 (MQ=255)
                                                                       ccTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACgg                                                                             <  1:644771/150‑1 (MQ=255)
                                                                        cTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGc                                                                            >  1:1082928/1‑150 (MQ=255)
                                                                           aGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAg                                                                         <  2:2474132/150‑1 (MQ=255)
                                                                               gggAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGc                                                                                                      >  2:4144398/1‑117 (MQ=255)
                                                                                  tacgCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTa                                                                                                   <  1:4359083/114‑1 (MQ=255)
                                                                                  tacgCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTa                                                                                                   >  2:4359083/4‑117 (MQ=255)
                                                                                    ggCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAg                                                                         >  1:3920435/1‑141 (MQ=255)
                                                                                     gCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTgtgagt                                                                                     >  1:2904356/1‑128 (MQ=255)
                                                                                         gTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGAc                                                           >  2:4700658/1‑150 (MQ=255)
                                                                                         gTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGAc                                                           >  2:4700676/1‑150 (MQ=255)
                                                                                            tACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCa                                                                                                        >  1:5421688/1‑102 (MQ=255)
                                                                                             aCCCGTGAGACGCCCTGCCAAGTGATGGAGACCGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTc                                                       <  2:2612025/150‑1 (MQ=255)
                                                                                               ccGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGa                                                     <  1:2788745/150‑1 (MQ=255)
                                                                                               ccGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGa                                                     >  2:2991816/1‑150 (MQ=255)
                                                                                                cGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGAt                                                    <  2:1082928/150‑1 (MQ=255)
                                                                                                 gTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATg                                                   <  1:678402/150‑1 (MQ=255)
                                                                                                  tGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATgc                                                  >  2:3218566/1‑150 (MQ=255)
                                                                                                        gCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAgg                                            <  2:2318123/150‑1 (MQ=255)
                                                                                                        gCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAgg                                            <  2:147443/150‑1 (MQ=255)
                                                                                                          ccTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGcc                                          >  1:2164498/1‑150 (MQ=255)
                                                                                                            tGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTg                                        >  2:118910/1‑150 (MQ=255)
                                                                                                              gCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGc                                      <  2:5198134/150‑1 (MQ=255)
                                                                                                                aaGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGccc                                    >  2:3278263/1‑150 (MQ=255)
                                                                                                                   tGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCt                                 >  1:4688179/1‑150 (MQ=255)
                                                                                                                    gATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCtt                                <  1:4582609/150‑1 (MQ=255)
                                                                                                                       ggAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGAt                                                                                  <  1:2712200/97‑1 (MQ=255)
                                                                                                                       ggAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGAt                                                                                  >  2:2712200/1‑97 (MQ=255)
                                                                                                                       ggAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGAt                             >  1:1303003/1‑150 (MQ=255)
                                                                                                                         agaACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCa                           >  1:3833184/1‑150 (MQ=255)
                                                                                                                            aCGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGcc                        <  1:5003848/150‑1 (MQ=255)
                                                                                                                             cGGCCAAGTGCTCGATGCCGGGCGGGGCCGGCCCGCCAGCGAGCAGATCAGGAAGCTGATCGGTGCCCCGCGGCTCGAGGGTGGGAGTGATCACGGCGCGGTCGACGGGCTGACCTGCGATGCGCCAGGCCTGGCCCCCGGGATCAGCCg                       >  1:4230210/1‑150 (MQ=255)
                                                                                                                                ccAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGtt                                                                                           <  2:3782267/79‑1 (MQ=255)
                                                                                                                                ccAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGtt                                                                                           >  1:3782267/1‑79 (MQ=255)
                                                                                                                                ccAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGgc                    >  1:1765227/1‑150 (MQ=255)
                                                                                                                                ccAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGgc                    <  1:4726572/150‑1 (MQ=255)
                                                                                                                                    gTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGa                <  1:1782577/150‑1 (MQ=255)
                                                                                                                                     tGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGaa               <  1:2818028/150‑1 (MQ=255)
                                                                                                                                       cTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTTCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGc                                      <  1:1648926/125‑1 (MQ=255)
                                                                                                                                       cTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGc                                      >  2:1648926/1‑125 (MQ=255)
                                                                                                                                       cTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGc             >  2:5213188/1‑150 (MQ=255)
                                                                                                                                        tCGATGGCGGGCAGGCCCCGCCCGGCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGca            >  2:5404294/1‑150 (MQ=255)
                                                                                                                                        tCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGca            >  2:4306730/1‑150 (MQ=255)
                                                                                                                                          gATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGcaca          >  1:5450889/1‑150 (MQ=255)
                                                                                                                                           aTGGCGGGCAGGCCCCGCCCGGCATCGAGCAGATCAGGAAGTTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACAt         <  1:3913067/150‑1 (MQ=255)
                                                                                                                                           aTGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACAt         >  2:1974623/1‑150 (MQ=255)
                                                                                                                                              ccGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGACGCACATGTc      >  1:5037880/1‑150 (MQ=255)
                                                                                                                                                  gcgggGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACATGTCggtg  >  1:1634340/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGCATACAGGCCATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACATGTCGGTG  >  NC_002947/3003638‑3003933

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