Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,797,501 A→T 6.2% 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,5010AT 6.2% 228.6 / 2.7 75C165S (TGC→AGC) vdhvanillin dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (37/33);  new base T (2/3);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

GCTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCTG  >  NC_002947/3797354‑3797645
                                                                                                                                                   |                                                                                                                                                
gCTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAg                                                                                                                                                 >  7:78297/1‑149 (MQ=255)
 cTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGCGCACCACGGTGTTGCCGCAgg                                                                                                                                                >  4:141873/1‑149 (MQ=255)
  tCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGc                                                                                                                                               >  6:417933/1‑149 (MQ=255)
    gacgaTGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGc                                                                                                                                               <  8:204625/147‑1 (MQ=255)
      cgaTGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGc                                                                                                                                           <  2:72581/149‑1 (MQ=255)
       gaTGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCTGCg                                                                                                                                          >  5:393769/1‑149 (MQ=255)
         tGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCgg                                                                                                                                         <  4:263201/148‑1 (MQ=255)
           gCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCTGGCCTGCGGCaa                                                                                                                                      <  6:393769/149‑2 (MQ=255)
            cGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTTGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCTGGCCTGCGGCaa                                                                                                                                      <  4:260036/148‑2 (MQ=255)
             gggggCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCTGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCatgg                                                                                                                                    <  1:380305/149‑1 (MQ=255)
                 gCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCatggcaat                                                                                                                                <  2:77910/149‑1 (MQ=255)
                       tGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCatggcaatgg                                                                                                                              <  6:75784/145‑1 (MQ=255)
                       tGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCatggcaatgg                                                                                                                              >  5:75784/1‑145 (MQ=255)
                         cggcgCATTGCTGAGGACAATGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGAGCAGCCGCTGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGCTGCCGCAGGCCAGCGGCATGGCAATGGCAcgc                                                                                                                         <  5:342128/148‑1 (MQ=255)
                         ccgcgcATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCAcgc                                                                                                                         <  5:348003/146‑1 (MQ=255)
                           gcgcATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGcc                                                                                                                                              <  2:190633/125‑1 (MQ=255)
                           gcgcATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGcc                                                                                                                                              >  1:190633/1‑125 (MQ=255)
                           gcgcATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTg                                                                                                                      >  5:68509/1‑149 (MQ=255)
                               aTTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGc                                                                                                                    >  6:415944/1‑147 (MQ=255)
                                       tgaCATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCAcc                                                                                                          >  3:273187/1‑149 (MQ=255)
                                            ttGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGc                                                                                                     <  2:188174/149‑1 (MQ=255)
                                            ttGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCGTCAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCCTGGCCAGTATCACCGGGGc                                                                                                     >  5:75372/1‑149 (MQ=255)
                                               aCCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGtt                                                                                                  <  5:415944/149‑1 (MQ=255)
                                                  aCGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCa                                                                                               <  8:50707/149‑1 (MQ=255)
                                                    gCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACg                                                                                             <  7:20957/149‑1 (MQ=255)
                                                       gTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTg                                                                                          <  6:68509/149‑1 (MQ=255)
                                                       gTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTg                                                                                          >  1:313287/1‑149 (MQ=255)
                                                       gTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAACCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTg                                                                                          <  3:39768/149‑1 (MQ=255)
                                                           ccGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGAt                                                                                      <  5:358405/149‑1 (MQ=255)
                                                               tGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGccc                                                                                  <  4:178487/149‑1 (MQ=255)
                                                                  cTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAac                                                                               >  4:60692/1‑149 (MQ=255)
                                                                   tGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAaca                                                                              >  7:219559/1‑149 (MQ=255)
                                                                          tGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCAcgcc                                                                       >  3:259892/1‑149 (MQ=255)
                                                                             agcaCCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCa                                                                    <  5:95205/149‑1 (MQ=255)
                                                                               caCCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAgg                                                                  <  6:121390/149‑1 (MQ=255)
                                                                                aCCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAggg                                                                 <  4:272525/149‑1 (MQ=255)
                                                                                      ccGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGgcgcgcg                                                           >  3:374672/1‑149 (MQ=255)
                                                                                       cGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGgcgcgcg                                                           >  6:255537/1‑148 (MQ=255)
                                                                                       cGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGgcgcgcg                                                           >  7:222235/1‑148 (MQ=255)
                                                                                         aTCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCAGCAGGGCGCGCGCAg                                                        >  6:147638/1‑149 (MQ=255)
                                                                                            aGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACCCCACGCCGCAGGGCGCGCGCAGGGc                                                     >  6:77845/1‑149 (MQ=255)
                                                                                              ccGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCa                                                   >  7:57447/1‑149 (MQ=255)
                                                                                                    gACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAAGGGCCCGCGGGCCCGGAAAAACCGGGGCATACCCCGGGCCGATGCCCCCCCCCCCCCCGCCCGGCGCCCGCACGGCCAtgccga                                             >  8:151394/1‑144 (MQ=255)
                                                                                                       cgccgGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGc                                          >  3:259354/1‑149 (MQ=255)
                                                                                                       cgccgGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGc                                          >  8:174624/1‑149 (MQ=255)
                                                                                                               tCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCggg                                   >  2:198481/1‑148 (MQ=255)
                                                                                                                cAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGaa                                 >  3:119084/1‑149 (MQ=255)
                                                                                                                   cTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGt                              >  4:113729/1‑149 (MQ=255)
                                                                                                                      cGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAg                                                        <  1:204281/120‑1 (MQ=255)
                                                                                                                      cGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAg                                                        >  2:204281/1‑120 (MQ=255)
                                                                                                                      cGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTcc                            >  2:327447/1‑148 (MQ=255)
                                                                                                                        aGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCg                                                                                        <  8:280258/86‑1 (MQ=255)
                                                                                                                        aGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCg                                                                                        >  7:280258/1‑86 (MQ=255)
                                                                                                                        aGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAg                         <  5:417933/149‑1 (MQ=255)
                                                                                                                         ggCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAgg                        <  3:113729/149‑1 (MQ=255)
                                                                                                                          gCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAggg                       <  8:219559/149‑1 (MQ=255)
                                                                                                                          gCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAggg                       <  8:57447/149‑1 (MQ=255)
                                                                                                                            cTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGa                     >  5:118345/1‑149 (MQ=255)
                                                                                                                             ttGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGAt                    >  4:426240/1‑149 (MQ=255)
                                                                                                                               gAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGAt                    <  4:119084/147‑1 (MQ=255)
                                                                                                                                  caccacGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACtt               >  5:32671/1‑149 (MQ=255)
                                                                                                                                    ccacGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCa             >  7:316309/1‑149 (MQ=255)
                                                                                                                                     cacGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCAc            >  2:130175/1‑149 (MQ=255)
                                                                                                                                     cacGGTGCTGCCGCAGTCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCAc            <  6:75372/149‑1 (MQ=255)
                                                                                                                                       cGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCg          <  4:273187/149‑1 (MQ=255)
                                                                                                                                       cGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCg          <  8:222235/149‑1 (MQ=255)
                                                                                                                                           catGCCGCTGGCCTGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGgcgcgc                                                            >  7:111224/3‑95 (MQ=255)
                                                                                                                                           catGCCGCTGGCCTGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGgcgcgc                                                            <  8:111224/93‑1 (MQ=255)
                                                                                                                                            atGCCGCTGGCCTGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGAt     >  1:289604/2‑149 (MQ=255)
                                                                                                                                               ccGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAggg                       <  7:52524/128‑1 (MQ=255)
                                                                                                                                               ccGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCt   <  6:32671/148‑1 (MQ=255)
                                                                                                                                               ccGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCTg  >  5:89995/1‑149 (MQ=255)
                                                                                                                                               acGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAggg                       >  8:52524/2‑128 (MQ=255)
                                                                                                                                                 gCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCa             <  1:51266/136‑1 (MQ=255)
                                                                                                                                                 gCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCa             >  2:51266/1‑136 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCTCGACGATGGCGGGGGCATCCTGCGGCGCATTGCTGATGACATTGACCACGCCGTCGCCGATGCCTGCATCGTGCAGCACCTGGCCGATCAGCCGATGGACCGCCGGGCTCAGCTCCGAGGCCTTGAGCACCACGGTGTTGCCGCAGGCCAGCGGCATGGCAATGGCACGCGTGGCCAGTATCACCGGGGCGTTCCACGGTGCGATGCCCAACACCACGCCGCAGGGCGCGCGCAGGGCCATTGCGAAGCTGCCGGGAACGTCCGAGGGGATCACTTCACCGGTGATCTG  >  NC_002947/3797354‑3797645

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