Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,224,235 G→C 10.3% G104G (GGC→GGG PP_1962 ← site‑specific recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,224,2350GC10.3% 187.5 / 10.3 68G104G (GGC→GGGPP_1962site‑specific recombinase
Reads supporting (aligned to +/- strand):  ref base G (28/33);  new base C (3/4);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

GCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCT  >  NC_002947/2224087‑2224379
                                                                                                                                                       |                                                                                                                                                
gctaTAAGGAACTCCAAAATCGCAACATGAAAAACATGCCTGAATATTTACTCTTTTTTCCACTTGCTTTTGCTTCGTCGAGGGCCATATTAAGCATCTTATTTGTTATGAGATCCAAGTTCTCAGGAATCGCTGAAAGCCTATTgggg                                                                                                                                                  <  8:354797/146‑1 (MQ=255)
  gATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGc                                                                                                                                                 <  7:130980/148‑1 (MQ=255)
  gATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGAGTCGTCGAGGGCCCGATTAAGCATCGGATTAGTTATTAGAAAAAAGTGATAATGAATCGCTGAAAGCATATTGGGGc                                                                                                                                                 >  3:206355/1‑148 (MQ=255)
    tAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCAGTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGcccc                                                                                                                                               >  8:159422/1‑148 (MQ=255)
        gAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAg                                                                                                                                          <  8:240764/149‑1 (MQ=255)
        gAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAg                                                                                                                                          >  1:5921/1‑149 (MQ=255)
          aaGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAgg                                                                                                                                         >  6:21358/1‑148 (MQ=255)
             ccAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGCaaa                                                                                                                                     >  6:271382/1‑149 (MQ=255)
               aaaaTCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGCCCCAATAggctttca                                                                                                                                   <  7:159422/149‑7 (MQ=255)
                aaaTCCCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGACGGCCCGATTAACCATCGCCATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGCCCCAATAggctttca                                                                                                                                   <  1:303557/149‑7 (MQ=255)
                  atagCAACATGAAAAACATTCCTTTATCATGACTAGGTTGTCCAATTGCTTTTGATGCTTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCg                                                                                                                                <  4:206355/146‑1 (MQ=255)
                     gCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGAt                                                                                                                              <  3:3304/148‑1 (MQ=255)
                     gCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATa                                                                                                                             <  3:188642/149‑1 (MQ=255)
                       aaCATCAAAAAAAGGCCTGAATAATGCCTCGGATGAACACTTGCTTTTGCGTCGTCGAGGTCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGTAATCGCTGAATGCCTATTGGTGCCCAATAGGTAAATTTCGATaaa                                                                                                                           >  6:281240/1‑149 (MQ=255)
                        aCATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAAc                                                                                                                          >  6:75924/1‑149 (MQ=255)
                        aCATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGCAAATTTCGATAAAc                                                                                                                          <  5:271382/149‑1 (MQ=255)
                          aTGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAACCTATTGGGGCCCAATAGGTAAATTTCGATAAAcc                                                                                                                         >  3:251756/1‑148 (MQ=255)
                              aaaaCTGGCCTTAATTATTACTCTTTTGTACATTTTTTTTTGCGTTGTCGTGGCCCATCTTAAAGAATCGTATCTCCATATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGAt                                                                                                                      <  1:307970/149‑1 (MQ=255)
                                      ccTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAAtt                                                                                                                                   <  2:13018/126‑1 (MQ=255)
                                      ccTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAAtt                                                                                                                                   >  1:13018/1‑126 (MQ=255)
                                       cTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGAt                                                                                                           >  4:366547/1‑149 (MQ=255)
                                        tGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTga                                                                                                              >  1:345377/1‑145 (MQ=255)
                                        tGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTga                                                                                                              <  2:345377/145‑1 (MQ=255)
                                        tGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGAtt                                                                                                          >  7:260060/1‑149 (MQ=255)
                                          aataatGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTc                                                                                                        >  4:333461/1‑149 (MQ=255)
                                              atGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCtt                                                                                                    >  8:128147/1‑149 (MQ=255)
                                               tGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCttt                                                                                                   <  2:210397/149‑1 (MQ=255)
                                                  cTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGga                                                                                                >  6:52604/1‑149 (MQ=255)
                                                   tCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGCGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGgat                                                                                               >  7:261883/1‑149 (MQ=255)
                                                     cgtgggccACCTGCTGTTGCGGCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTAATGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGgat                                                                                               <  8:77809/141‑1 (MQ=255)
                                                      gTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGgatgat                                                                                            <  8:260060/149‑1 (MQ=255)
                                                       ttGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATa                                                                                           <  5:21358/149‑1 (MQ=255)
                                                            cACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAg                                                                                      >  6:249670/1‑149 (MQ=255)
                                                               ttGCATTTTCGTCGTCGAGGGCCCGATTTAGAATCGGATTTGTTATTAGATCAAAGTGCTCATGATTCGCACAAAGAATATTGGGGCCCAATAGGTAAATTTCGATAAACCGATATTCGTGAGATTTCGCTTTGTATGATACGAAGGaa                                                                                   <  8:344165/149‑1 (MQ=255)
                                                                 gCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCTGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTAATgtgg                                                                                                                                                  <  8:236730/84‑1 (MQ=255)
                                                                 gCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAAc                                                                                  >  6:348421/1‑148 (MQ=255)
                                                                 gCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAAc                                                                                  <  5:348421/148‑1 (MQ=255)
                                                                 gCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTAATgggg                                                                                                                                                  >  7:236730/1‑84 (MQ=255)
                                                                   ttttGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGaa                                                                               >  5:142265/1‑149 (MQ=255)
                                                                   ttttGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGaa                                                                               >  5:166110/1‑149 (MQ=255)
                                                                       gCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATAc                                                                            <  7:128147/148‑1 (MQ=255)
                                                                         gtcgtcGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGAt                                                                         <  3:366547/149‑1 (MQ=255)
                                                                         gtcgtcGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGAt                                                                         >  3:304724/1‑149 (MQ=255)
                                                                             tcgAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTg                                                                       >  3:212202/1‑147 (MQ=255)
                                                                               gAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAgg                                                                   <  5:75924/149‑1 (MQ=255)
                                                                               gAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAgg                                                                   <  7:283633/149‑1 (MQ=255)
                                                                                 gggCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGa                                                                 <  2:5921/149‑1 (MQ=255)
                                                                                  ggCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAg                                                                <  5:249670/149‑1 (MQ=255)
                                                                                             aGCATCGGATTTGTTATGTGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAg                                                       <  8:261883/148‑1 (MQ=255)
                                                                                             aGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAAc                                                                                                                          >  2:235876/1‑81 (MQ=255)
                                                                                             aGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAAc                                                                                                                          <  1:235876/81‑1 (MQ=255)
                                                                                              gCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAggg                                                     <  6:166110/149‑1 (MQ=255)
                                                                                                       ttttttATTAGATCTAATTGGTCAGGAATCGCTGAAAGCCTATTGTGGCCAAATAGGTAAATTTCTATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAAt                                              <  1:170569/148‑1 (MQ=255)
                                                                                                          tGTTATGAGAACCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATgg                                            <  2:108723/148‑1 (MQ=255)
                                                                                                            ttATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTa                                          <  6:142265/148‑1 (MQ=255)
                                                                                                             tATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTaaa                                        >  4:163534/1‑149 (MQ=255)
                                                                                                                  gaTCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAAtt                                   >  4:78771/1‑149 (MQ=255)
                                                                                                                     ccAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAAAAGATAGTCGTGAGAATACGCTTTGGATGATACGAAGGAACGAATACGATAGAAGGGAGTAAATGAAGGCAGCAAATGGTAAACAATTGca                                >  6:280691/1‑148 (MQ=255)
                                                                                                                            cTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTc                          <  4:304724/148‑1 (MQ=255)
                                                                                                                             tCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAAtt                                   <  1:107641/138‑1 (MQ=255)
                                                                                                                             tCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAAtt                                   >  2:107641/1‑138 (MQ=255)
                                                                                                                                         atttaccTATTGGGCCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCAcc            <  3:131368/144‑1 (MQ=255)
                                                                                                                                          tttaccTATTGGGCCCCAATAGGTAAATTTCGATAGACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGa                                                    >  2:14719/5‑108 (MQ=255)
                                                                                                                                          aaaGCCTATTGGTGGACAATAGTTAAATTTCTATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCAcc            <  5:281240/148‑1 (MQ=255)
                                                                                                                                          aaaGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCt           >  4:275748/1‑149 (MQ=255)
                                                                                                                                           ttaccTATTGGGCCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTg          >  6:327783/4‑149 (MQ=255)
                                                                                                                                           aaGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTg          <  5:52604/149‑1 (MQ=255)
                                                                                                                                             gCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTgg         >  5:207722/1‑148 (MQ=255)
                                                                                                                                               atATTGGGGCCGAATAGGTAAATTTCGATAAACCGATACTCGTGAGATTTCGCTTTTGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTAt      <  4:251756/148‑1 (MQ=255)
                                                                                                                                                tATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAgg                                                                   <  2:361987/87‑1 (MQ=255)
                                                                                                                                                 aTTGGGCCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGAt                                                                                                           >  4:303712/1‑46 (MQ=255)
                                                                                                                                                 aTTGGGCCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGAt                                                                                                           <  3:303712/46‑1 (MQ=255)
                                                                                                                                                    ggggCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTAtctct  <  3:275748/148‑1 (MQ=255)
                                                                                                                                                    gcggCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTAtctct  <  3:163534/146‑1 (MQ=255)
                                                                                                                                                       gCCCAATAGGTAAATTTCGATAAACCGATAGTCATGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGt                                           <  1:100953/104‑1 (MQ=255)
                                                                                                                                                       gCCCAATAGGTAAATTTCGATAAACCGATAGTCATGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGt                                           >  2:100953/1‑104 (MQ=255)
                                                                                                                                                       |                                                                                                                                                
GCGATAAGGAAAGCCAAAATCGCAACATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGGTAAATTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCT  >  NC_002947/2224087‑2224379

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