Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,253,815 G→A 6.3% G45D (GGT→GAT)  leuB → 3‑isopropylmalate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,253,8150GA 6.3% 192.9 / 2.1 63G45D (GGT→GAT) leuB3‑isopropylmalate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (26/33);  new base A (2/2);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

GGATAGAGGAAAGCATGAGCAAGCAGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGAA  >  NC_002947/2253668‑2253962
                                                                                                                                                    |                                                                                                                                                   
ggATAGAGGAAAGCATGAGCAAGCAGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGt                                                                                                                                                    >  8:215645/1‑149 (MQ=255)
  aTAGAGGAAAGCATGAGCAAGCAGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGATAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTgc                                                                                                                                                  <  2:252989/149‑1 (MQ=255)
  aTAGAGGAAAGCATGAGCAAGCAGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTgc                                                                                                                                                  >  6:94379/1‑149 (MQ=255)
   tAGAGGAAAGCATGAGCAAGCAGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTgcc                                                                                                                                                 <  7:338822/149‑1 (MQ=255)
   tAGAGGAAAGCATGAGCAAGCAGAATCTGATTCTCCCAGCTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACACCGTGCTTGTCGATggc                                                                                                                                                 >  6:132154/1‑147 (MQ=255)
        gAAAGCATGAGCAAGCAGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCAt                                                                                                                                            <  6:125823/149‑1 (MQ=255)
          aaGCATGAGCAAGCAGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCg                                                                                                                                          <  1:59017/149‑1 (MQ=255)
          aaGCATGAGCAAGCAGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCg                                                                                                                                          <  5:128415/149‑1 (MQ=255)
           aGCATGAGCAAGCATATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAATGTACTGGAGCTGGCCAACAGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCg                                                                                                                                          <  1:108996/149‑1 (MQ=255)
               tGAGCAAGCAGATTCTGATTCTCCCAGGGGACGGTATCGGCCCGGAACTCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAg                                                                                                                                     <  4:190824/149‑1 (MQ=255)
                      gcaGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGAACGgtg                                                                                                                              >  5:182765/1‑149 (MQ=255)
                       caGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGgtgt                                                                                                                             <  5:290967/149‑1 (MQ=255)
                       caGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACACCGTGCTTGTCGATGGCGGCACCGCCAATCACgtcgt                                                                                                                             <  7:98848/149‑5 (MQ=37)
                        agcttctgattctCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCCGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGCCCAGGCCGGtgtg                                                                                                                            >  4:96356/4‑149 (MQ=255)
                        aGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGtgtg                                                                                                                            <  5:94379/149‑1 (MQ=255)
                         gattctgattctCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGc                                                                                                                           <  8:257637/149‑1 (MQ=255)
                          attctgattctCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACACCGTGCTTGTCGATGGCGGCACCGCCAATCACGTCGTGcg                                                                                                                          >  2:226196/1‑148 (MQ=25)
                                attctCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGc                                                                                                                    <  3:96211/149‑1 (MQ=255)
                                attctCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGc                                                                                                                    >  8:222907/1‑149 (MQ=255)
                                attctCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACACCGTGCTTGTCGATGGCGGCACCGCCAATCACGTCGTGcgcca                                                                                                                       <  5:132154/146‑5 (MQ=25)
                                 ttctCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGcc                                                                                                                   <  6:180132/149‑1 (MQ=255)
                                    acccAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTACCGCCATCGACAAGCACGGTGTGCCGCTGGCcga                                                                                                                 <  1:303901/147‑1 (MQ=255)
                                             aCGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCt                                                                                                        <  6:92929/148‑1 (MQ=255)
                                             aCGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTg                                                                                                       >  7:331556/1‑149 (MQ=255)
                                             aCGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTg                                                                                                       <  7:215645/149‑1 (MQ=255)
                                                   tCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGt                                                                                                 >  3:196291/1‑149 (MQ=255)
                                                     ggCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTgc                                                                                               >  2:355378/1‑149 (MQ=255)
                                                          ggAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTgcgcgc                                                                                           >  3:11940/1‑148 (MQ=255)
                                                              aTCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGc                                                                                      <  3:96356/149‑1 (MQ=255)
                                                                  tGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGAt                                                                                  >  5:60280/1‑149 (MQ=255)
                                                                     ccGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATg                                                                                 >  6:183037/1‑147 (MQ=255)
                                                                        aGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGtg                                                                            <  4:196291/149‑1 (MQ=255)
                                                                              tCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGctgctg                                                                      >  4:79662/1‑149 (MQ=255)
                                                                                aaGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTggg                                                                    >  6:148187/1‑149 (MQ=255)
                                                                                   gTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGGGCTGCTGGGGGc                                                                 >  8:179489/1‑149 (MQ=255)
                                                                                          aGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTgggc                                                          <  3:79662/149‑1 (MQ=255)
                                                                                           gCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTgggcg                                                         <  5:161345/149‑1 (MQ=255)
                                                                                            cTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTgggcgg                                                        <  5:148187/149‑1 (MQ=255)
                                                                                                ccAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGcc                                                     >  4:80563/1‑148 (MQ=255)
                                                                                                  aaCGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCg                                                    >  6:247482/1‑147 (MQ=255)
                                                                                                  aaCGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGaa                                                  <  1:355378/149‑1 (MQ=255)
                                                                                                    cGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGt                                                <  8:331556/149‑1 (MQ=255)
                                                                                                       aCAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTgg                                              <  3:80563/149‑1 (MQ=255)
                                                                                                       aCAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTgg                                              <  5:247482/149‑1 (MQ=255)
                                                                                                                            gCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTgg                                              >  6:187678/1‑128 (MQ=255)
                                                                                                                            gCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTgg                                              <  5:187678/128‑1 (MQ=255)
                                                                                                                            gCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGCGACAtt                         >  5:177924/1‑149 (MQ=255)
                                                                                                                            gCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACAtt                         >  4:376494/1‑149 (MQ=255)
                                                                                                                             ccTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATcgag                                  <  8:329521/139‑1 (MQ=255)
                                                                                                                             ccTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATcgag                                  >  7:329521/1‑139 (MQ=255)
                                                                                                                               tGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGc                      >  6:45506/1‑149 (MQ=255)
                                                                                                                                 gcgcACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGcc                     <  3:376494/148‑1 (MQ=255)
                                                                                                                                  cgcACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCg                   >  5:127199/1‑149 (MQ=255)
                                                                                                                                     acgacgTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGcc                                                                                                                   <  7:128714/50‑1 (MQ=255)
                                                                                                                                     acgacgTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGcc                                                                                                                   >  8:128714/1‑50 (MQ=255)
                                                                                                                                       gacgTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAgcg              <  6:60280/149‑1 (MQ=255)
                                                                                                                                       gacgTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAgcg              >  6:151310/1‑149 (MQ=255)
                                                                                                                                        acgTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATcgagcg                                >  5:107857/1‑130 (MQ=255)
                                                                                                                                        acgTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATcgagcg                                <  6:107857/130‑1 (MQ=255)
                                                                                                                                              ttGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCctg       <  6:182765/149‑1 (MQ=255)
                                                                                                                                                ggcggcACCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCctgc      <  6:155065/148‑1 (MQ=255)
                                                                                                                                                 gCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCctgct     >  6:84011/1‑148 (MQ=255)
                                                                                                                                                  cGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGa   <  5:45506/149‑1 (MQ=255)
                                                                                                                                                  cGGTGCCGCCATCGACAAGCACGGTGTGCCGCGGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCGGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGa   <  7:151731/149‑1 (MQ=255)
                                                                                                                                                   ggTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCg                       >  4:6865/1‑128 (MQ=255)
                                                                                                                                                   ggTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCg                       <  3:6865/128‑1 (MQ=255)
                                                                                                                                                   ggTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGaa  <  7:222907/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGATAGAGGAAAGCATGAGCAAGCAGATTCTGATTCTCCCAGGTGACGGTATCGGCCCGGAAATCATGGCCGAGGCAGTCAAGGTACTGGAGCTGGCCAACGACAAGTTCCAGCTCGGCTTCAGCCTGGCGCACGACGTGATTGGCGGTGCCGCCATCGACAAGCACGGTGTGCCGCTGGCCGACGAAACCCTGGAGCGTGCGCGCAAGGCCGATGCAGTGCTGCTGGGGGCTGTGGGCGGGCCGAAGTGGGACAAGATCGAGCGAGACATTCGCCCGGAGCGCGGCCTGCTGAA  >  NC_002947/2253668‑2253962

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: