Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,571,097 A→G 14.7% L160P (CTC→CCC)  yjiR ← putative DNA‑binding transcriptional regulator/putative aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,571,0970AG14.7% 162.3 / 2.0 73L160P (CTC→CCC) yjiRputative DNA‑binding transcriptional regulator/putative aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (9/50);  new base G (10/2);  total (21/52)
Fisher's exact test for biased strand distribution p-value = 9.88e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTCCAGCGCCAGTTCCAGCGCTTCAACGCTGATGCCAGTTTCTGGATCGGTTGGGATTTCAATCACTTTCACGCCCATGCCGCGCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCG  >  NC_000913/4570863‑4571322
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                      
tGTTCCAGCGCCAGTTCCAGCGTTTCAACGCTGATGCCAGTTTCTGGATCGGTTGGGATTTCAATCACTTTCACGCCCATGCCGCGCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCAt                                                                                                                                                                                                                                    <  1:117612/239‑1 (MQ=255)
    ccAGCGCCAGTTCCAGCGCTCCAACGCTGATGCCAGTTTCTGGATCGGTTGGGATTTCAATCACTTTCACGCCCATGCCGCGCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAgg                                                                                                                                                                                                                                <  2:794641/239‑1 (MQ=255)
         gcCAGTTCCAGCGCTTCAACGCTGATGCCAGTTTCTGGATCGGTTGGGATTTCAATCACTTTCACGCCCATGCCGCGCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCggg                                                                                                                                                                                                                            <  1:248347/238‑1 (MQ=255)
              ttCCAGCGCTTCAACGCTGATGCCAGTTTCTGGATCGGTTGGGATTTCAATCACTTTCACGCCCATGCCGCGCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGGGACCACCGAACCCCCGTCTGAGCTCAAGCGGGTGATCt                                                                                                                                                                                                                      >  2:971087/1‑239 (MQ=255)
               tCCAGCGCTTCAACGCTGATGCCAGTTTCTGGATCGGTTGGGATTCCAATCACTTTCACGCCCATGCCGCGCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCGACAATATCGCCCGGTTTACCCACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTg                                                                                                                                                                                                                     <  1:1130997/239‑1 (MQ=255)
                ccAGCGCTTCAACGCTGATGCCAGTTTCTGGATCGGTTGGGATTTCAATCACTTTCACGCCCATGCCGCGCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGc                                                                                                                                                                                                                    <  1:971087/239‑1 (MQ=255)
                     gcTTCAACGCTGATGCCAGTTTCTGGATCGGTTGGGATTTCAATCACTTTCACGCCCATGCCGCGCAGCGTCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGctc                                                                                                                                                                                                                  <  1:708304/236‑1 (MQ=255)
                     gcTTCAACGCTGATGCCAGTTTCTGGATCGGTTGGGATTTCAATCACTTTCACGCCCATGCCGCGCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGctct                                                                                                                                                                                                                 <  1:880428/237‑1 (MQ=255)
                        tCACCGCGGATGCCAGTTTCTGGATCGGTTGGGATTCCAATCACTTTCACGCCCATGCCGCGCAGCATCTGCATCGAACCGTAATAACCGGGGGATCCGCCCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCACCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCaa                                                                                                                                                                                                            <  1:527408/239‑1 (MQ=255)
                                   gCCAGTTTCTGGATCGGTTGGGATTTCAATCACTTTCACGCCCATGCCGCGCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGAc                                                                                                                                                                                                 <  1:936388/239‑1 (MQ=255)
                                    caaGTTTCTGGATGGGTTGTAATTGCAATCACCGTCAGCCCCATGCGGCGCAGCATCTGCATAGAACCGTAATACCGGGGGGATTCGACCGCGACAATATCGCCCGGTTTACACGCTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGAc                                                                                                                                                                                                 <  2:3980/236‑1 (MQ=255)
                                          tCTGGATCGGTTGGGATTTCAATCACTTTCACGCCCATGCCGCGCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCGACAATATCGCCGGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTaa                                                                                                                                                                                          <  1:617409/239‑1 (MQ=255)
                                          tCTGGATCGGTTGGGATTTCAATCACTTTCACCCCCATGCCCCGCAGCATCTGCATCGAACCGTAATACCAGGGGGATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTaa                                                                                                                                                                                          <  1:627659/239‑1 (MQ=255)
                                           cTGGATCGGTTGGATTTTCAATCCCTTTCACCCCCATGCCGCGCAGCATCTGCATCGAACCGTAATACCAGGGGGATTCGCCCGCGACAATATCGCCGGTTTTACACACTGCCATTAACGCCAGCGCCATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGCCCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAAc                                                                                                                                                                                         <  2:193331/239‑1 (MQ=255)
                                                 cGGTTGGGATTTCAATCACTTTCACGCCCATGCCGCGCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCa                                                                                                                                                                                   <  1:1087123/239‑1 (MQ=255)
                                                 cGGTTGGGATTTCAACCCCTTTCACGCCCATGCCGCGCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCa                                                                                                                                                                                   <  2:845431/239‑1 (MQ=255)
                                                                      tCACGCCCATGCCGCGCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCCCCCAGCCGCCGGCGGGCATCAGGCGGGCAATCTGCTCTCGCCATACCCGCTGACCGGCTCACAACCCATAACCGCGAACGAGTTGCCGa                                                                                                                                                             >  2:109078/1‑239 (MQ=255)
                                                                      tCACGCCCATGCCGCGCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGa                                                                                                                                                             <  1:530030/239‑1 (MQ=255)
                                                                      tCACGCCCATGCCACCCAGCATCTGAACCGAACCGTAATAACAGGGGGATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGa                                                                                                                                                             <  2:182407/239‑1 (MQ=255)
                                                                           cccATGCCGCGCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTAtg                                                                                                                                                        <  2:198418/239‑1 (MQ=255)
                                                                            ccagtccGCGCAGCATCTGCATCGAACGGTAATAACAGGGGGATTCGCCGGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTAtg                                                                                                                                                        <  1:261147/233‑1 (MQ=255)
                                                                                     gcagcaTCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGcaccac                                                                                                                                               <  1:747423/238‑1 (MQ=255)
                                                                                           tCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAAACTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAAAGGGTTTCAGATTATGCTGCACCACCCGGCtt                                                                                                                                        >  1:795671/1‑239 (MQ=255)
                                                                                                   gAACCGTAATAACAGGGGGATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGCGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCCATACTCGCTGACCCGCTAACAAATCATAAACCAGAACGGGTTGCAGATTATGCTGCCCCACCCGGCTTATCcccc                                                                                                                                 >  1:1100266/1‑238 (MQ=255)
                                                                                                    acccgaaaaatgggggggATGCGCCCGCGACAATACCGCCCGTTTTAAACCCTGCCTTTAACGCCGGCGACACCGAGTTAGGGCAGCGGCGGGTGATGAGGATGTAATCGGCGGTGACCACCGGGCCGCTGTCGAGCATCAGGCGGGCAACCGGCTCTCGAAATACTCGCTGACCGGCTAACAAGTCATACCCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGTTTAGCTCCCGc                                                                                                                               <  2:402716/226‑1 (MQ=255)
                                                                                                    aaCCGTAATAACGGGGGGATCCGCCCGCGACAATATCGCCGGTTTTACACCCGGCCATAACCGCCAGCGCCACCGGGTTATGGCGGCCGCTGGGGATGATGATGTCATCGGCGGTGACCCCCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGc                                                                                                                               <  2:495721/239‑1 (MQ=255)
                                                                                                     cccGAAAAAACGGGGGGATCGGCCCGGGCAAAAACCGCCGGTTTTCCACCCGGCCATTACCGCAAGCGACTCCGGTTTATGGCAGCCGCGGGTGAGGATGATGTCATCGGCGGTGACCCCCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGcc                                                                                                                              <  2:124549/238‑1 (MQ=255)
                                                                                                         taataaCAGGGGGATTCACCCGCGACAATTCCTCCCGGTTTACACCCTGCCATTAACGCCAGCGACCTGAAGTGAGGGCGGCGGGGGGTGATGAGTATGTCATGGGCGGTGCCCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACCCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAg                                                                                                                          <  2:625546/239‑1 (MQ=255)
                                                                                                                 gggggATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGGCCACCGCGCCGCCGTCGAGCATCCGGCGGGGCATCTGCTCTCGCACTACTCGCCGACCGGCTAACAAATCATAACCCCGAACGGTTTGCAGATTATGCTGCACCACCCCGCTTCGCTCCCGCCCCCGCGGgtttc                                                                                                                >  2:489465/1‑237 (MQ=255)
                                                                                                                   gggATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCCGGCCGCCGTCGGGCATCAGGCGGGGAAACTGCTCTCGCCAAACCCCGTGACCCGCCAAAAAATCATAAACCAGAACCGGTTGCAGATTTTTGTGCACCACCCGGCGTTGGTccccaccacccggggtttga                                                                                                                >  2:859060/1‑222 (MQ=255)
                                                                                                                   gggATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGGGACCACCCCGCCGCCGGCGGGCATCAGGCGGGCAAACTCCTCTCGCACTACTCCCTCGCCGGCCAACAAATCATAACCGAGAACGGGTTTCAGATTATGCTGGACCACCCGCCCTAGCCCCCCCCACCGCGGGTTctgt                                                                                                              >  2:5292/1‑236 (MQ=255)
                                                                                                                    ggATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCCAGCCGCTGTCGAGCATCAGGGGGGCAATCTGCTCTCGCCATACTCGCCGACCGGCTAACAAGACATAACCCAGAAAGGTTTGCAGATTATGCTGCACCACCCCGCTTAGCCCCCCCCACAGCGGGTTTAgg                                                                                                              >  1:1122433/1‑238 (MQ=255)
                                                                                                                        tCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCCCCCAGCCGCCCTCGAGCATCAGGGGGGCAAACTTCTCTCGCCATACTCCCCCGCCGGCCAACAAATCCTAAACGGCAACGGTTTGCCGCATATTCTGCCCCAACCCGCCTAGCTCCCGCCCCCGCGGTTTCAAGCCggg                                                                                                         >  2:639500/1‑239 (MQ=255)
                                                                                                                            ccGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGGCCAACCAGCCGCCGTCGGGCAATAGGCGGGGAACATTCCCTCGCAATACTCCCCGGCCGGCCAAAAAAACATAAACCAGAAAGGGGTGGAGAATTATCTGGACCACCCGGCCTATCCTCCGCCCCAGACGGTTTAGGGTGGgcgc                                                                                                      >  2:125632/1‑238 (MQ=255)
                                                                                                                            ccGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGCCCCCCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCGCCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCt                                                                                                     <  1:1002972/239‑1 (MQ=255)
                                                                                                                                  ggtggataagagacagacacacTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGCTTCAGGCTGGGCGCTTCGAc                                                                                                <  1:435394/222‑1 (MQ=255)
                                                                                                                                   aaTATCGCCCGGTTTACCCACTGCATTTACCGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACAt                                                                                              <  1:859060/239‑1 (MQ=255)
                                                                                                                                      tccGCCCGGTTTACACACTGCCATTAACGCCTGCGACACCGAGTTATGGCAGCCGCTGTTGATGATGATGTCACCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCg                                                                                           <  1:976472/237‑1 (MQ=255)
                                                                                                                                               tttACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGct                                                                                                                                                                                                                   <  1:740758/113‑1 (MQ=255)
                                                                                                                                               tttACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGaa                                                                                                                                                                         <  1:637221/155‑1 (MQ=255)
                                                                                                                                                       cTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTTACCGGCTAAAAAGTCATAAACCAGAACGGGTTTCAGATTATGCTGCACCACCCCGCTTAGCTCCCCCCACAGCGGGTTTAGGCTGGGGGCTTCGCCATCCGGGGGGCTTTTGCTTAAc                                                                          >  1:729899/1‑239 (MQ=255)
                                                                                                                                                        tGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCCCCCCGCCGCCGTCGAGCATCAGGCGGGGAATCTGCTCTCGCACTACTCCCTGACCGGCTAACAAATCATAACCGCGAAAGGTTTTCAGATTATGCTGCACCACCCCGCTTAGCTCCCCCCCCCGCGGGTTCAGGCTGGGGGCCTCGCCATCCCGTGTGCTTTTTCTTTACg                                                                         >  2:376436/1‑239 (MQ=255)
                                                                                                                                                        tGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGGCGGGCATCAGGCGGGGAAACTGCTCTCGCCATACTCGCCGACCGGCTAACAAGTCATAAACCAGAAAGGGTTGCAGATTATTCTGCACCCCCCCGCCTTGCCCCCGCCCACGCGGGTTCCGGACCGGGGCTTCTACATCCCGTGTGCgtttttgttacg                                                                         >  2:340887/1‑231 (MQ=255)
                                                                                                                                                        tGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCCCCGTCGGGCAACCGGGGGGCAATCTGCCTTCTCAATACTCGCTGACCCGCTAAAAAGTCATAACCGAGAACGGGGTGTAGAGTATGTTGCGCCCCCCGGCCTTGATCCCGCCCCAGCGGTTTCCGGCTGTGGGCTTTCACATCCGGGGGGCgttttgttt                                                                            >  2:461335/1‑231 (MQ=255)
                                                                                                                                                         gCCATTAACGCCGGCGACATCGAGTTAGGGCGGCCGCTGGTGATGATGATGTCATCGGCGGTGCCCCCCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACgg                                                                       <  1:1192610/239‑1 (MQ=255)
                                                                                                                                                         gCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGGCCACCCCCCCGCTGTTGTGCCTCCGGCGGGCCAATTGCTCCTTCCATACTCGCGTCCCCGCTGACAATACAGAAACAGGAACGGTTTTTAAATTATGCAGCCACACCCGGCCTAGCTCCCGCCACACCGGGTTCCGGCTGGGTGCTTCGGCATCCGGTGTGCtttttgtttagt                                                                         >  2:336974/1‑230 (MQ=255)
                                                                                                                                                         gCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGGGACCACCGAGGCCCTGTCGGGCATCAGGGGGGCAATCCGCCCTCCCCAAACTCCCCGGCCCGCGAACACATCATAAACCAGAACAGGTTGCAGATTATGCTGCACCACCCCGCTTTGCTCCCGCCACACCGGGTTCCGGCGGGGCGCCGCCCCACCCGGCGTGCTTTTGTTTTACgg                                                                       >  2:21963/1‑239 (MQ=255)
                                                                                                                                                            aTTACCGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGCCCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGgcg                                                                                                       <  1:1112594/205‑1 (MQ=255)
                                                                                                                                                                aCGCCAGCGACATCGATTTATGGCAGCCGCTGGTGATGATGATTTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTggg                                                                                                         <  1:274380/199‑1 (MQ=255)
                                                                                                                                                                  gCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCCAGCCCCCGTCGCGCATCAGGGGGGCAATCTCTCCTCCCACTACGCGCCTAGCCGCCAACCAGACCTAAACGCGAACGGGTTTCAGGTTATGCTGCGCCACCCCGCGTTGCGCCCGCCCCAACGGGTTTAAGCCGGGGCGTCCGACCTCCGGGGCGCCTTTTCTTAAAGCGAACAACgg                                                              >  2:1030559/1‑239 (MQ=255)
                                                                                                                                                                   ccAGCGACCTCGATTTAGGGCGGCCGTTGGTGATGATGATGTCACCGGCGGGACCCCCCGAGCCGCGGTCGAGCATCAGGCGGCCAATCTGCTCTCGCAATACTCGCTGACGGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCCCCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGaa                                                            <  2:323064/239‑1 (MQ=255)
                                                                                                                                                                    cAGCGACCTCGAGTTATGGCAGCCGCTGGTGATGATGTGGTCATCGGCGGTGACCCCCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAAc                                                           <  1:21963/239‑1 (MQ=255)
                                                                                                                                                                      gCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCCCCCAGCCGCCGTCGGGCATCAGGCGGGCAATCTGCTCTCGCCATACTCGCTGACCGGCCAAAAAGTCATAACCGAGAACGGGTTGCAGATTATTCTGCCCCACCCGGCCTTGCTCCCCCCCCAGCGGGTTCCGGCTGGGCGCTTCTACATc                                                                                             >  2:1112594/1‑205 (MQ=255)
                                                                                                                                                                               gtgttgggCAGCCGCGGGTGATGATGAGTTCGCCGGCGGGGCCCCCCGGGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATg                                                 <  2:119947/232‑1 (MQ=255)
                                                                                                                                                                                  tATGGCAGCCGCGGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGcc                                             <  1:269984/239‑1 (MQ=255)
                                                                                                                                                                                   aggcgcgCCGCGGGTGAGTAAGATGCCATCGCCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAACCTGCCCTCCCAATACTCGCTGCCCGCCTACCAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCt                                            <  2:900013/233‑1 (MQ=255)
                                                                                                                                                                                   aTGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCCGTCGGGCATCCAGCGGGCAAAATGCTCTCGCAAAACTCGCTGACCGGCCAACAAATGATAAACCAGGACGGTTTGCCAGTTATGTTGCACCACCCGGCTTAACTCCCCCCACCGGGGCGTTCGGCTGGGGGCTTCGACATCCGGTGTGCTTTTGCTTTACGGAAACAATGCaatggcgccactcgcc                                              >  2:875365/1‑224 (MQ=255)
                                                                                                                                                                                          ccGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCa                                      <  1:474660/238‑1 (MQ=255)
                                                                                                                                                                                           cGCTGGTGATGATGATGTCAGCGGCGGTGACCCCCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGc                                       <  1:67149/236‑1 (MQ=255)
                                                                                                                                                                                                gTGATGATAATGTCATCGGCGGTGACCCCCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCCCCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCa                                <  2:795671/238‑1 (MQ=255)
                                                                                                                                                                                                 tgatgagaaTGTCATCGGCGGTGACCCCCGGCCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAg                               <  2:805923/238‑1 (MQ=255)
                                                                                                                                                                                                   atgatgatgTCCTCGGCGGTGACCCCCGAGCCGTTGTCGGGCATCAGGCGGGCAATCTGCTCGCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGTTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCa                             <  1:303067/238‑1 (MQ=255)
                                                                                                                                                                                                       tgatgTAATCGGCGGTGACCCCCGAGCCGCTGTCGAGCATAAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCCCCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGa                        <  2:1100266/239‑1 (MQ=255)
                                                                                                                                                                                                            gcttCGGCGGTGACCTCCGAGCCGTGGTGGAGCAACAGGCGGGCAATCTGCTCTCGCAATACTCGCTAACCGGCTAACAAGTCACAACCGAAAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCa                   <  1:412499/236‑1 (MQ=255)
                                                                                                                                                                                                            ccacgggcggTGCCCACGGACCCGCGGCCGACCACCAGGCGGGAAATCTGCTCCCGCAATCCGCGCTGACCGGCTAAAAAGTCATACCCGAGAACGGTTTGCAGATTATGCTGCACCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATccc                    <  2:1198226/233‑1 (MQ=255)
                                                                                                                                                                                                             cagcggcggTGACCCCCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCAc                  <  1:760871/236‑1 (MQ=255)
                                                                                                                                                                                                             cATCGGCGTTGCCCCCCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCa                   <  1:109078/238‑1 (MQ=255)
                                                                                                                                                                                                               tCGGCGGTGACCACCCAGCCCCTGTCGAGCATCAGGCGGGCAATCTGCTCTCCCAATACTCGCTGACCGGCTAACAAATCATAAACCAGAACGGGTTGCCGATTATGCTGCACCACCCGGCTTTGCTCCCCCCACAGCGGTTTCCGGCTCGGGGCCTCGACATCCCGTGTGCGTTTTCTTTACGGAACAATAGAACAGCCCCTATGCGCCCCCCGCATATCCACCACCTGATCCCCCtc                >  2:976393/1‑238 (MQ=255)
                                                                                                                                                                                                                 ggcggTGACCCCCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTgg               <  1:340887/238‑1 (MQ=255)
                                                                                                                                                                                                                      tgcccaccGGGCCGCTGGCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAAt          <  2:905308/235‑1 (MQ=255)
                                                                                                                                                                                                                            cgcgcccGCGGTCGGGCATGAGGCGGCCAACCTGCTTCCGAAATCCCCGCTGCCCGGTTAACAAGTCAAAGCCGAGAACGGTTTGCATATTATGCTGCCCCACCCGGCTTAGCTCCCGCCACAGCGGTTTCGGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCAcc   <  2:788940/234‑1 (MQ=255)
                                                                                                                                                                                                                            ccGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCCCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATCCTGCACCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCAcc   <  1:990772/239‑1 (MQ=255)
                                                                                                                                                                                                                             cGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCg  >  2:1061488/1‑239 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                      
TGTTCCAGCGCCAGTTCCAGCGCTTCAACGCTGATGCCAGTTTCTGGATCGGTTGGGATTTCAATCACTTTCACGCCCATGCCGCGCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCG  >  NC_000913/4570863‑4571322

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