Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 370,550 C→T 6.5% Q11* (CAG→TAG) ‡ dgcA → dimethylglycine dehydrogenase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947370,5500CT 6.5% 232.6 / 3.9 77Q11* (CAG→TAG) ‡dgcAdimethylglycine dehydrogenase subunit
Reads supporting (aligned to +/- strand):  ref base C (33/39);  new base T (2/3);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

GCTGCCCGCGGCAGCAAGATCCGCACCGTGGCCGAGCAGGTCTATGGCGGGTCGGAAGAAGGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTG  >  NC_002947/370406‑370684
                                                                                                                                                |                                                                                                                                      
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                                                                           <  1:188461/149‑1 (MQ=255)
                                                         gaagGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGc                                                                           <  5:387054/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                                                                <  8:148885/149‑1 (MQ=255)
                                                                    aGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGAcggcg                                                                <  3:91575/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                                                  >  4:254507/1‑149 (MQ=255)
                                                                                  ttACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGAcgaccgac                                                  >  2:317934/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                             <  8:351237/149‑1 (MQ=255)
                                                                                                       gagaTAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGaa                             <  4:96795/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                             <  1:115369/129‑1 (MQ=255)
                                                                                                                           gAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGTACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGaa                             >  2:115369/1‑129 (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             >  7:12182/4‑134 (MQ=255)
                                                                                                                                      tttgCCGATCTGGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGggcg             <  8:12182/131‑1 (MQ=255)
                                                                                                                                        aGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTAt                                            >  2:358128/1‑101 (MQ=255)
                                                                                                                                        aGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTAt                                            <  1:358128/101‑1 (MQ=255)
                                                                                                                                         gCCGATCTGGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGggcggc           >  6:294206/1‑133 (MQ=255)
                                                                                                                                         gCCGATCTGGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGggcggc           <  5:294206/133‑1 (MQ=255)
                                                                                                                                         gCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTAtgtg                                         <  7:289310/103‑1 (MQ=255)
                                                                                                                                         gCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTAtgtg                                         >  8:289310/1‑103 (MQ=255)
                                                                                                                                         gCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGa                              <  4:206826/114‑1 (MQ=255)
                                                                                                                                         gCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGa                              >  3:206826/1‑114 (MQ=255)
                                                                                                                                         gCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTg  >  7:392269/1‑142 (MQ=255)
                                                                                                                                         gCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTg  <  8:392269/142‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                      
GCTGCCCGCGGCAGCAAGATCCGCACCGTGGCCGAGCAGGTCTATGGCGGGTCGGAAGAAGGCCACGAAGATGGCCTGTTCGTTACAAAGCCGTTGTAAGCCGGAGATAGCGTCATGGCATTCGAAGCAATGTTCCAGCCGATCCAGATCGGCAAACTGACCATCCGCAACCGTGTGCTCAGCACCGCGCATGCCGAGGTCTACGCCACTGACGGCGGCATGACGACCGACCGCTATGTGAAGTACTACGAAGAGAAGGCCAAGGGCGGCATCGGCCTG  >  NC_002947/370406‑370684

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: