Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,544,531 C→A 17.8% intergenic (+48/+66) mdfA → / ← ampG multidrug MFS transporter/muropeptide permease AmpG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,544,5310CA17.8% 239.0 / 32.1 90intergenic (+48/+66)mdfA/ampGmultidrug MFS transporter/muropeptide permease AmpG
Reads supporting (aligned to +/- strand):  ref base C (38/36);  new base A (8/8);  total (46/44)
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.58e-01

CTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAG  >  NC_002947/1544382‑1544677
                                                                                                                                                     |                                                                                                                                                  
ctgcATTTTGCCGTGGCGGCGGGGGGGGCGGGGGCGTTGGCGTTGTGGCCCCTGTGGTTTGTGGCGAGTGGGCGGCCTGTGGAAGGGGCTGTAGCTCGATAATCGGTATTCCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTacac                                                                                                                                                    <  1:481992/150‑1 (MQ=255)
 tgcATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAcaca                                                                                                                                                   <  2:66184/150‑1 (MQ=255)
 tgcATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAcaca                                                                                                                                                   >  1:66184/1‑150 (MQ=255)
   cATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGGGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAgg                                                                                                                                                 <  2:1007501/150‑1 (MQ=255)
     ttttGCCCTGGCGGCGGGGGTGGCGGGGGGGTTGGCGCTGTGTGTGGTATGTTTTTGGGGGGGTGTGCGGCGAGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGgc                                                                                                                                               >  1:989233/1‑150 (MQ=37)
      tttGCCCTGGGGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGgcc                                                                                                                                              <  1:322879/150‑1 (MQ=255)
            cTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTAGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTCTGCGTATAAACCCGCTCCTACACAGGgccgcgccg                                                                                                                                        >  2:612047/1‑150 (MQ=255)
             tGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGgccgcgccgc                                                                                                                                       >  2:210659/1‑150 (MQ=255)
             tGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCggcgcggccc                                                                                                                                       >  2:1588077/1‑146 (MQ=255)
               gcggcgGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCtt                                                                                                                                     <  1:1510080/150‑1 (MQ=255)
                gggcggGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCttt                                                                                                                                    <  1:2304726/149‑1 (MQ=255)
                   cggGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTa                                                                                                                                 <  1:1530625/150‑1 (MQ=255)
                    ggggggggCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAg                                                                                                                                <  1:620533/150‑1 (MQ=255)
                      ggggggCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTTGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGa                                                                                                                              <  1:1256447/150‑1 (MQ=25)
                      ggggggCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGa                                                                                                                              <  1:280871/150‑1 (MQ=255)
                       ggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAg                                                                                                                             <  1:517795/150‑1 (MQ=255)
                       ggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAg                                                                                                                             >  2:435362/1‑150 (MQ=255)
                          ggCGGGGGCGTTGGCGTTGTGGCCGCTGTGGTTTGTGGTGGGGGGGCGGCGTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCgg                                                                                                                          <  2:2193512/150‑1 (MQ=255)
                          ggCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCgg                                                                                                                          <  1:1320446/150‑1 (MQ=255)
                          ggCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCgg                                                                                                                          <  1:1273730/150‑1 (MQ=37)
                          ggCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTAGTGGTCGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCgg                                                                                                                          <  1:612047/150‑1 (MQ=255)
                           gCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCggg                                                                                                                         <  1:435362/150‑1 (MQ=255)
                           gCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCggg                                                                                                                         <  1:1588077/150‑1 (MQ=25)
                                   gTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCcg                                                                                                                 >  1:2064407/1‑150 (MQ=255)
                                    ttGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCcgc                                                                                                                >  1:351894/1‑150 (MQ=255)
                                       gcgcTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGaa                                                                                                             <  1:210659/150‑1 (MQ=255)
                                         gcTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAaga                                                                                                           <  2:351894/150‑1 (MQ=255)
                                              ggCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGt                                                                                                                        <  2:835550/132‑1 (MQ=21)
                                              ggCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGt                                                                                                                        >  1:835550/1‑132 (MQ=21)
                                               gCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCgg                                                                                                                          <  1:1669249/129‑1 (MQ=255)
                                               gCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCgg                                                                                                                          >  2:1669249/1‑129 (MQ=255)
                                                    tgtgGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAg                                                                                                <  1:696287/150‑1 (MQ=255)
                                                    tgtgGTTGGTGGGGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAg                                                                                                <  1:2185012/150‑1 (MQ=255)
                                                                 gggtgggCTGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTATTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCt                                                                                   >  1:207299/1‑150 (MQ=255)
                                                                      ggcggcCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCCATACCAAGCTTCAgg                                                                              >  2:943967/1‑150 (MQ=255)
                                                                             tGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTc                                                                       >  2:647869/1‑150 (MQ=255)
                                                                               gggAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTAc                                                                                                                    <  2:2190629/103‑1 (MQ=255)
                                                                               gggAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTAc                                                                                                                    >  1:2190629/1‑103 (MQ=255)
                                                                                ggAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGc                                                                    <  1:647869/150‑1 (MQ=255)
                                                                                 gAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAagag                                                                                                          >  2:723537/1‑111 (MQ=255)
                                                                                 gAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAagag                                                                                                          <  1:723537/111‑1 (MQ=255)
                                                                                 gAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGaa                                                                                                             <  1:715626/108‑1 (MQ=14)
                                                                                 gAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGaa                                                                                                             >  2:715626/1‑108 (MQ=14)
                                                                                       gCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTc                                                                                 <  1:1167106/130‑1 (MQ=255)
                                                                                       gCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTc                                                                                 >  2:1167106/1‑130 (MQ=255)
                                                                                       gCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAg                                                             >  1:882790/1‑150 (MQ=255)
                                                                                       gCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAg                                                             >  1:1102373/1‑150 (MQ=255)
                                                                                         tGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTGCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCGATACAAAGCTCCAGGGCCGCTCAGCCACCGAgcg                                                           >  2:415702/1‑148 (MQ=255)
                                                                                         tGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACACGGCCGCGCCGCGTTGTAGGAGCGGGGTGCGGCGCGAAGAGGCCAGCACAGGGAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTg                                                           >  1:1552791/1‑150 (MQ=38)
                                                                                                  aTAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTa                                                                                                                                 >  1:849541/1‑71 (MQ=255)
                                                                                                  aTAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTa                                                                                                                                 <  2:849541/71‑1 (MQ=255)
                                                                                                     aTCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCa                                                                      <  1:1427203/127‑1 (MQ=255)
                                                                                                     aTCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCa                                                                      >  2:1427203/1‑127 (MQ=255)
                                                                                                      tttgTATTGCCTGTGCTGGCCTCTTCGCTGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGc                                                                                                                            <  1:1151187/69‑1 (MQ=14)
                                                                                                      tttgTATTGCCTGTGCTGGCCTCTTCGCTGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGc                                                                                                                            >  2:1151187/4‑72 (MQ=14)
                                                                                                      tCGGTATTGCCGGGGCCGGCCTCTTCGCGTGTAAACCCGCTCCTACACAGAGCCGCGCCGCTTTGTAGGAACGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCATCAGGGCCGCTCCGCCACCGAGTGGTCGTCCGCCAcc                                              >  1:2192794/1‑150 (MQ=37)
                                                                                                          tATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCg                                          <  2:882790/150‑1 (MQ=255)
                                                                                                                ctgtgctggccTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTccccc                                    >  1:1978315/8‑150 (MQ=17)
                                                                                                                ctgtgctggccTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTccccc                                    >  1:1774773/8‑150 (MQ=17)
                                                                                                                  gggCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCcccgc                                  <  2:2064407/150‑1 (MQ=255)
                                                                                                                     ccggacTCTTGGCGGGTATACCCGTGACTGCACAGGGCCGGGCCGCTGTGTAGGAGAGGGTTTAGCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCTGCGAGTGGTCGTCCGCCACCGGCGTACcccgcccg                               <  2:989233/150‑1 (MQ=37)
                                                                                                                      tggccTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCcccgcccgc                              <  2:1978315/149‑2 (MQ=21)
                                                                                                                       ggccTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCTCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCCATCCGTCCTCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGc                             >  1:1895306/1‑150 (MQ=37)
                                                                                                                          ctattCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTc                          <  2:207299/147‑1 (MQ=255)
                                                                                                                           tctTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAg                                                                                     <  2:378637/90‑1 (MQ=255)
                                                                                                                           tctTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAg                                                                                     >  1:378637/1‑90 (MQ=255)
                                                                                                                           tctTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTcc                         >  1:233376/1‑150 (MQ=21)
                                                                                                                             ttCGCGGGTAAACCCGCTCCTGCACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTg                       >  1:2180309/1‑150 (MQ=255)
                                                                                                                             ttCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGCTTACCTACGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTg                       >  2:1735487/1‑150 (MQ=255)
                                                                                                                               cgcgGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGCAATACCAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGAt                     <  2:889773/150‑1 (MQ=17)
                                                                                                                                gcgGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGcaca                                                                                                 <  2:901068/73‑1 (MQ=255)
                                                                                                                                gcgGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGcaca                                                                                                 >  1:901068/1‑73 (MQ=255)
                                                                                                                                gcgGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATg                    <  2:1102373/150‑1 (MQ=255)
                                                                                                                                 cgGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTcc                                                    <  2:201398/117‑1 (MQ=255)
                                                                                                                                 cgGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTcc                                                    >  1:201398/1‑117 (MQ=255)
                                                                                                                                        aCCCGCCCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTgg                                                          >  2:2106761/1‑104 (MQ=255)
                                                                                                                                        aCCCGCCCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTgg                                                          <  1:2106761/104‑1 (MQ=255)
                                                                                                                                         cccGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGCTTACCCGCGAAGAGGCCAGCACACGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTGCGGCACCGGCGTCTCCCGACCGGCGTCGTGATGCCACTGCac           >  2:1086346/1‑149 (MQ=38)
                                                                                                                                         cccGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCt                        <  1:1692459/137‑1 (MQ=25)
                                                                                                                                         cccGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCt                        >  2:1692459/1‑137 (MQ=25)
                                                                                                                                          ccGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCcgc                                                                                                                <  1:1466502/48‑1 (MQ=255)
                                                                                                                                          ccGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCcgc                                                                                                                >  2:1466502/1‑48 (MQ=255)
                                                                                                                                          ccGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACATAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg          >  2:299706/1‑150 (MQ=255)
                                                                                                                                          ccGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCCgg          >  2:1442304/1‑150 (MQ=255)
                                                                                                                                          ccGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg          >  1:969437/1‑150 (MQ=255)
                                                                                                                                            gCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGc        >  1:2264834/1‑150 (MQ=255)
                                                                                                                                            gCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGc        >  1:546321/1‑150 (MQ=255)
                                                                                                                                            gCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGc        >  2:2056829/1‑150 (MQ=37)
                                                                                                                                              tCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCcccgcc                                 <  1:1703229/123‑1 (MQ=255)
                                                                                                                                              tCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCcccgcc                                 >  2:1703229/1‑123 (MQ=255)
                                                                                                                                              tCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCga      >  1:1075842/1‑150 (MQ=255)
                                                                                                                                                cTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAg                                                                                     <  1:1721668/69‑1 (MQ=255)
                                                                                                                                                cTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAg                                                                                     >  2:1721668/1‑69 (MQ=255)
                                                                                                                                                  acacAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGc                                                                                                    <  1:724708/52‑1 (MQ=255)
                                                                                                                                                  acacAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGc                                                                                                    >  2:724708/1‑52 (MQ=255)
                                                                                                                                                  acacAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcag  >  1:1555399/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAG  >  NC_002947/1544382‑1544677

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