Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,098,367 A→G 5.7% N244S (AAC→AGC)  PP_2707 → exodeoxyribonuclease III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,098,3670AG 5.7% 209.0 / 2.5 71N244S (AAC→AGC) PP_2707exodeoxyribonuclease III
Reads supporting (aligned to +/- strand):  ref base A (30/36);  new base G (2/2);  total (33/38)
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.00e-01

GTACCCGGACGAGCGTGTCTACACCTTCTGGGACTATTTTCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAGATTGCATTCAGATTCGGCCCCG  >  NC_002947/3098221‑3098514
                                                                                                                                                  |                                                                                                                                                   
gTACCCGGACGAGCGTGTCTACACCTTCTGGGACTATTTTCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAAc                                                                                                                                                    >  1:184339/1‑148 (MQ=255)
 tACCCGGACGAGCGTGTCTACACCTTCTGGGACTATTTTCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCa                                                                                                                                                  <  4:14547/149‑1 (MQ=255)
   cccTGACGAGCGTGTCTACACCTTATGGGACTATTATCGTAACCACTGTGCCCTCAATGCCGGGTTGCGAATCTATCACCTGTTGCTCAATCTCGAGCTTGCGCCTTATCTGAAACGGGCTGGTGTCGATGCATGGGTGCGCAACCag                                                                                                                                                 <  4:167189/148‑1 (MQ=255)
   cccGGACGAGCGTGTCTACACCTTCTGGGACTATTTTCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCagg                                                                                                                                                >  5:156607/1‑149 (MQ=255)
   cccGGACGAGCGTGTCTACACCTTCTGGGACTATTTTCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCagg                                                                                                                                                >  4:48233/1‑149 (MQ=255)
    ccGGACGAGCGTGTCTACACCTTCTGGGACTATTTTCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCaggc                                                                                                                                               >  4:320740/1‑149 (MQ=255)
       gccGAGCGTGTCTACACCTTCTGGGACTATTTTCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCaggccaa                                                                                                                                            <  4:134694/147‑1 (MQ=255)
          gAGCGTGTCTACACCTTCTGGGACTATTTTCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCaggccaaggc                                                                                                                                         <  3:38876/149‑1 (MQ=255)
                gtCTACACCTTCTGGGACTATTTTCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGa                                                                                                                                   >  5:116691/1‑149 (MQ=255)
                   tACACCTTCTGGGACTATTTTCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCa                                                                                                                                >  7:209565/1‑149 (MQ=255)
                   tACACCTTCTGGGACTATTTTCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTCAAACGGGCGGGTGTGGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCa                                                                                                                                >  6:359660/1‑149 (MQ=255)
                    acacCTTCTGGGACTATTTTCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCAt                                                                                                                               >  5:209546/1‑149 (MQ=255)
                          tCTGGGACTATTTTCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGccc                                                                                                                         >  6:68285/1‑149 (MQ=255)
                                aCTATTTTCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTgg                                                                                                                   <  3:48233/149‑1 (MQ=255)
                                  tATTTTCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGAt                                                                                                                 >  1:198267/1‑149 (MQ=255)
                                    ttttCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATc                                                                                                                <  8:209565/148‑1 (MQ=255)
                                      ttCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGGGTCGATGCTTGGGTGCGCACCCCGGCAACGGCCAGCGCCCTGGCGCCCGAGTGGATCGAg                                                                                                             >  2:108309/1‑149 (MQ=255)
                                          tAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTgg                                                                                                         >  4:126202/1‑149 (MQ=255)
                                          tAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGACAGCGACCATGCTCCCGTCTGGATCGAGCTgg                                                                                                         <  2:184339/149‑1 (MQ=255)
                                            aCCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAGCCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTgg                                                                                                         <  1:374376/147‑1 (MQ=255)
                                            aCCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAGCCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTgg                                                                                                         <  1:380219/147‑1 (MQ=255)
                                            aCCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAGCCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTgg                                                                                                         >  2:380219/1‑147 (MQ=255)
                                            aCCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAGCCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTgg                                                                                                         >  2:374376/1‑147 (MQ=255)
                                            aCCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTggg                                                                                                        <  3:320740/148‑1 (MQ=255)
                                            aCCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTggg                                                                                                        <  6:116691/148‑1 (MQ=255)
                                            aCCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTgg                                                                                                         <  3:157332/147‑1 (MQ=255)
                                            aCCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTgg                                                                                                         >  4:157332/1‑147 (MQ=255)
                                             ccACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTa                                                                                                      <  6:363535/149‑1 (MQ=255)
                                             ccACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTa                                                                                                      <  6:156607/149‑1 (MQ=255)
                                              cACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTAc                                                                                                     >  3:86655/1‑149 (MQ=255)
                                              cACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTAc                                                                                                     >  4:192945/1‑149 (MQ=255)
                                                cTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGa                                                                                                                  <  3:307113/134‑1 (MQ=255)
                                                cTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGa                                                                                                                  >  4:307113/1‑134 (MQ=255)
                                                     cccGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGcgtcgt                                                                                              >  7:231798/1‑149 (MQ=255)
                                                          aaTGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCg                                                                                         <  5:164606/149‑1 (MQ=255)
                                                          aaTGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGGGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCAAGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCg                                                                                         <  2:105870/149‑1 (MQ=255)
                                                            tGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCTGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGt                                                                                       <  6:182221/149‑1 (MQ=255)
                                                                gggTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCa                                                                                     >  6:7503/1‑147 (MQ=255)
                                                                gggTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATa                                                                                   <  2:198267/149‑1 (MQ=255)
                                                                  gTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGt                                                                                 <  3:192945/149‑1 (MQ=255)
                                                                       gAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCGGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGtt                                                                            >  1:360220/1‑149 (MQ=255)
                                                                            gATCACCTGTTGCTCAATCGCGAGCTTGCGGCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTGAGc                                                                        <  2:95870/148‑1 (MQ=255)
                                                                             aTCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGca                                                                       >  4:66494/1‑148 (MQ=255)
                                                                                aCCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACAt                                                                    <  5:7503/148‑1 (MQ=255)
                                                                                     ttGCTCAATCGGGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCAAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAAc                                                                <  3:219329/147‑1 (MQ=255)
                                                                                     ttGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAAc                                                                >  5:335530/1‑147 (MQ=255)
                                                                                     ttGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGa                                                              <  4:112029/149‑1 (MQ=255)
                                                                                         tCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAgcgc                                                          >  7:97963/1‑149 (MQ=255)
                                                                                           aaTCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGt                                                                                 <  7:122204/124‑1 (MQ=255)
                                                                                           aaTCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGt                                                                                 >  8:122204/1‑124 (MQ=255)
                                                                                           aaTCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCg                                                        >  4:247367/1‑149 (MQ=255)
                                                                                             tCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCggg                                                      <  3:66494/149‑1 (MQ=255)
                                                                                               gcgAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGc                                                     <  2:360220/148‑1 (MQ=255)
                                                                                               gcgAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCt                                                    >  5:87472/1‑149 (MQ=255)
                                                                                                cgAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTg                                                   <  3:126202/149‑1 (MQ=255)
                                                                                                 gAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGa                                                  <  4:86655/149‑1 (MQ=255)
                                                                                                 gAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGa                                                  <  3:274925/149‑1 (MQ=255)
                                                                                                     ttGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGgcg                                              >  4:336196/1‑149 (MQ=255)
                                                                                                       gcgcCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCt                                            >  7:108028/1‑149 (MQ=255)
                                                                                                                   aaaCGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTa                                   <  6:209546/146‑1 (MQ=255)
                                                                                                                           gggTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAg                        <  5:200108/149‑1 (MQ=255)
                                                                                                                                tCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAGATTGc                   <  3:247367/149‑1 (MQ=255)
                                                                                                                                tCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAGATTGc                   <  5:68285/149‑1 (MQ=255)
                                                                                                                                     gCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAGATTGCATTc               <  8:108028/148‑1 (MQ=255)
                                                                                                                                     gCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGAGCTCCGACAGTAGCCGACTTCAGATTGCAtt                >  8:194247/1‑147 (MQ=255)
                                                                                                                                       aTGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAGATTGCATTCAGa            <  6:335530/149‑1 (MQ=255)
                                                                                                                                         gggTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAGATTGCATTCAGAtt          >  6:156932/1‑149 (MQ=255)
                                                                                                                                           gTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAGATTGCATTCAGATTCg        <  1:108309/149‑1 (MQ=255)
                                                                                                                                               gcAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAGATTGCATTCAGATTCGGccc    >  2:50229/1‑149 (MQ=255)
                                                                                                                                                cAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAGATTGCATTCAGATTCGGc      <  7:309364/146‑1 (MQ=255)
                                                                                                                                                cAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAGATTGCATTCAGATTCGGc      >  8:309364/1‑146 (MQ=255)
                                                                                                                                                 aaCCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAGATTGCATTCAGATTCGGCCCCg  <  8:231798/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTACCCGGACGAGCGTGTCTACACCTTCTGGGACTATTTTCGTAACCACTGGGCCCGCAATGCCGGGTTGCGAATCGATCACCTGTTGCTCAATCGCGAGCTTGCGCCTTATCTGAAACGGGCGGGTGTCGATGCATGGGTGCGCAACCAGGCCAAGGCCAGCGACCATGCGCCCGTCTGGATCGAGCTGGGTACGCGTCGTGCTCGGTCATAGTCAGTTCAGCACATGAACGAGCGCCGGGCTGAGGCGCTCCGACAGTAGCCGACTTCAGATTGCATTCAGATTCGGCCCCG  >  NC_002947/3098221‑3098514

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: ‑