Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 370,551 A→G 6.5% Q11R (CAG→CGG) ‡ dgcA → dimethylglycine dehydrogenase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947370,5510AG 6.5% 222.6 / 3.0 76Q11R (CAG→CGG) ‡dgcAdimethylglycine dehydrogenase subunit
Reads supporting (aligned to +/- strand):  ref base A (33/38);  new base G (2/3);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGCCCGCGGCAGCAAGATCCGCACCGTGGCCGAGCAGGTCTATGGCGGGTCGGAAGAAGGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGC  >  NC_002947/370406‑370699
                                                                                                                                                 |                                                                                                                                                    
gCTGCCCGCGGCAGCAAGATCCGCACCGTGGCCGAGCAGGTCTATGGCGGGTCGGAAGAAGGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCa                                                                                                                                                      <  8:151606/146‑1 (MQ=255)
gCTGCCCGCGGCAGCAAGATCCGCACCGTGGCCGAGCAGGTCTATGGCGGGTCGGAAGAAGGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGa                                                                                                                                                    <  7:268590/148‑1 (MQ=255)
gCTGCCCGCGGCAGAAAGATCCGCACCGTGGCCGAGCAGGTCTATGGCGGGTCGGAAGAAGGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGAt                                                                                                                                                   <  8:273533/149‑1 (MQ=255)
 cTGCCCGCGGCAGCAAGATCCGCACCGTGGCCGAGCAGGTCTATGGCGGGTCGGAAGAAGGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATc                                                                                                                                                  <  3:335795/149‑1 (MQ=255)
    cccGCGGCAGCAAGATCCGCACCGTGGCCGAGCAGGTCTATGGCGGGTCGGAAGAAGGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGc                                                                                                                                               <  2:349231/149‑1 (MQ=255)
     ccGCGGCAGCAAGATCCGCACCGTGGCCGAGCAGGTCTATGGCGGGTCGGAAGAAGGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCTGGATCGGct                                                                                                                                              <  2:259173/149‑2 (MQ=255)
            agcaAGATCCGCACCGTGGCCGAGCAGGTCTATGGCGGGTCGGAAGAAGGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGAc                                                                                                                                       <  4:256088/149‑1 (MQ=255)
              caAGATCCGCACCGTGGCCGAGCAGGTCTATGGCGGGTCGGAAGAAGGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCa                                                                                                                                     >  1:260560/1‑149 (MQ=255)
               aaGATCCGCACCGTGGCCGAGCAGGTCTATGGCGGGTCGGAAGAAGGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCAt                                                                                                                                    >  5:405513/1‑149 (MQ=255)
                aGATCCGCACCGTGGCCGAGCAGGTCTATGGCGGGTCGGAAGAAGGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATc                                                                                                                                   >  8:103688/1‑149 (MQ=255)
                   tCCGCACCGTGGCCGAGCAGGTCTATGGCGGGTCGGAAGAAGGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCg                                                                                                                                 >  5:215701/1‑148 (MQ=255)
                       cACCGTGGCCGAGCAGGTCTATGGCGGGTCGGAAGAAGGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAAcc                                                                                                                            <  6:405513/149‑1 (MQ=255)
                           gTGGCCGAGCAGGTCTATGGCGGGTCGGAAGAAGGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCgtgt                                                                                                                        <  8:370691/149‑1 (MQ=255)
                            tGGCCGAGCAGGTCTATGGCGGGTCGGAAGAAGGTCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGtgtg                                                                                                                       <  1:124225/149‑1 (MQ=255)
                              gCCGAGCAGGTCTATGGCGGGTCGGAAGAAGGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGc                                                                                                                      <  7:355128/148‑1 (MQ=255)
                                       gTCTATGGCGGGTCGGAAGAAGGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACcgc                                                                                                            >  2:338583/1‑149 (MQ=255)
                                           aTGGCGGGTCGGAAGAAGGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCa                                                                                                                <  6:181054/141‑1 (MQ=255)
                                           aTGGCGGGTCGGAAGAAGGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCa                                                                                                                >  5:181054/1‑141 (MQ=255)
                                            tGGCGGGTCGGAAGAAGGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCa                                                                                                         <  7:103688/147‑1 (MQ=255)
                                            tGGCGGGTCGGAAGAAGGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATg                                                                                                       >  7:148885/1‑149 (MQ=255)
                                            tGGCGGGTCGGAAGAAGGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATg                                                                                                       <  1:375296/149‑1 (MQ=255)
                                                         gaagGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGc                                                                                          <  5:387054/149‑1 (MQ=255)
                                                         gaagGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGc                                                                                          <  1:188461/149‑1 (MQ=255)
                                                         gaagGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGc                                                                                          >  5:183580/1‑149 (MQ=255)
                                                           agGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCa                                                                                        <  6:60016/149‑1 (MQ=255)
                                                               cACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGCCCATCCGCAACCGTGTGCTCAGCACAGAGCATGCCGAGGTCTACGCAACTGa                                                                                    >  8:347735/1‑149 (MQ=255)
                                                               cACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGa                                                                                    >  1:132746/1‑149 (MQ=255)
                                                               cACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCAGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGa                                                                                    >  7:351237/1‑149 (MQ=255)
                                                                aCGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGAc                                                                                   <  7:36469/149‑1 (MQ=255)
                                                                    aGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGAcggcg                                                                               <  3:91575/149‑1 (MQ=255)
                                                                    aGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGAcggcg                                                                               <  8:148885/149‑1 (MQ=255)
                                                                        ggCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCACCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCAt                                                                           >  3:185534/1‑149 (MQ=255)
                                                                          ccTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATga                                                                         >  5:143241/1‑149 (MQ=255)
                                                                           cTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATgac                                                                        <  4:405796/149‑1 (MQ=255)
                                                                           cTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATgac                                                                        >  5:329909/1‑149 (MQ=255)
                                                                                 gTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACcg                                                                                                             >  4:392381/1‑106 (MQ=255)
                                                                                 gTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACcg                                                                                                             <  3:392381/106‑1 (MQ=255)
                                                                                  ttACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGAcgaccgac                                                                 >  2:317934/1‑149 (MQ=255)
                                                                                  ttACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGAcgaccgac                                                                 >  4:254507/1‑149 (MQ=255)
                                                                                    aCAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGaccgaccg                                                               >  5:334559/1‑149 (MQ=255)
                                                                                      aaaGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGCCGACCGACCGc                                                              >  8:289577/1‑148 (MQ=255)
                                                                                      aaaGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGc                                                              >  8:295187/1‑148 (MQ=255)
                                                                                      aaaGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGc                                                              >  6:419975/1‑148 (MQ=255)
                                                                                        aGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTAt                                                           <  6:143241/149‑1 (MQ=255)
                                                                                          ccGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTAtgt                                                         <  3:254507/149‑1 (MQ=255)
                                                                                                tAAGCCGGAGATATCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAg                                                                                                  <  7:176177/102‑1 (MQ=255)
                                                                                                tAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAg                                                                                                  >  8:176177/1‑102 (MQ=255)
                                                                                                tAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGta                                                   >  8:105658/1‑149 (MQ=255)
                                                                                                  aGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGtac                                                  <  2:132746/148‑1 (MQ=255)
                                                                                                   gCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGtacta                                                <  1:338583/149‑1 (MQ=255)
                                                                                                     cGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACg                                              <  6:215701/149‑1 (MQ=255)
                                                                                                       gagaTAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGaa                                            <  4:96795/149‑1 (MQ=255)
                                                                                                       gagaTAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGaa                                            <  8:351237/149‑1 (MQ=255)
                                                                                                             gCGTCTTGGATTTGGATGAAATGTTCCAGCCTATCCATATGGGAAAATTTACCATCGGCAACCGTTTGCTCAGAACGGCGATTGCGTAGGAGTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGaa                                       <  1:317934/148‑1 (MQ=255)
                                                                                                             gCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGaa                                       <  6:329909/148‑1 (MQ=255)
                                                                                                                   tGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTAtgtg                                                        >  6:76877/1‑125 (MQ=255)
                                                                                                                        ttCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGggcgg                           >  7:204141/1‑149 (MQ=255)
                                                                                                                           gAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGTACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGaa                                            >  2:115369/1‑129 (MQ=255)
                                                                                                                           gAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGTACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGaa                                            <  1:115369/129‑1 (MQ=255)
                                                                                                                            aaGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATc                       >  6:47046/1‑149 (MQ=255)
                                                                                                                            aaGCAATGTTCCAGCCGATACAGATCGGCAAACTGACAATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCa                         >  5:43177/1‑147 (MQ=255)
                                                                                                                             aGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCg                      <  2:260560/149‑1 (MQ=255)
                                                                                                                              gCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCgg                     <  6:183580/149‑1 (MQ=255)
                                                                                                                                aaTGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTCTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAagag                                         >  6:387835/1‑127 (MQ=255)
                                                                                                                                 aTGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCt                  <  7:105658/149‑1 (MQ=255)
                                                                                                                                      tttgCCGATCTGGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGggcg                            <  8:12182/131‑1 (MQ=255)
                                                                                                                                      tttgCCGATCTGGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGggcg                            >  7:12182/4‑134 (MQ=255)
                                                                                                                                        aGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTAt                                                           >  2:358128/1‑101 (MQ=255)
                                                                                                                                        aGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTAt                                                           <  1:358128/101‑1 (MQ=255)
                                                                                                                                         gCCGATCTGGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGggcggc                          <  5:294206/133‑1 (MQ=255)
                                                                                                                                         gCCGATCTGGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGggcggc                          >  6:294206/1‑133 (MQ=255)
                                                                                                                                         gCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTAtgtg                                                        <  7:289310/103‑1 (MQ=255)
                                                                                                                                         gCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTAtgtg                                                        >  8:289310/1‑103 (MQ=255)
                                                                                                                                         gCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGa                                             >  3:206826/1‑114 (MQ=255)
                                                                                                                                         gCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGa                                             <  4:206826/114‑1 (MQ=255)
                                                                                                                                         gCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTg                 >  7:392269/1‑142 (MQ=255)
                                                                                                                                         gCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTg                 <  8:392269/142‑1 (MQ=255)
                                                                                                                                                 aGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTggcggc  <  7:295187/149‑1 (MQ=255)
                                                                                                                                                 aGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTggcggc  <  7:289577/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCTGCCCGCGGCAGCAAGATCCGCACCGTGGCCGAGCAGGTCTATGGCGGGTCGGAAGAAGGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTGGCCATCTGTGGCGGC  >  NC_002947/370406‑370699

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