Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,001,876 C→G 6.9% A186G (GCC→GGC)  PP_2625 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,001,8760CG 6.9% 304.6 / 12.9 101A186G (GCC→GGC) PP_2625hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (45/49);  new base G (2/5);  total (47/54)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.15e-01

Only 100 of 102 total aligned reads displayed.

ATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACATGTCGGTGGTAGCGCA  >  NC_002947/3001728‑3002019
                                                                                                                                                    |                                                                                                                                               
aTCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATTGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGcc                                                                                                                                                <  1:2255909/150‑1 (MQ=255)
 tctcCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCa                                                                                                                                               <  1:1712399/150‑1 (MQ=255)
   tcCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATc                                                                                                                                             <  1:1042054/150‑1 (MQ=255)
   tcCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATc                                                                                                                                             >  1:1987879/1‑150 (MQ=255)
   tcCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATc                                                                                                                                             >  1:1432306/1‑150 (MQ=255)
   tcCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATc                                                                                                                                             >  2:1532391/1‑150 (MQ=255)
     ccGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGGCGGGCAGGCCCCGCCCGGCATCGa                                                                                                                                           <  2:1245473/150‑1 (MQ=255)
     ccGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGa                                                                                                                                           >  2:2295449/1‑150 (MQ=255)
     ccGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGa                                                                                                                                           >  2:1469708/1‑150 (MQ=255)
     ccGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGa                                                                                                                                           >  1:1856880/1‑150 (MQ=255)
     ccGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGAAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGa                                                                                                                                           >  2:1025550/1‑150 (MQ=255)
     ccGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCGACCCGTGAGAAGCCCTGGAAAGTGTTAGATAACTTCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGa                                                                                                                                           >  1:1742091/1‑150 (MQ=255)
             aGCTTATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGGCGGGCAGGCCCCGCCCGGCATCGAGCacttgg                                                                                                                                   <  2:569007/150‑6 (MQ=38)
               cTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCAGCGAGCAGATCAgg                                                                                                                                 >  1:1977221/1‑150 (MQ=255)
                 gATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGaa                                                                                                                               >  2:119101/1‑150 (MQ=255)
                  aTTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAg                                                                                                                              <  1:785166/150‑1 (MQ=255)
                  aTTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAg                                                                                                                              >  1:1578587/1‑150 (MQ=255)
                    tGCCAGGCATCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGATCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCt                                                                                                                            <  2:1983231/150‑1 (MQ=255)
                       cAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGAt                                                                                                                         >  2:118164/1‑150 (MQ=255)
                            gTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGCTCAGGAGGCTGATCGGTg                                                                                                                    >  2:39682/1‑150 (MQ=255)
                            gTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTg                                                                                                                    <  2:1469081/150‑1 (MQ=255)
                             tCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCATGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGc                                                                                                                   >  2:188698/1‑150 (MQ=255)
                               cGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGccc                                                                                                                 >  2:1442674/1‑150 (MQ=255)
                                 aTCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAg                                                                                                               >  1:2372678/1‑150 (MQ=255)
                                 aTCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAg                                                                                                               <  2:2419015/150‑1 (MQ=255)
                                  tCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGc                                                                                                              <  1:1532391/150‑1 (MQ=255)
                                      ggCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTAt                                                                                                          >  1:216570/1‑150 (MQ=255)
                                       gCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATc                                                                                                         <  1:139278/150‑1 (MQ=255)
                                       gCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATc                                                                                                         >  2:2152952/1‑150 (MQ=255)
                                       gCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATc                                                                                                         <  2:2141062/150‑1 (MQ=255)
                                       gCAAGCTGGGTCAAGGTACGCCCGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATc                                                                                                         <  2:1199379/150‑1 (MQ=255)
                                          aGCTGGGTTAAGGTTCGCCTGAGCGGGGAGGGCTTGTATACCCGGGAGACGCCCTGGCAAGTGATGGATAACGGCCAAGTGCTCGATGCCGGTCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAg                                                                                                      <  2:588790/150‑1 (MQ=255)
                                            cTGGGCCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGt                                                                                                    <  2:1110795/150‑1 (MQ=255)
                                              gggTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAg                                                                                                                                          >  1:177043/1‑110 (MQ=255)
                                              gggTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAg                                                                                                                                          <  2:177043/110‑1 (MQ=255)
                                              gggTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTtg                                                                                                  <  1:2152952/150‑1 (MQ=255)
                                               ggTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTTATCGGTGCCCAGCGTATCGAGGTtgt                                                                                                 <  1:2037240/150‑1 (MQ=255)
                                                 tCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAgg                                                                                                                                 <  2:1743860/116‑1 (MQ=255)
                                                 tCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAgg                                                                                                                                 >  1:1743860/1‑116 (MQ=255)
                                                 tCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTgtga                                                                                               <  2:2226955/150‑1 (MQ=255)
                                                 tCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCCGGCGGGTGGTGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTgtga                                                                                               <  2:64864/150‑1 (MQ=255)
                                                 ctaaGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTgtga                                                                                               <  1:2355505/148‑1 (MQ=255)
                                                   aaGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTgtgagt                                                                                             <  2:1432306/150‑1 (MQ=255)
                                                   aaGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGGTgtgagt                                                                                             >  2:1956177/1‑150 (MQ=255)
                                                        aCGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCa                                                                                        <  2:792090/150‑1 (MQ=255)
                                                            cTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGc                                                                                    >  1:1604676/1‑150 (MQ=255)
                                                            cTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGc                                                                                    >  2:1361052/1‑150 (MQ=255)
                                                               aGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGtt                                                                                                   >  1:1159085/1‑132 (MQ=255)
                                                               aGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGtt                                                                                                   <  2:1159085/132‑1 (MQ=255)
                                                                gCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAgg                                                                                >  2:775801/1‑150 (MQ=255)
                                                                    ggAGGGCTGGTCTCCCCGTGAGAAGCTCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGTGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGa                                                                            <  2:499369/150‑1 (MQ=255)
                                                                    ggAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCa                                                                                                                                   <  2:2146477/95‑1 (MQ=255)
                                                                    ggAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCa                                                                                                                                   >  1:2146477/1‑95 (MQ=255)
                                                                    ggAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGa                                                                            >  1:17375/1‑150 (MQ=255)
                                                                    ggAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCAATGCCGGGCGGGGCCTGCCCGCCATCGAGc                                                                                                                                         <  1:2017844/89‑1 (MQ=255)
                                                                    ggAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCAATGCCGGGCGGGGCCTGCCCGCCATCGAGc                                                                                                                                         >  2:2017844/1‑89 (MQ=255)
                                                                        ggCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTAGATGCCGGGGGGGGTCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTCCGTGATCACGGCGAGGTCGACggg                                                                        <  1:1057038/150‑1 (MQ=255)
                                                                             gTCTACCCGTGAGACACCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGAc                                                                   >  2:2299248/1‑150 (MQ=255)
                                                                                tACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCATGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCtt                                                                >  1:2358185/1‑150 (MQ=255)
                                                                                 aCCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCCGCGTAACGAGGTTGCAAGTGAGCACGGCGAGGTCGACGGGCTGACCTTc                                                               <  1:1871731/150‑1 (MQ=255)
                                                                                   ccGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGa                                                             >  2:1959644/1‑150 (MQ=255)
                                                                                   ccGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGa                                                             <  1:1959644/150‑1 (MQ=255)
                                                                                    cGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGAt                                                            <  1:1802890/150‑1 (MQ=255)
                                                                                     gTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATg                                                           <  1:2338568/150‑1 (MQ=255)
                                                                                      tGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATgc                                                          <  2:1856880/150‑1 (MQ=255)
                                                                                          ccGCCCTGGCAAGTGAGGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCaa                                                      <  1:119101/149‑1 (MQ=255)
                                                                                          aCGCCCTGGCAAGTGGTGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCaa                                                      >  1:621637/1‑150 (MQ=255)
                                                                                                tGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCTCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTg                                                >  1:1240773/1‑150 (MQ=255)
                                                                                                tGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTg                                                <  2:1578587/150‑1 (MQ=255)
                                                                                                 ggCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTgg                                               >  1:162856/1‑150 (MQ=255)
                                                                                                  gCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGc                                              <  2:235020/150‑1 (MQ=255)
                                                                                                    aaGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACg                                                                          <  2:124585/120‑1 (MQ=255)
                                                                                                    aaGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGccc                                            >  1:1443807/1‑150 (MQ=255)
                                                                                                    aaGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGACTTCGACg                                                                          >  1:124585/1‑120 (MQ=255)
                                                                                                     aGTGGTGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGcccc                                           <  2:621637/150‑1 (MQ=255)
                                                                                                       tGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCt                                         <  2:2372678/150‑1 (MQ=255)
                                                                                                        gATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTACCACTCGTGACGCAATTCCTGGCCCCCtt                                        <  2:1402026/150‑1 (MQ=255)
                                                                                                           ggAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGTGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGAt                                     >  1:551975/1‑150 (MQ=255)
                                                                                                             agaACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCa                                   >  1:1402097/1‑150 (MQ=255)
                                                                                                                   gCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCggg                             <  1:775801/150‑1 (MQ=255)
                                                                                                                    ccAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCt                                         <  1:1815188/137‑1 (MQ=255)
                                                                                                                    ccAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCt                                         >  2:1815188/1‑137 (MQ=255)
                                                                                                                          gCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAg                      <  1:118164/150‑1 (MQ=255)
                                                                                                                            tCGATGGCGGGCAGGCCCCGCCCGGCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGTGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGgc                            >  2:1765793/1‑142 (MQ=38)
                                                                                                                            tCGATGGCGGGCAGGCCCCGCCCGGCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGTGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGgc                            <  1:1765793/142‑1 (MQ=38)
                                                                                                                            tCGATGGCGGGCAGGCCCCGCCCGGCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGca                    >  2:1873969/1‑150 (MQ=255)
                                                                                                                            tCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGc                                                                                                              >  1:658170/1‑60 (MQ=255)
                                                                                                                            tCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGc                                                                                                              <  2:658170/60‑1 (MQ=255)
                                                                                                                            tCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGca                    >  2:2072937/1‑150 (MQ=255)
                                                                                                                              gATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGcaca                  >  1:1720486/1‑150 (MQ=255)
                                                                                                                                tggcgggcaggCCCCGCCCGGCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGgc                            <  1:1755664/135‑1 (MQ=255)
                                                                                                                                      gcgggGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCTCCTTGATCAGCCGGGCGCGAAGCACATGTCggtg          <  1:188698/150‑1 (MQ=255)
                                                                                                                                       cgggGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACATGTCggtgg         <  2:1443807/150‑1 (MQ=255)
                                                                                                                                            ccTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGAGCCCCTTGATCAGCCGGGCGCGAAGCACATGTCGGTGGTAgcg    >  1:513094/1‑150 (MQ=255)
                                                                                                                                             cTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGcac                   <  2:509748/134‑1 (MQ=255)
                                                                                                                                             cTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGcac                   >  1:509748/1‑134 (MQ=255)
                                                                                                                                              tGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACATGTCGGTGGTAGCGCa  <  2:1402097/150‑1 (MQ=255)
                                                                                                                                              tGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACATGTCGGTGGTAGCGCa  <  1:2299248/150‑1 (MQ=255)
                                                                                                                                              tGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGACGCGAAGCACATGTCGGTGGTa       <  1:606440/145‑1 (MQ=255)
                                                                                                                                              tGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGACGCGAAGCACATGTCGGTGGTa       >  2:606440/1‑145 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
ATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACATGTCGGTGGTAGCGCA  >  NC_002947/3001728‑3002019

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: