Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 848,507 T→C 5.6% K21K (AAA→AAG PP_0730 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947848,5070TC 5.6% 312.8 / 6.5 107K21K (AAA→AAGPP_0730hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (44/57);  new base C (3/3);  total (47/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

Only 100 of 113 total aligned reads displayed.

AATTCGCCGCCTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAAGGTGCGCC  >  NC_002947/848361‑848653
                                                                                                                                                  |                                                                                                                                                  
aaTTCGCCGCCTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGttt                                                                                                                                                  >  4:114931/1‑149 (MQ=255)
aaTTCGCCGCCTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGttt                                                                                                                                                  <  6:140220/149‑1 (MQ=255)
aaTTCGCCGCCTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGtt                                                                                                                                                   <  7:18334/148‑1 (MQ=255)
aaTTCGCCGCCTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGt                                                                                                                                                    <  3:119618/147‑1 (MQ=255)
 aTTCGCCGCCTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGttt                                                                                                                                                  <  2:300954/148‑1 (MQ=255)
   tCGCCGCCTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGc                                                                                                                                               >  6:190980/1‑149 (MQ=255)
    cgccgcCTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCATGATCGAGATCGTCAGGAACGAGCCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCt                                                                                                                                              >  1:56374/1‑149 (MQ=255)
     gccgccTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGCGGCGTCGCGCTGATGAATGATGGTGACGTGGCAGCCTAGCCGTTTGGCTg                                                                                                                                             >  6:307410/1‑149 (MQ=255)
     gccgccTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTg                                                                                                                                             <  5:422084/149‑1 (MQ=255)
      ccgccTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGc                                                                                                                                            <  6:410309/149‑1 (MQ=255)
      ccgccTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGc                                                                                                                                            >  6:266747/1‑149 (MQ=255)
      ccgccTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGc                                                                                                                                            <  8:378460/149‑1 (MQ=255)
      ccgccTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGc                                                                                                                                            <  2:231944/149‑1 (MQ=255)
      ccgccTTCCAGATTATCGAAACGCAAATAAGCATCCACATAACGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGCAACGAGGCGTCGCGCTGATGAATGAAGATGACGTCGCCGCCAAGCAGTTTGGCTGc                                                                                                                                            >  4:7377/1‑149 (MQ=255)
       cgccTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGcc                                                                                                                                           <  5:9294/149‑1 (MQ=255)
       cgccTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAACGGCTTGGCTGcc                                                                                                                                           <  2:100138/149‑1 (MQ=255)
         ccTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGAGGCAGCCAAACGGCTTGGCTGccac                                                                                                                                         <  6:68701/149‑3 (MQ=255)
               aGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGc                                                                                                                                   >  3:177903/1‑149 (MQ=255)
               aGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGc                                                                                                                                   <  5:220783/149‑1 (MQ=255)
               aGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGAGTCGTCCTTGATCGAGATCGTCAGGAACGAGCCGTCGAGCTGATGAATGAAGTTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAga                                                                                                                                   >  5:335990/1‑148 (MQ=255)
                  ttATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCcag                                                                                                                                <  6:171699/149‑1 (MQ=255)
                    aTCGATACTCAAATAAGCATCCAGATATCGCAGGTAAGGCTCGAGACGGTAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCcagc                                                                                                                               <  8:105292/148‑1 (MQ=255)
                         tACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCTGTTTGGCTGCCTCAAAAGCCAGCAGGCCg                                                                                                                         <  1:253597/149‑1 (MQ=255)
                         tACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCg                                                                                                                         >  3:128217/1‑149 (MQ=255)
                           cGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCga                                                                                                                        <  3:251433/148‑2 (MQ=255)
                              aaaTAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGc                                                                                                                           <  8:377389/142‑1 (MQ=255)
                              aaaTAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGc                                                                                                                           >  7:377389/1‑142 (MQ=255)
                              aaaTAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTAc                                                                                                                    >  1:301074/1‑149 (MQ=255)
                                     cATACACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGGCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTg                                                                                                             <  8:341244/149‑1 (MQ=255)
                                      aTCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGt                                                                                                            >  6:128610/1‑149 (MQ=255)
                                          acaTATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTcgac                                                                                                        >  8:192667/1‑149 (MQ=255)
                                              atCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGaa                                                                                                     >  6:379657/1‑148 (MQ=255)
                                              atCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGaa                                                                                                     >  4:72006/1‑148 (MQ=255)
                                                   aGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCAc                                                                                               <  5:190980/149‑1 (MQ=255)
                                                     gTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTcga                                                                                                         >  1:58084/1‑137 (MQ=255)
                                                     gTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTcga                                                                                                         <  2:58084/137‑1 (MQ=255)
                                                     gTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGAAGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGCCGAAAACCACgt                                                                                             >  2:111179/1‑149 (MQ=255)
                                                       aaGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAGGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACgtgt                                                                                           >  8:226114/1‑149 (MQ=255)
                                                       aaGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGACAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGACGCTTACGGTGCTGTCGACGAAAACCACgtgt                                                                                           >  1:78194/1‑149 (MQ=255)
                                                        aGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACgtgt                                                                                           <  2:402442/148‑1 (MQ=255)
                                                          gCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGttt                                                                                         >  5:217407/1‑148 (MQ=255)
                                                           cTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCa                                                                                       <  5:142378/149‑1 (MQ=255)
                                                           cTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCa                                                                                       >  5:431840/1‑149 (MQ=255)
                                                              gagaCGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCAtt                                                                                     <  3:255935/148‑1 (MQ=255)
                                                              gagaCGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTg                                                                                    >  1:186472/1‑149 (MQ=255)
                                                                  cGGGAGGGGTCGTCCTTGATCTTGATCGGCAGGAACGAGGCGTCGAGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCATCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAAc                                                                                <  8:151369/149‑1 (MQ=255)
                                                                     gAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTc                                                                             <  5:56278/149‑1 (MQ=255)
                                                                     gAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTc                                                                             <  3:72006/149‑1 (MQ=255)
                                                                      aGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCg                                                                            <  2:301074/149‑1 (MQ=255)
                                                                       ggTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCgg                                                                           <  5:266747/149‑1 (MQ=255)
                                                                        gtgtCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGa                                                                          <  5:128610/149‑1 (MQ=255)
                                                                        gtgtCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGa                                                                          >  1:105829/1‑149 (MQ=255)
                                                                              tcCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTc                                                                                        <  1:350489/129‑1 (MQ=255)
                                                                              tcCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTc                                                                                        >  2:350489/1‑129 (MQ=255)
                                                                                cTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAgg                                                                  >  4:271526/1‑149 (MQ=255)
                                                                                cTTGATCGACATCGTCACGATCGAGGGGGCGCGCTGATGAAGGATCGTGACGCGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAAGCACGTGTTTCATTGGAACCTCGGATTACCAgg                                                                  <  4:100198/149‑1 (MQ=255)
                                                                                  tGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGa                                                                >  8:334225/1‑149 (MQ=255)
                                                                                   gATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTAc                                                                      >  8:203778/1‑142 (MQ=255)
                                                                                   gATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTAc                                                                      <  7:203778/142‑1 (MQ=255)
                                                                                   gATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGa                                                                <  4:128217/148‑1 (MQ=255)
                                                                                       gagaTCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGaaa                                                           <  2:186472/149‑1 (MQ=255)
                                                                                       gagaTCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGaaa                                                           <  3:159982/149‑1 (MQ=255)
                                                                                                aGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGt                                                  <  5:147626/149‑1 (MQ=255)
                                                                                                  gAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCAt                                                                                      >  7:317302/1‑111 (MQ=255)
                                                                                                  gAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCAt                                                                                      <  8:317302/111‑1 (MQ=255)
                                                                                                   aaCGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCATCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGc                                                <  1:226590/148‑1 (MQ=255)
                                                                                                    aCGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTg                                              <  3:88884/149‑1 (MQ=255)
                                                                                                    aCGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTg                                              <  6:372148/149‑1 (MQ=255)
                                                                                                          cGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCg                                        <  1:202712/149‑1 (MQ=255)
                                                                                                           gTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGGGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCg                                        >  8:330580/1‑148 (MQ=255)
                                                                                                              gcgcTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGaa                                      >  3:363911/1‑147 (MQ=255)
                                                                                                               cgcTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAACCCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAg                                   >  6:417468/1‑149 (MQ=255)
                                                                                                                  tgatgaATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCAt                                >  3:393801/1‑149 (MQ=255)
                                                                                                                    atgtatgaAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTg                              >  7:154263/1‑149 (MQ=255)
                                                                                                                    atgaatgaAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTg                              >  1:5982/1‑149 (MQ=255)
                                                                                                                     tgaatgaaGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTg                              <  7:184172/148‑1 (MQ=255)
                                                                                                                      gaatgaaGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGc                             >  2:91128/1‑148 (MQ=255)
                                                                                                                        atgaaGTTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGGAGGGCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCTGGGATGAAAGTAAATAGTGCTGTGCCCGAACATCATTGCGa                           <  6:335990/148‑1 (MQ=255)
                                                                                                                        atgaaGGTGACTTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaa                          <  4:177903/149‑1 (MQ=255)
                                                                                                                        atgaaGGTGACGTGGCATCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaa                          <  4:347518/149‑1 (MQ=255)
                                                                                                                        atgaaGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaa                          <  7:192667/149‑1 (MQ=255)
                                                                                                                        atgaaGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaa                          <  2:105829/149‑1 (MQ=255)
                                                                                                                         tgaaGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaag                         <  6:217407/149‑1 (MQ=255)
                                                                                                                         tgaaGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaag                         >  8:153018/1‑149 (MQ=255)
                                                                                                                          gaaGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaa                          >  6:365342/1‑147 (MQ=255)
                                                                                                                          gaaGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaa                          <  5:365342/147‑1 (MQ=255)
                                                                                                                              gTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGACCGAACAGCATTGCGaagcaagc                    >  7:163516/1‑149 (MQ=255)
                                                                                                                                tttgaggCAGCCAAACGGCTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaagcaagc                    >  7:143385/6‑147 (MQ=255)
                                                                                                                                tttgaggCAGCCAAACGGCTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaagcaagc                    <  8:143385/142‑1 (MQ=255)
                                                                                                                                gACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGaa                                      <  6:196786/129‑1 (MQ=255)
                                                                                                                                gACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGaa                                      >  5:196786/1‑129 (MQ=255)
                                                                                                                                 aCGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGTGAc                 <  7:58859/149‑1 (MQ=255)
                                                                                                                                    tGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGt              <  4:363911/149‑1 (MQ=255)
                                                                                                                                     ggCAGCCAAACGGCTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGATTCAAGCGCCCGTc             >  1:152093/1‑149 (MQ=255)
                                                                                                                                        aGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGaa                                      >  5:24076/1‑121 (MQ=255)
                                                                                                                                        aGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGaa                                      <  6:24076/121‑1 (MQ=255)
                                                                                                                                        aGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaag                         >  5:69583/1‑134 (MQ=255)
                                                                                                                                        aGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaag                         <  6:69583/134‑1 (MQ=255)
                                                                                                                                         gCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAAg         <  2:5982/149‑1 (MQ=255)
                                                                                                                                               ccGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTTAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAAGGTgcgc   <  2:56374/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AATTCGCCGCCTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAAGGTGCGCC  >  NC_002947/848361‑848653

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