Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,560,216 A→T 5.9% L111Q (CTG→CAG)  hisZ ← ATP phosphoribosyltransferase regulatory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,560,2160AT 5.9% 260.7 / 4.4 84L111Q (CTG→CAG) hisZATP phosphoribosyltransferase regulatory subunit
Reads supporting (aligned to +/- strand):  ref base A (47/32);  new base T (2/3);  total (49/35)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

TGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAGGCGGCCGGACTGCGGGTCCACTACCTTGAAGGTGCGCTGATCCAGG  >  NC_002947/5560067‑5560365
                                                                                                                                                     |                                                                                                                                                     
tgtgAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTgca                                                                                                                                                       <  2:906105/150‑1 (MQ=255)
tgtgAAAATCCGGCACATGCGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTgca                                                                                                                                                       >  2:151106/1‑150 (MQ=255)
    aaCATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCAc                                                                                                                                                   >  1:1754390/1‑150 (MQ=255)
    aaCATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCAc                                                                                                                                                   >  2:286699/1‑150 (MQ=255)
       aTCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACACGGCACGTGGCTGAGCGTGCAGCACGCt                                                                                                                                                >  1:942800/1‑150 (MQ=255)
        tCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTg                                                                                                                                               >  2:2089546/1‑150 (MQ=255)
         ccGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCTGCACGCtca                                                                                                                                              <  1:1392705/150‑3 (MQ=255)
                aTCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATa                                                                                                                                       <  1:2235765/150‑1 (MQ=255)
                aTCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATa                                                                                                                                       >  2:2374866/1‑150 (MQ=255)
                 tCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATCTGGAGTACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAg                                                                                                                                      <  1:126406/150‑1 (MQ=255)
                  cGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGc                                                                                                                                     <  1:2405499/150‑1 (MQ=255)
                   ggCCAGTTACAGCGTCGAGAGCATAAGGCTGATGACTTCGCCATCGCTGGTAGGGCTGGCAGCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGca                                                                                                                                    >  2:398423/1‑150 (MQ=255)
                     ccAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGcaca                                                                                                                                  >  1:1519925/1‑150 (MQ=255)
                     ccAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGcaca                                                                                                                                  >  1:1310182/1‑150 (MQ=255)
                              gCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCg                                                                                                                         >  2:430355/1‑150 (MQ=255)
                               cGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCgg                                                                                                                        <  1:520620/150‑1 (MQ=255)
                                 tGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGc                                                                                                                      <  1:2086721/150‑1 (MQ=255)
                                   gatctcATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGAACAGTTCGGCACCCAGCTGGATCGCACAGCACGCGCACAAAGGGCCACAGTGCCCACCGCTCACCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCt                                                                                                                    <  1:1252706/145‑1 (MQ=255)
                                   gCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCt                                                                                                                    >  1:2260661/1‑150 (MQ=255)
                                    cGAGCCACAGGCAGAGGATTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCTAGGGCTGCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCtt                                                                                                                   <  1:151106/150‑1 (MQ=255)
                                        cATCAGGCTGATGACTTCGACATCGCTGGCAGGGCTGGCATCGCCGGACAGTTCGGCACCCGGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACg                                                                                                               >  2:637562/1‑150 (MQ=255)
                                         aTCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCttagatc                                                                                                              >  2:1252706/1‑145 (MQ=255)
                                            aGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGgcgc                                                                                                           <  2:2073842/150‑1 (MQ=255)
                                                tgatgaCTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAgtg                                                                                                       >  2:1527148/1‑150 (MQ=255)
                                                 gatgaCTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAgtgt                                                                                                      >  2:412431/1‑150 (MQ=255)
                                                     aCTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGGTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTgc                                                                                                  >  2:945390/1‑150 (MQ=255)
                                                       ttCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCTGCACGCtca                                                                                                                                              <  1:1212356/104‑3 (MQ=255)
                                                       ttCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCTGCACGCtca                                                                                                                                              >  2:1212356/1‑102 (MQ=255)
                                                       ttCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTgcgc                                                                                                >  1:2271579/1‑150 (MQ=255)
                                                         ccacacCGCTAGTAGCGCTCGCATAGCCGTACAGTGCGGCTCCCAGCTGGATCGGGCTGCCCGAGCTGGACAATGCACGTGGCTGAGCGTGCAGCACGTTCTCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCAt                                                                                              <  2:942800/148‑1 (MQ=255)
                                                              tCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATa                                                                                         >  1:1767018/1‑150 (MQ=255)
                                                              tCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATa                                                                                         >  2:1879039/1‑150 (MQ=255)
                                                                gCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATAcg                                                                                       >  2:757339/1‑150 (MQ=255)
                                                                gCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATAcg                                                                                       >  1:447828/1‑150 (MQ=255)
                                                                 cTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTTACGGCGCAGTGTCTGCGCATCGCTAcgc                                                                                      >  2:988567/1‑150 (MQ=255)
                                                                  tGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATAcgcg                                                                                     >  2:1488420/1‑150 (MQ=255)
                                                                  tGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATAcgcg                                                                                     <  2:1486157/150‑1 (MQ=255)
                                                                  tGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATAGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGACGTCATAGCACAGCCGGCTCGGGCCGTCACGGCGCAGTGTGTGCGCATAGATAcgcg                                                                                     >  1:398200/1‑150 (MQ=255)
                                                                            ggCATCGCCGTAGAGTGCCGTGTACAGCTGGATCGGGCTGTGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGGTCGGGCCTTGACGGCGCAGTGTGTGCGCATCGATACGCGCGACCTGCgg                                                                           <  2:26169/150‑1 (MQ=255)
                                                                             gCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAgtgt                                                                                                      <  1:732472/122‑1 (MQ=255)
                                                                             gCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAgtgt                                                                                                      >  2:732472/1‑122 (MQ=255)
                                                                               aTCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCTGCACGCtca                                                                                                                                              >  2:688397/1‑78 (MQ=255)
                                                                               aTCGCCGCACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCTGCACGCtca                                                                                                                                              <  1:688397/80‑3 (MQ=255)
                                                                                  gCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGaa                                                                     <  2:731702/150‑1 (MQ=255)
                                                                                  gCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGaa                                                                     >  1:19797/1‑150 (MQ=255)
                                                                                         aGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGAt                                                                                          >  1:1673735/1‑122 (MQ=255)
                                                                                         aGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGAt                                                                                          <  2:1673735/122‑1 (MQ=255)
                                                                                         aGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACgcgc                                                                                    <  2:1369307/128‑1 (MQ=255)
                                                                                         aGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACgcgc                                                                                    >  1:1369307/1‑128 (MQ=255)
                                                                                         aGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTACGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGgcg                                                              >  2:452174/1‑150 (MQ=255)
                                                                                           ttCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGcgcg                                                            >  2:2344496/1‑150 (MQ=255)
                                                                                               gCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAg                                                        >  2:2283420/1‑150 (MQ=255)
                                                                                               gCACCCAGCTGGATCGGGCTGCGCGAGCTAAACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAg                                                        >  1:1642796/1‑150 (MQ=255)
                                                                                                   ccAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATAc                                                                                        <  2:547228/114‑1 (MQ=255)
                                                                                                   ccAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATAc                                                                                        >  1:547228/1‑114 (MQ=255)
                                                                                                   ccAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCa                                                    >  2:812476/1‑150 (MQ=255)
                                                                                                   ccAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCa                                                    >  2:1596637/1‑150 (MQ=255)
                                                                                                   ccAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCAAGGAGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCa                                                    >  1:1998340/1‑150 (MQ=255)
                                                                                                           aTCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAggcgg                                            >  2:2138169/1‑150 (MQ=255)
                                                                                                              gggCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAGGCGGCCg                                         >  2:2094995/1‑150 (MQ=255)
                                                                                                              gggCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAGGCGGCCg                                         <  2:1310182/150‑1 (MQ=255)
                                                                                                                gCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAGGCGGCCgtt                                       <  1:696387/150‑3 (MQ=255)
                                                                                                                gCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAGGCGGCCGGa                                       >  2:1216482/1‑150 (MQ=255)
                                                                                                                    cgcgAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAGGCGGCCGGACTGc                                   >  1:2178510/1‑150 (MQ=255)
                                                                                                                     gcgAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAGGCGGCCGGACTGCg                                  >  1:2354451/1‑150 (MQ=255)
                                                                                                                      cgAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAGGCGGCCGGACTGCgg                                 <  1:398423/150‑1 (MQ=255)
                                                                                                                      cgAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGAATTGAAGCTCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAGGCGGCCGGACTGCgg                                 <  1:149737/150‑1 (MQ=255)
                                                                                                                          gtagACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAGGCGGCCGGACTGCGGGTcc                             >  1:2407008/4‑150 (MQ=255)
                                                                                                                            ggACAAGGCACGTGGCTGAGCGTGCAGCACGATGGCGGCATAGCACAGTCGGCTCGCGGCTTGACGGAGGAGTGTGTGCGCATCGATACCCGCTACCGGGGGGGTGCAGTCGGCGCGGAAGCCCATCAGGCGGCCGGACTGCGGGTccgt                           <  2:684256/150‑3 (MQ=255)
                                                                                                                             gACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAGGCGGCCGGACTGCGGGTCCACt                          >  2:727441/1‑150 (MQ=255)
                                                                                                                              aCAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAGGCGGCCg                                         <  2:2366756/134‑1 (MQ=255)
                                                                                                                              aCAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAGGCGGCCg                                         >  1:2366756/1‑134 (MQ=255)
                                                                                                                               cAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAGGCGGCCGGACTGCGGGTCCACTAc                        <  1:1527148/150‑1 (MQ=255)
                                                                                                                                 aGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAGGCGGCCGGACTGCGGGTCCACTACCt                      <  1:430355/150‑1 (MQ=255)
                                                                                                                                  tgCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAGGCGGCCGGACTGCGGGTCCACTACCtt                     >  1:1909344/2‑150 (MQ=255)
                                                                                                                                  ggCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAGGCGGCCGGACTGCGGGTCCACTACCtt                     <  2:324302/150‑1 (MQ=255)
                                                                                                                                    cACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAGGCGGCCGGACTGCGGGTCCACTACCTTGa                   >  2:1798241/1‑150 (MQ=255)
                                                                                                                                     aCGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAGGCGGCCGGACTGCGGGTCCACTACCTTGaa                  <  2:2407008/150‑1 (MQ=255)
                                                                                                                                      cGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAGGCGGCCGGACTGCGGGTCCACTACCTTGAAg                 <  1:1596637/150‑1 (MQ=255)
                                                                                                                                      cGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAGGCGGCCGGACTGCGGGTCCACTACCTTGAAg                 <  1:757339/150‑1 (MQ=255)
                                                                                                                                      cGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACACCCGGTTCGCACCTTCACGGGGCAGCGTGTCCGCATCGACACGCGCTACCTGTGGCGTGAAGTCAGCGCCGACGCCAAGCATGCGGCCGAAATGCGGGTGCACTACCGTGAAg                 <  2:897483/150‑1 (MQ=255)
                                                                                                                                       gTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGgc                                                               <  1:1662577/103‑1 (MQ=255)
                                                                                                                                       gTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGgc                                                               >  2:1662577/1‑103 (MQ=255)
                                                                                                                                               gCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAGGCGGCCGGACTGCGGGTCCACTACCTTGAAGGTGCGCTGa        <  1:945390/150‑1 (MQ=255)
                                                                                                                                               gCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCACGGAAGCCCATCAGGCGGCCGGACTGCGGGTCCACTACCTGGGAGGTGCGCTGa        >  1:41021/1‑150 (MQ=255)
                                                                                                                                                 gTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAGGCGGCCGGACTGCGGGTCCACTACCTTGAAGGTGCGCTGATc      <  1:988567/150‑1 (MQ=255)
                                                                                                                                                   gcagcaCGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAGGCGGCCGGACTGCGGGTCCACTACCTTGAAGGTGCGCTGATCCa    >  1:797996/1‑150 (MQ=255)
                                                                                                                                                    cagcaCGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAGGCGGCCGGACTGCGGGTCCACTACCTTGAAGGTGCGCTGATCCAg   >  2:165428/1‑150 (MQ=255)
                                                                                                                                                     agcaCGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGAGCCAGATGCGGCCGGACTGCGGGTCCACTACCTTGAAGGTGCGCTGATCCAgg  >  1:1352748/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTGAACATCCGGCACATCGGCCAGTTGCAGCGTGGCGAGCATCAGGCTGATGACTTCGACATCGCTGGTAGGGCTGGCATCGCCGTACAGTTCGGCACCCAGCTGGATCGGGCTGCGCGAGGTGGACAAGGCACGTGGCTGAGCGTGCAGCACGCTGCCGGCATAGCACAGCCGGCTCGGGCCTTCACGGCGCAGTGTGTGCGCATCGATACGCGCTACCTGCGGGGTGAAGTCGGCGCGGAAGCCCATCAGGCGGCCGGACTGCGGGTCCACTACCTTGAAGGTGCGCTGATCCAGG  >  NC_002947/5560067‑5560365

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