Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,151,986 Δ1 bp 5.7% 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,9860G. 5.7% 224.4 / 6.9 70intergenic (+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 (29/37);  new base . (2/2);  total (31/39)
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.73e-01

GCAGCGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCACTCTGAATAACA  >  NC_002947/4151840‑4152128
                                                                                                                                                   |                                                                                                                                              
gCAGCGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAggcc                                                                                                                                             >  2:1750412/1‑150 (MQ=255)
gCAGCGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAggcc                                                                                                                                             <  1:7067/150‑1 (MQ=255)
   gCGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGACTATTTGGCGAAAAATGGGCCGGCTTCAGAAggccggc                                                                                                                                          <  1:1095737/150‑1 (MQ=255)
   gCGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATACATGGCGAAAAATGGGCCGGCTTCAGAAggccggc                                                                                                                                          <  2:421283/150‑1 (MQ=255)
   gCGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAggccggc                                                                                                                                          <  1:955394/150‑1 (MQ=255)
      aTGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCg                                                                                                                                       >  1:483439/1‑150 (MQ=255)
       tGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCCTTCTGAAGCCGGcccat                                                                                                                                      <  1:1616518/150‑3 (MQ=255)
       tGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGt                                                                                                                                      >  2:1213551/1‑150 (MQ=255)
           ttAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCg                                                                                                                                  >  2:91497/1‑150 (MQ=255)
            tAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGGACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGc                                                                                                                                 >  2:1292263/1‑150 (MQ=255)
                 ggTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTAttt                                                                                                                            <  2:143199/150‑1 (MQ=255)
                     cTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGATTCAGAAGGCCGGCCCGTTTCGCTATTTAGaa                                                                                                                        <  2:1069644/150‑1 (MQ=255)
                       ggCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGt                                                                                                                      <  1:91497/150‑1 (MQ=255)
                        gCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGtt                                                                                                                     >  1:108874/1‑150 (MQ=255)
                         ccTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTg                                                                                                                    >  1:1275972/1‑150 (MQ=255)
                            ttCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAg                                                                                                                 <  1:627570/150‑1 (MQ=255)
                               gTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGgtg                                                                                                              >  1:1437327/1‑150 (MQ=255)
                                  gCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAg                                                                                                           <  1:1213551/150‑1 (MQ=255)
                                   cAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAg                                                                                                                          >  2:48189/1‑134 (MQ=255)
                                   cAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAg                                                                                                                          <  1:48189/134‑1 (MQ=255)
                                    aaaaGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAgcg                                                                                                         >  2:777925/1‑150 (MQ=255)
                                    aaaaGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAgcg                                                                                                         >  2:1616900/1‑150 (MQ=255)
                                      aagaagAAACGGCTGGGTGCAGGCGTCCGTACGCCTGCTGCAATTCAGGCGGAAAACAGCGAACCCTTCTATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCa                                                                                                       <  1:1251788/150‑1 (MQ=255)
                                         aTGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGa                                                                                                    >  2:194994/1‑150 (MQ=255)
                                             aaCGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGTCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAg                                                                                                <  2:1112603/150‑1 (MQ=255)
                                             aaCGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAg                                                                                                >  2:1706871/1‑150 (MQ=255)
                                             aaCGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAg                                                                                                <  2:293482/150‑1 (MQ=255)
                                                  cTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGAt                                                                                           <  1:777925/150‑1 (MQ=255)
                                                  cTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTAGATTTGGCATAAATGGGTGAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGAt                                                                                           <  1:1175804/150‑1 (MQ=255)
                                                    gggTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGAttt                                                                                         <  2:1234659/150‑1 (MQ=255)
                                                     ggTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCAAtttt                                                                                        <  2:24074/150‑1 (MQ=255)
                                                       tGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAg                                                                                                >  1:970931/1‑140 (MQ=255)
                                                       tGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAg                                                                                                <  2:970931/140‑1 (MQ=255)
                                                       accAGGCGTTCGGACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCa                                                                                      >  1:371217/3‑150 (MQ=255)
                                                        gCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCAt                                                                                     <  2:1436012/150‑1 (MQ=255)
                                                          aGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTa                                                                                   >  2:1212049/1‑150 (MQ=255)
                                                           ggCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAggcc                                                                                                                                             >  2:211770/1‑91 (MQ=255)
                                                           ggCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAggcc                                                                                                                                             <  1:211770/91‑1 (MQ=255)
                                                            gCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAgg                                                                                               >  1:1185022/1‑136 (MQ=255)
                                                            gCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAgg                                                                                               <  2:1185022/136‑1 (MQ=255)
                                                                      gCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCg                                                                       <  2:371217/150‑1 (MQ=255)
                                                                      gCCTGGTGCAATTCAGGCGGAAAAAAGCGAACCCTTCGAATGGGAATAAATGGCAAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCg                                                                       <  2:670111/150‑1 (MQ=255)
                                                                       ccTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGt                                                                      >  2:675752/1‑150 (MQ=255)
                                                                        cTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTa                                                                     <  1:1192858/150‑1 (MQ=255)
                                                                        cTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTa                                                                     <  2:1164341/150‑1 (MQ=255)
                                                                            tGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAg                                                                                                                 <  2:429749/102‑1 (MQ=255)
                                                                            tGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAg                                                                                                                 >  1:429749/1‑102 (MQ=255)
                                                                            tGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTc                                                                 <  1:1750412/150‑1 (MQ=255)
                                                                                aTTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGc                                                             <  2:1148564/150‑1 (MQ=255)
                                                                                aTTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGc                                                             >  1:1598237/1‑150 (MQ=255)
                                                                                 ttCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAggccggc                                                                                                                                          <  2:68240/72‑1 (MQ=255)
                                                                                 ttCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAggccggc                                                                                                                                          >  1:68240/1‑72 (MQ=255)
                                                                                    aGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAAAAGGCAGGCTGCAGAAGGCCGGCCCGGGTCGCTATATAGAAGTCGAAGGCGTGGCGCATGATCAGGCGATTTTCATTAGCGGTATCAGCGTAGATCGaaagatc                                                         >  2:926761/1‑144 (MQ=255)
                                                                                      gCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCt                                                                              <  1:106483/127‑1 (MQ=255)
                                                                                      gCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCt                                                                              >  2:106483/1‑127 (MQ=255)
                                                                                        ggAAAACAGCGATCCCTTCGATTGGGAATAACTGGCGAAAAATGGGCCGGCTTCAGATGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGAGCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTg                                                     <  2:669657/150‑1 (MQ=255)
                                                                                          aaaaCAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGa                                                   <  1:276609/150‑1 (MQ=255)
                                                                                          aaaaCAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGa                                                   <  2:1664105/150‑1 (MQ=255)
                                                                                              cAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAgcg                                                                                >  2:333030/1‑117 (MQ=255)
                                                                                              cAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAgcg                                                                                <  1:333030/117‑1 (MQ=255)
                                                                                                   aaCCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTGCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCa                                          <  1:573319/150‑1 (MQ=255)
                                                                                                       cTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTga                                      >  1:1165621/1‑150 (MQ=255)
                                                                                                        ttCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGt                                                                            <  2:1122784/111‑1 (MQ=255)
                                                                                                        ttCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGt                                                                            >  1:1122784/1‑111 (MQ=255)
                                                                                                        ttCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTgaa                                     >  2:884158/1‑150 (MQ=255)
                                                                                                              tGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTgaaacgaaa                               <  1:884158/150‑1 (MQ=255)
                                                                                                                ggAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGGGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAaca                             <  2:629402/150‑1 (MQ=255)
                                                                                                                            gAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGa                 <  2:1275972/150‑1 (MQ=255)
                                                                                                                              gaaacgggCCGGCCTTCCGAAGCCGGCCCATTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGAcc               >  1:1595808/6‑150 (MQ=255)
                                                                                                                                  tGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGCTACGACGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCActc           <  1:1572716/150‑1 (MQ=255)
                                                                                                                                    ggccggccTTCTGAAGCCGGCCCGTTTCGCTATTTAGAAGTTGa                                                                                                                   >  1:882414/1‑44 (MQ=14)
                                                                                                                                    ggccggccTTCTGAAGCCGGCCCGTTTCGCTATTTAGAAGTTGa                                                                                                                   <  2:882414/44‑1 (MQ=14)
                                                                                                                                      ccGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCACTCTGaa       >  1:189621/1‑150 (MQ=255)
                                                                                                                                           cTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGATCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCACTCTGAATAAc   <  1:1212049/150‑1 (MQ=255)
                                                                                                                                            ttCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCACTCTGAATAACa  >  2:401887/1‑150 (MQ=255)
                                                                                                                                                gAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTgaaacgaaa                               <  2:509622/117‑1 (MQ=255)
                                                                                                                                                gAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTgaaacgaaa                               >  1:509622/1‑117 (MQ=255)
                                                                                                                                                 aaGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTgaaacgaaa                               >  2:1386377/1‑116 (MQ=255)
                                                                                                                                                 aaGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTgaaacgaaa                               <  1:1386377/116‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
GCAGCGATGGCTTAATCGGTACTGGCCTTTCGTTGCAAAAGATGAAACGGCTGGGTGCAGGCGTCCGAACGCCTGGTGCAATTCAGGCGGAAAACAGCGAACCCTTCGATTGGGAATAAATGGCGAAAAATGGGCCGGCTTCAGAAGGCCGGCCCGTTTCGCTATTTAGAAGTTGAAGGTGTAGCGCATGATCAGGCGATTTTCATTAGCGCTGTCAGCGTAGGTCGAGCGGTTCGTGGAGTTACGCCACTGAAACGAAACACCCTTCGCCGGACCACTCTGAATAACA  >  NC_002947/4151840‑4152128

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