Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,797,544 A→T 5.8% C165S (TGC→AGC)  vdh ← vanillin dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,797,5440AT 5.8% 365.8 / 10.3 122C165S (TGC→AGC) vdhvanillin dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (62/51);  new base T (4/3);  total (67/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

Only 100 of 124 total aligned reads displayed.

CCGCTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCTGCGTGG  >  NC_002947/3797395‑3797693
                                                                                                                                                      |                                                                                                                                                     
ccGCTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCa                                                                                                                                                       <  1:1858571/150‑1 (MQ=255)
 cGCTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAg                                                                                                                                                      <  1:2207020/150‑1 (MQ=255)
    tCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGcc                                                                                                                                                   >  1:1572411/1‑150 (MQ=255)
      gacgaTGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAg                                                                                                                                                 >  2:2248133/1‑150 (MQ=255)
       acgaTGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCACACGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGc                                                                                                                                                >  1:25072/1‑150 (MQ=255)
          aTGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTCAGCACCACGGTGTTGCCGCAGGCCAGCGGc                                                                                                                                             <  2:1410651/150‑1 (MQ=255)
           tGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCa                                                                                                                                            <  1:1771883/150‑1 (MQ=255)
           tGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGATCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCa                                                                                                                                            <  2:1902564/150‑1 (MQ=255)
            ggCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCTGGCCTGCGGCaa                                                                                                                                           <  1:66612/150‑2 (MQ=255)
            ggCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCat                                                                                                                                           <  2:35369/150‑1 (MQ=255)
             gCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCatg                                                                                                                                          <  2:1866927/150‑1 (MQ=255)
              cGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCatgg                                                                                                                                         <  1:1035523/150‑1 (MQ=255)
              cGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCatgg                                                                                                                                         <  2:1803408/150‑1 (MQ=255)
                     aTCCTGCGGCGCATTGCTGATGACATTGACCACACCGACCACGTTAACCGCGTCGTGCAGCACCTGGCCTATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCatggcaatggc                                                                                                                                  >  1:46359/1‑150 (MQ=255)
                      tCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCatggcaatggcc                                                                                                                                 >  1:933118/1‑149 (MQ=255)
                      tCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCatggcaatggca                                                                                                                                 >  2:1845408/1‑150 (MQ=255)
                        cTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCAcg                                                                                                                               >  2:2047664/1‑150 (MQ=255)
                             gcgcATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTgg                                                                                                                          >  1:1281450/1‑150 (MQ=255)
                                cATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCa                                                                                                                       >  1:2021280/1‑150 (MQ=255)
                                  ttGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGt                                                                                                                     >  1:304128/1‑150 (MQ=255)
                                           aCATTGACCACGCCGTCGCCGATGTCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCggg                                                                                                            <  2:282643/150‑1 (MQ=255)
                                           aCATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGAGCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCATGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCggg                                                                                                            <  1:2297142/150‑1 (MQ=255)
                                             aTTNACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTTGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGc                                                                                                          >  2:1498117/1‑150 (MQ=255)
                                             aTTGACCACGCCGTCGCCGATGCCTGCATCGTTCAGCACCTTGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCAGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGc                                                                                                          >  2:76500/1‑150 (MQ=255)
                                              ttGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCTGGCCTGCGGCaacac                                                                                                                                        <  2:977973/119‑5 (MQ=255)
                                              ttGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCTGGCCTGCGGCaacac                                                                                                                                        >  1:976506/1‑115 (MQ=255)
                                              ttGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCTGGCCTGCGGCaacac                                                                                                                                        >  1:977973/1‑115 (MQ=255)
                                              ttGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCTGGCCTGCGGCaacac                                                                                                                                        <  2:976506/119‑5 (MQ=255)
                                               tGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGAACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGt                                                                                                        <  1:1689579/150‑1 (MQ=255)
                                                gACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCgtggaacg                                                                                                       >  1:562805/1‑146 (MQ=255)
                                                   cACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCa                                                                                                    >  1:1851451/1‑150 (MQ=255)
                                                    aCGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCCGCGGAAGTCCTGTGGCCCGCCCGGCCAGTATCACCGGGGCGTTCCAc                                                                                                   >  1:625224/1‑150 (MQ=255)
                                                    aCGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAg                                                                                                                      >  1:593530/1‑131 (MQ=255)
                                                    aCGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAg                                                                                                                      <  2:593530/131‑1 (MQ=255)
                                                    aCGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCAACACGGTGTTGCGAGCAGCGTCGATCTGGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCAc                                                                                                   >  1:1742389/1‑150 (MQ=255)
                                                         gTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGTAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGATCCACGGTGc                                                                                              <  2:92114/150‑1 (MQ=255)
                                                          tCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCg                                                                                             >  1:872406/1‑150 (MQ=255)
                                                            gCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCTGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGAt                                                                                           <  1:1227714/150‑1 (MQ=255)
                                                               gATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGcc                                                                                        >  1:448545/1‑150 (MQ=255)
                                                                aTGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGccc                                                                                       >  2:1564612/1‑150 (MQ=255)
                                                                aTGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGccc                                                                                       <  1:1845408/150‑1 (MQ=255)
                                                                  gCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCaaa                                                                                     <  2:1807119/150‑1 (MQ=255)
                                                                  gCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCaa                                                                                     <  2:534193/150‑1 (MQ=255)
                                                                    cTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAaca                                                                                   >  2:171830/1‑150 (MQ=255)
                                                                    cTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAaca                                                                                   >  2:2270475/1‑150 (MQ=255)
                                                                     tGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAacac                                                                                  >  2:2189057/1‑150 (MQ=255)
                                                                     tGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAacac                                                                                  >  1:2391009/1‑150 (MQ=255)
                                                                          cGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCCGGCCCGCGGCATGGCAATGGCACGCGTAGCCAGTATCACCTGAGCGTTCCACTGGGCGATGGCCAACACCAcgc                                                                             <  2:25072/150‑1 (MQ=255)
                                                                             gcagcaCCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCg                                                                                                                                               <  2:332148/81‑1 (MQ=255)
                                                                             gcagcaCCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCg                                                                                                                                               >  1:332148/1‑81 (MQ=255)
                                                                                gcaCCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGAGTGTGCAGTAGCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAgg                                                                       >  2:666013/1‑150 (MQ=255)
                                                                                 caCCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAggg                                                                      >  1:2283361/1‑150 (MQ=255)
                                                                                 caCCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAggg                                                                      >  1:2291952/1‑150 (MQ=255)
                                                                                     tGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGgcgcg                                                                  >  1:985672/1‑150 (MQ=255)
                                                                                       gCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCg                                                                                                                                               <  2:1129415/71‑1 (MQ=255)
                                                                                       gCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCAc                                                                                                   >  2:2154758/1‑115 (MQ=255)
                                                                                       gCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCAc                                                                                                   <  1:2154758/115‑1 (MQ=255)
                                                                                          gATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAg                                                             <  2:2095253/150‑1 (MQ=255)
                                                                                            tCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAggg                                                           <  2:11029/150‑1 (MQ=255)
                                                                                                 cGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCCTGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCAtt                                                      <  2:448545/150‑1 (MQ=255)
                                                                                                 cGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCgg                                                                                                             <  2:2298599/95‑1 (MQ=255)
                                                                                                 cGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCgg                                                                                                             >  1:2298599/1‑95 (MQ=255)
                                                                                                   aTGGACCGCCGGGCTCACCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTa                                                                                                                    <  2:1377321/86‑1 (MQ=255)
                                                                                                       aCCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCa                                                                         <  1:445014/125‑1 (MQ=255)
                                                                                                         cgccgGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCt                                              >  2:1399365/1‑150 (MQ=255)
                                                                                                          gccgGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGGGTGGCCAGTATCAGCGGGGCGGGCCCCGGTGCGATGCGCAACACCACGCCGCAGCGCGCGCGGAGGTCCATGGCgcacatg                                             >  2:501110/1‑144 (MQ=255)
                                                                                                              ggCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCgg                                         <  2:879027/150‑1 (MQ=255)
                                                                                                                 tCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCCGGACTGCGGCATGGCAACGTCACGCGTGGCCAGTAACACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGaa                                      >  1:1003956/1‑150 (MQ=255)
                                                                                                                 tCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGaa                                      <  2:1851451/150‑1 (MQ=255)
                                                                                                                 tCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGCCGCCAGTATCACCGGGGCGTTCCACGGTGCGACGCCCAACACCACGCCGCAGGGCGAGCGCAGGGGCATTGCCAAGCTGCCGGGaa                                      >  2:610088/1‑150 (MQ=255)
                                                                                                                  cAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCACCACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAAc                                     <  1:1992472/150‑1 (MQ=255)
                                                                                                                   aGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGTACg                                    >  2:37213/1‑150 (MQ=255)
                                                                                                                     cTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTc                                  >  2:953334/1‑150 (MQ=255)
                                                                                                                        cGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGaa                                                 <  1:1078880/132‑1 (MQ=255)
                                                                                                                        cGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCg                                >  1:843354/1‑149 (MQ=255)
                                                                                                                        cGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCg                                <  2:843354/149‑1 (MQ=255)
                                                                                                                         gAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAg                              >  1:243575/1‑150 (MQ=255)
                                                                                                                         gAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAg                              <  1:1564612/150‑1 (MQ=255)
                                                                                                                          aGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAgg                             <  1:1006806/150‑1 (MQ=255)
                                                                                                                            gCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACAAAGCCGATGGGCGCGCGCGTGGCTCTGGGTATGCTTCCGGGAACGTCCGAgggg                           >  2:1204613/1‑150 (MQ=255)
                                                                                                                             ccTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGa                          >  1:849856/1‑150 (MQ=255)
                                                                                                                             ccTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACAGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGa                          <  2:380124/150‑1 (MQ=255)
                                                                                                                                       cacGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCaa                                                                                     >  1:150032/1‑81 (MQ=255)
                                                                                                                                       cacGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCaa                                                                                     <  2:150032/81‑1 (MQ=255)
                                                                                                                                       cacGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAgg                                                                       <  2:1014524/95‑1 (MQ=255)
                                                                                                                                       cacGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAgg                                                                       >  1:1014524/1‑95 (MQ=255)
                                                                                                                                       cacGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCAcc                >  2:1895033/1‑150 (MQ=255)
                                                                                                                                       cacGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCAcc                >  2:1175917/1‑150 (MQ=255)
                                                                                                                                       cacGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGAGTAACACTTCAcc                >  2:662792/1‑150 (MQ=255)
                                                                                                                                        acGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGAATCTCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCg               <  1:1204613/150‑1 (MQ=255)
                                                                                                                                         cGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCGCCGGGGCGTTCCGCGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGACCGAGCGGAGCACTTCACCgg              >  2:1322222/1‑150 (MQ=255)
                                                                                                                                           gtgtTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTg            <  1:2189057/150‑1 (MQ=255)
                                                                                                                                              ttGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATc         >  1:1765598/1‑150 (MQ=255)
                                                                                                                                              atGCCGCTGGCCTGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATc         >  2:779883/2‑150 (MQ=255)
                                                                                                                                               tGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCt        <  1:1175917/150‑1 (MQ=255)
                                                                                                                                                   cGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCGGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCTGCg     >  2:1279258/1‑150 (MQ=255)
                                                                                                                                                   cGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCAtt                                                      >  2:528265/1‑101 (MQ=255)
                                                                                                                                                   cGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCAtt                                                      <  1:528265/101‑1 (MQ=255)
                                                                                                                                                   cGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCTGCg     >  1:2189264/1‑150 (MQ=255)
                                                                                                                                                    gCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCTGCGt    >  2:1561369/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CCGCTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCTGCGTGG  >  NC_002947/3797395‑3797693

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: