Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,870,190 A→C 8.0% D113A (GAC→GCC)  cobD → Subunit of adenosylcobinamide‑phosphate synthase beta component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,870,1900AC 8.0% 169.8 / 7.8 64D113A (GAC→GCC) cobDSubunit of adenosylcobinamide‑phosphate synthase beta component
Reads supporting (aligned to +/- strand):  ref base A (29/29);  new base C (3/2);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.66e-01

GTGGTGCCGCTGACCCTGGTGGCGCTGGTCCTTTCGTGGCTGCCGTATATCGGCTGGCTGGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGACCGAGTCGGTGCTGGAGAACGGCAGCGACGCAGTGTTCGCCGCGTTGTTCTGGTTTGTGGTG  >  NC_002947/1870042‑1870338
                                                                                                                                                    |                                                                                                                                                    
gtggggCCGCTGACCCTGGTGGCGCTGGTCCTTTCGTGGCTGCCGTATCACGGCTGGCTGGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGAc                                                                                                                                                     <  1:73372/150‑1 (MQ=255)
    tgCCGCTGACCCTGGTGGCGCTGGTCCTTTCGTGGCTGCCGTATATCGGCTGGCTGGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTgg                                                                                                                                                 >  2:1201948/1‑150 (MQ=255)
     gCCGCTGACCCTGGTGGCGCTGGTCCTTTCGTGGCTGCCGTATATCGGCTGGCTGGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTgg                                                                                                                                                 <  1:116516/149‑1 (MQ=255)
     gCCGCTGACCCTGGTGGCGCTGGTCCTTTCGTGGCTGCCGTATATCGGCTGGCTGGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTgg                                                                                                                                                 >  2:116516/1‑149 (MQ=255)
       cGCTGACCCTGGTGGCGCTGGTCCTTTCGTGGCTGCCGTATATCGGCTGGCTGGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCgg                                                                                                                                              >  2:181963/1‑150 (MQ=255)
        gCTGACCCTGGTGGCGCTGGTCCTTTCGTGGCTGCCGTATATCGGCTGGCTGGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGa                                                                                                                                             <  2:1205578/150‑1 (MQ=255)
         cTGACCCTGGTGGCGCTGGTCCTTTCGTGGCTGCCGTATATCGGCTGGCTGGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGTCGCCCTGGCgcaa                                                                                                                                            <  1:715926/150‑4 (MQ=255)
                tggtggCGCTGGTCCTTTCGTGGCTGCCGTATATCGGCTGGCTGGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAggcgcggc                                                                                                                                     <  1:690364/150‑1 (MQ=255)
                      cgcTGGTCCTTTCGTGGCTGCCGTATATCGGCTGGCTGGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCggg                                                                                                                               <  1:181963/150‑1 (MQ=255)
                             ccTTTCGTGGCTGCCGTATATCGGCTGGCTGGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTa                                                                                                                        <  1:4377/150‑1 (MQ=255)
                               tttCGTGGCTGCCGTATATCGGCTGGCTGGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTAcc                                                                                                                      >  2:421740/1‑150 (MQ=255)
                               tttCGTGGCTGCCGTATATCGGCTGGCTGGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTAcc                                                                                                                      <  1:421740/150‑1 (MQ=255)
                                   gTGGCTGCCGTATATCGGCTGGCTGGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTACTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGt                                                                                                                  >  1:1065514/1‑150 (MQ=255)
                                          ccGTATATCGGCTGGCTGGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAggcgcggc                                                                                                                                     >  2:573860/1‑124 (MQ=255)
                                          ccGTATATCGGCTGGCTGGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAggcgcggc                                                                                                                                     <  1:573860/124‑1 (MQ=255)
                                            gTATATCGGCTGGCTGGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTAcc                                                                                                                      >  2:813445/1‑137 (MQ=255)
                                            gTATATCGGCTGGCTGGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTAcc                                                                                                                      <  1:813445/137‑1 (MQ=255)
                                            gTATATCGGCTGGCTGGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGcc                                                                                                             >  1:1081713/1‑146 (MQ=255)
                                            gTATATCGGCTGGCTGGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGcc                                                                                                             <  2:1081713/146‑1 (MQ=255)
                                            gTATATCGGCTGGCTGGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGa                                                                                                         >  2:1139357/1‑150 (MQ=255)
                                             tataTCGGCTGGCTGGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGAAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGaa                                                                                                        >  2:347805/1‑150 (MQ=255)
                                              atatCGGCTGGCTGGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGaaa                                                                                                       >  1:128769/1‑150 (MQ=255)
                                                 tCGGCTGGCTGGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAAcgc                                                                                                    >  2:718483/1‑150 (MQ=255)
                                                   ggctggctGGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAAcgcgc                                                                                                  <  1:1201948/150‑1 (MQ=255)
                                                    gctggctgGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAAcgcgcg                                                                                                 >  2:727148/1‑150 (MQ=255)
                                                      tggctggTTGCTGTGCTGGCGCTTTACTGCGTCATGGGGCTGCGCAGCCTCGGCAAACATGTGCTGCCGGCGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAg                                                                                               >  2:140423/1‑150 (MQ=255)
                                                         ctggTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTg                                                                                            >  1:984087/1‑150 (MQ=255)
                                                           ggTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGGGTCAGGGGGAACTGTCTGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGa                                                                                          >  1:445278/1‑150 (MQ=255)
                                                            gTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGAc                                                                                         <  1:700182/150‑1 (MQ=255)
                                                            gTTGATGTGCTGGCGCTTTACTGCGCCGTGGAGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCATTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGAc                                                                                         >  2:286623/1‑150 (MQ=255)
                                                              tGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGa                                                                                       <  2:128769/150‑1 (MQ=255)
                                                                        gcgcTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTg                                                                             >  1:512151/1‑150 (MQ=255)
                                                                               aCTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAg                                                                      <  2:15598/150‑1 (MQ=255)
                                                                                            gCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGACCGAGTc                                                         >  2:899239/1‑150 (MQ=255)
                                                                                              tGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGACCGAGTCgg                                                       <  1:1100692/150‑1 (MQ=255)
                                                                                                 gcAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAAGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGACCGAGTCGGTGc                                                    >  1:479738/1‑150 (MQ=255)
                                                                                                   aGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGACCGAGTCGGTGCTg                                                  <  2:86005/150‑1 (MQ=255)
                                                                                                   aGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGACCGAGTCGGTGCTg                                                  <  1:286623/150‑1 (MQ=255)
                                                                                                    gCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGACCGAGTCGGTGCTgg                                                 <  1:793782/150‑1 (MQ=255)
                                                                                                           cGAGCATGTGCTGCGGGTGGCCAATGGATCGCGCGAGGGCGGCCTGGCGGAGGCGCGGCGCGGGGTTGGCTACCTGGTAGGCCGCGAAACGCGAGCGCTGGACGGGCCCGCGCTGGCACGAGCGGCGACAGGGTCGGTGCTGGAGAACgg                                          >  2:563335/1‑150 (MQ=255)
                                                                                                              gCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGACCGAGTCGGTGCTGGAGAACGGCAg                                       >  1:97448/1‑150 (MQ=255)
                                                                                                                  gtgCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGACCGAGTCGGTGCTGGAGAACGGCAGCGAc                                   <  1:727148/150‑1 (MQ=255)
                                                                                                                     ctgcCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGACCGAGTCGGTGCTGGAGAACGGCAGCGACGCa                                <  1:718483/150‑1 (MQ=255)
                                                                                                                         cGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGAcc                                                              <  1:199783/116‑1 (MQ=255)
                                                                                                                         cGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGAcc                                                              >  2:199783/1‑116 (MQ=255)
                                                                                                                           gTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGAc                                                               <  1:522110/113‑1 (MQ=255)
                                                                                                                           gTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGAc                                                               >  2:522110/1‑113 (MQ=255)
                                                                                                                           gTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGACCGAGTCGGTGCTGGAGAACGGCAGCGACGCAGTGTTc                          <  1:900373/150‑1 (MQ=255)
                                                                                                                               ccAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAAcgc                                                                                                    <  1:678981/72‑1 (MQ=255)
                                                                                                                               ccAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAAcgc                                                                                                    >  2:678981/1‑72 (MQ=255)
                                                                                                                                cAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGACCGAGTCGGTGCTGGAGAACGGCAGCGACGCAGTGTTcgccgc                     >  1:715216/1‑150 (MQ=255)
                                                                                                                                 aaTGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGACCGAGTCGGTGCTGGAGAACGGc                                         >  1:1123894/1‑129 (MQ=255)
                                                                                                                                 aaTGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGACCGAGTCGGTGCTGGAGAACGGc                                         <  2:1123894/129‑1 (MQ=255)
                                                                                                                                 aaTGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGACGCGAAACGCGGGGGCTGGACGAGCCCGCGGAGGGCCTAGCGGCCACTGAGTTGGTTCGGTAGAACGGCATCGTCGCAGTGTTcgccgcc                    >  2:1156587/1‑149 (MQ=255)
                                                                                                                                      cctccgcCAGGTCGCCCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCTCGAGCGGCGACCGAGTCGGTGCTGGAGAACGGCAGCGACGCAGTGTTAGACGCgtagtt               >  2:526378/5‑150 (MQ=255)
                                                                                                                                      cTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGACCGAGTCGGTGCTGGAGAACGGCAGCGACGCAGTGTTCGCCGCgttgtt               >  1:444257/1‑150 (MQ=255)
                                                                                                                                      cTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGACCGAGTCGGTGCTGGAGAACGGCAGCGACGCAGTGTTCGCCGCgttgtt               <  1:140423/150‑1 (MQ=255)
                                                                                                                                       ttGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGCTGTCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGACCGAGTCGGTGCTGGAGAACGGCAGCGACGCAGTGTTCGCCGCGTTGTTc              <  1:1156587/150‑1 (MQ=255)
                                                                                                                                        tccgcCAGGTCGCCCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGTTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGAc                                                               <  2:1198978/98‑1 (MQ=255)
                                                                                                                                        tccgcCAGGTCGCCCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGAc                                                               >  1:1198978/3‑100 (MQ=255)
                                                                                                                                        tccgcCAGGTCGCCCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGACCGAGTCGGTGCTGGAGAACGGCAGCGACGCAGTGTTCGCCGCGTTGTTCt             >  2:491835/3‑150 (MQ=255)
                                                                                                                                        tccgcCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGACCGAGTCGGTGCTGGAGAACGGCAGCGACGCAGTGTTCGCCGCGTTGTTCt             >  1:1007161/3‑150 (MQ=255)
                                                                                                                                                ggCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGACCGAGTCGGTGCTGGAGAACGGCAGCGACGCAGTGTTCGCCGCGTTGTTCTGGTTtgtg     <  1:1139357/150‑1 (MQ=255)
                                                                                                                                                   gACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGACCGAGTCGGTGCTGGAGAACGGCAGCGACGCAGTGTTCGCCGCGTTGTTCTGGTTTgtggtg  <  1:899239/150‑1 (MQ=255)
                                                                                                                                                    aCCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGAcc                                                              >  1:785633/1‑89 (MQ=255)
                                                                                                                                                    aCCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGAcc                                                              <  2:785633/89‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTGGTGCCGCTGACCCTGGTGGCGCTGGTCCTTTCGTGGCTGCCGTATATCGGCTGGCTGGTTGATGTGCTGGCGCTTTACTGCGCCGTGGGGCTGCGCAGCCTGGGCGAGCATGTGCTGCCGGTGGCCAATGCCTTGCGCCAGGGCGACCTGGCGGAGGCGCGGCGCCGGGTTGGCTACCTGGTCAGCCGCGAAACGCGCGAGCTGGACGAGCCCGCGGTGGCCCGAGCGGCGACCGAGTCGGTGCTGGAGAACGGCAGCGACGCAGTGTTCGCCGCGTTGTTCTGGTTTGTGGTG  >  NC_002947/1870042‑1870338

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