Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,151,938 A→T 9.4% intergenic (+128/+23) 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,9380AT 9.4% 179.8 / 6.1 63intergenic (+128/+23)PP_3655/PP_3656cytosine/purine/uracil/thiamine/allantoin permease family protein/aromatic compound‑specific porin
Reads supporting (aligned to +/- strand):  ref base A (27/30);  new base T (2/4);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CTGGCAGCGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCACTCTGAATAAC  >  NC_002947/4151793‑4152083
                                                                                                                                                   |                                                                                                                                                  
cTGGCAGCGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGCAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGaa                                                                                                                                                 <  7:257041/149‑1 (MQ=255)
      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                                                                                                          >  5:37242/1‑149 (MQ=255)
                                       aaaaGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAgc                                                                                                          >  2:267946/1‑149 (MQ=255)
                                       aaaaGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAgc                                                                                                          >  8:55310/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                                                         >  1:39073/1‑149 (MQ=255)
                                                                                       aGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGt                                                         <  3:136632/149‑1 (MQ=255)
                                                                                       aGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGt                                                         <  2:287966/149‑1 (MQ=255)
                                                                                         gCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTc                                                       >  3:69170/1‑149 (MQ=255)
                                                                                            gAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTg                                                    <  6:37242/149‑1 (MQ=255)
                                                                                             aaaaCAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTg                                                    >  4:293460/1‑148 (MQ=255)
                                                                                             aaaaCAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTg                                                    <  6:32497/148‑1 (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                  >  6:255743/1‑145 (MQ=255)
                                                                                                                                   aaaCGGGCCGGCCTTCTGAAGCCGGCCCATTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCg                  <  5:255743/145‑1 (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)
                                                                                                                                                   |                                                                                                                                                  
CTGGCAGCGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCACTCTGAATAAC  >  NC_002947/4151793‑4152083

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: