Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,151,942 Δ1 bp 8.5% intergenic (+132/+19) PP_3655 → / ← PP_3656 cytosine/purine/uracil/thiamine/allantoin permease family protein/aromatic compound‑specific porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,151,9420G. 8.5% 179.8 / 10.8 59intergenic (+132/+19)PP_3655/PP_3656cytosine/purine/uracil/thiamine/allantoin permease family protein/aromatic compound‑specific porin
Reads supporting (aligned to +/- strand):  ref base G (26/28);  new base . (2/3);  total (28/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 = 8.93e-01

GCGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCACTCTGAATAACATAAG  >  NC_002947/4151799‑4152088
                                                                                                                                                 |                                                                                                                                                   
gCGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAggccgg                                                                                                                                                <  2:201280/149‑1 (MQ=255)
 cGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAggccggc                                                                                                                                               >  7:272842/1‑149 (MQ=255)
 cGATGGCTTAATCGGTACTGGCCTTTCGTTACAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAggccggc                                                                                                                                               <  3:315253/149‑1 (MQ=255)
  gATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAggccggcc                                                                                                                                              <  5:251256/149‑1 (MQ=255)
  gATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAgg                                                                                                                                                    >  2:295226/1‑143 (MQ=255)
  gATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAgg                                                                                                                                                    <  1:295226/143‑1 (MQ=255)
        ttAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTc                                                                                                                                        >  1:321622/1‑149 (MQ=255)
          aaTCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTc                                                                                                                                        <  1:308840/147‑1 (MQ=255)
          aaTCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCg                                                                                                                                       <  3:63983/148‑1 (MQ=255)
          aaTCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATAGCG‑‑AAACGGGCCGGCCTTCTGAAGCCGGCCCAtttttc                                                                                                                                      <  6:216678/147‑4 (MQ=37)
                     gCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGATTAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGt                                                                                                                           >  1:296203/1‑149 (MQ=255)
                     gCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGaa                                                                                                                             <  7:36752/147‑1 (MQ=255)
                     gCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGaa                                                                                                                             <  8:305794/147‑1 (MQ=255)
                      ccTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGaa                                                                                                                             <  4:80984/146‑1 (MQ=255)
                      ccTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGaa                                                                                                                             >  3:80984/1‑146 (MQ=255)
                          tctttGCAAAATATGAAACGGCTGGGTGCAGGCGTCCGACCGCCTGGTGCAATTCAGGCGGAAAACATCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGaa                                                                                                                       <  8:272842/145‑1 (MQ=255)
                           cGTTGCAAAATATGAAACGGCTGGGGGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAAACAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGaa                                                                                                                       <  1:256623/148‑1 (MQ=255)
                           cGTTGCAAAAGATGAAACGTCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAg                                                                                                                      <  8:268136/148‑1 (MQ=255)
                               gCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATTGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGgtg                                                                                                                   >  6:172554/1‑147 (MQ=255)
                                cAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAg                                                                                                                >  8:294861/1‑149 (MQ=255)
                                 aaaaGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCCTTCTGAAGCCGGccca                                                                                                                                            >  8:270581/1‑119 (MQ=255)
                                 aaaaGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCCTTCTGAAGCCGGccca                                                                                                                                            <  7:270581/120‑2 (MQ=255)
                                 aaaaGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAgc                                                                                                               >  2:267946/1‑149 (MQ=255)
                                 aaaaGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAgc                                                                                                               >  8:55310/1‑149 (MQ=255)
                                 aaaaGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAgc                                                                                                               >  5:37242/1‑149 (MQ=255)
                                 aaaaGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGACGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAgc                                                                                                               >  5:32497/1‑149 (MQ=255)
                                     gATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCAt                                                                                                           <  6:148796/149‑1 (MQ=255)
                                        gAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGAt                                                                                                        >  2:244841/1‑149 (MQ=255)
                                         aaaCGGCTTGGTGCAGGCGGCCGTACTCCTGGTGCAATTCTGGCTGAAAAAAGCGTACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTGAGAAGGCCTGCCCGTTTCGCTATTTAGAAGTTGAAGTTGTAGCGCATGATc                                                                                                       <  2:212609/149‑1 (MQ=255)
                                         aaaCGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCTGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATc                                                                                                       <  1:244841/149‑1 (MQ=255)
                                                    accAGGCGTTCGGACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGAtttt                                                                                             >  8:279648/3‑148 (MQ=255)
                                                         gCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGAtttt                                                                                             >  3:180016/1‑143 (MQ=255)
                                                         gCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGAtttt                                                                                             <  4:180016/143‑1 (MQ=255)
                                                            tCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAgcgc                                                                                    >  4:217997/1‑149 (MQ=255)
                                                               gAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTg                                                                                  >  1:287966/1‑148 (MQ=255)
                                                                   gCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGc                                                                             >  3:73332/1‑149 (MQ=255)
                                                                   gCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGc                                                                             >  2:44542/1‑149 (MQ=255)
                                                                     cTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCg                                                                            >  6:118126/1‑148 (MQ=255)
                                                                       ggTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAg                                                                        >  2:48460/1‑149 (MQ=255)
                                                                        gTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAg                                                                        <  4:73332/148‑1 (MQ=255)
                                                                           cAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTc                                                                     <  3:217997/148‑1 (MQ=255)
                                                                              ttCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGc                                                                 <  7:279648/149‑1 (MQ=255)
                                                                              ttCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGc                                                                 <  7:55310/149‑1 (MQ=255)
                                                                                 aGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGt                                                              <  2:287966/149‑1 (MQ=255)
                                                                                 aGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGt                                                              <  3:136632/149‑1 (MQ=255)
                                                                                 aGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGt                                                              >  1:39073/1‑149 (MQ=255)
                                                                                   gCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTc                                                            >  3:69170/1‑149 (MQ=255)
                                                                                      gAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTg                                                         <  6:37242/149‑1 (MQ=255)
                                                                                       aaaaCAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTg                                                         <  6:32497/148‑1 (MQ=255)
                                                                                       aaaaCAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTg                                                         >  4:293460/1‑148 (MQ=255)
                                                                                                 aCCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGcc                                               <  5:172554/148‑1 (MQ=255)
                                                                                                   ccTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACt                                            <  1:48460/149‑1 (MQ=255)
                                                                                                        gATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTTAGGTCGAGCGGTTCGTGGAGTTACGCCACTgaaa                                        <  1:44542/149‑1 (MQ=255)
                                                                                                              gAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAaca                                 >  2:81580/1‑149 (MQ=255)
                                                                                                                 aTAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACAcc                               >  2:306282/1‑149 (MQ=255)
                                                                                                                             aaaCGGGCCGGCCTTCTGAAGCCGGCCCATTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCg                       <  5:255743/145‑1 (MQ=255)
                                                                                                                             aaaCGGGCCGGCCTTCTGAAGCCGGCCCATTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCg                       >  6:255743/1‑145 (MQ=255)
                                                                                                                                  ggCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCActct              <  1:306282/149‑1 (MQ=255)
                                                                                                                                   gCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCActct              >  3:16980/1‑148 (MQ=255)
                                                                                                                                       gCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCACTCTGAAt          >  4:316354/1‑148 (MQ=255)
                                                                                                                                         cTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATCATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCACTCTGAAt          <  1:267946/147‑1 (MQ=255)
                                                                                                                                            cAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCTTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCACTCTGAATAAc       <  1:81580/147‑1 (MQ=255)
                                                                                                                                              gAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCACTCTGAATAACATaa   <  3:293460/149‑1 (MQ=255)
                                                                                                                                               aaGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCACTCTGAATAACAt     >  8:174191/1‑146 (MQ=255)
                                                                                                                                                aGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCACTCTGAATAACATAag  <  4:69170/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCACTCTGAATAACATAAG  >  NC_002947/4151799‑4152088

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