Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 794,943 T→G 5.7% intergenic (+21/‑146) acrZ → / → modA AcrAB‑TolC efflux pump accessory protein, membrane‑associated/molybdate ABC transporter periplasmic binding protein; chlorate resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913794,9430TG 5.7% 300.2 / 4.3 97intergenic (+21/‑146)acrZ/modAAcrAB‑TolC efflux pump accessory protein, membrane‑associated/molybdate ABC transporter periplasmic binding protein; chlorate resistance protein
Reads supporting (aligned to +/- strand):  ref base T (47/44);  new base G (2/4);  total (49/48)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTATTCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTACCAATG  >  NC_000913/794795‑795091
                                                                                                                                                    |                                                                                                                                                    
tGGTATTCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGt                                                                                                                                                      >  2:224733/1‑149 (MQ=255)
  gTATTCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCg                                                                                                                                                    >  1:130078/1‑149 (MQ=255)
   tATTCGCCTTAATCATGGTACCTTTCGTGATGTCCATCATCCTGGGTCTGATTTACGCTCTTGGTGAAGTATTCAACATCTTTTCTGTTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCgg                                                                                                                                                   <  8:184304/149‑1 (MQ=255)
     ttCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCggg                                                                                                                                                  >  3:115703/1‑148 (MQ=255)
        gCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATTTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGACTAAGCGGGCAtt                                                                                                                                              <  2:360514/149‑1 (MQ=255)
        gCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAtt                                                                                                                                              >  6:194649/1‑149 (MQ=255)
        gCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAtt                                                                                                                                              <  6:288178/149‑1 (MQ=255)
        gCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAtt                                                                                                                                              <  7:341809/149‑1 (MQ=255)
        gCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAt                                                                                                                                               <  3:54896/148‑1 (MQ=255)
              aTCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttctt                                                                                                                                         <  2:130078/148‑1 (MQ=255)
               tCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAATAAGACCAGCCCGTACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttctttt                                                                                                                                       <  7:151742/149‑1 (MQ=255)
                 aTGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGACTGCCCGCTTAGGCGGGCAttttctttttc                                                                                                                                     <  7:243243/149‑1 (MQ=255)
                 aTGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttctttttc                                                                                                                                     <  4:318416/149‑1 (MQ=255)
                 aTGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttcttttt                                                                                                                                      <  6:201130/148‑1 (MQ=255)
                 aTGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttcttttt                                                                                                                                      <  7:269089/148‑1 (MQ=255)
                      aCCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAAc                                                                                                                                >  5:180003/1‑149 (MQ=255)
                           tCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTcc                                                                                                                            <  2:364941/148‑1 (MQ=255)
                           tCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCt                                                                                                                           <  4:115703/149‑1 (MQ=255)
                             gTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGc                                                                                                                         <  5:194649/149‑1 (MQ=255)
                                     catcatCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTg                                                                                                                 >  2:252341/1‑149 (MQ=255)
                                     catcatCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTg                                                                                                                 >  6:177108/1‑149 (MQ=255)
                                      atcatcCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGc                                                                                                                <  5:177108/149‑1 (MQ=255)
                                             tGGGTCTGATTTACGGTCTGGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATAtt                                                                                                         >  2:306763/1‑149 (MQ=255)
                                                gTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATAttttt                                                                                                      >  1:216808/1‑149 (MQ=255)
                                                 tCTGATTTACGGTCTTGGTGAAGTATTAAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCTCACAAAATCATTGATTCCCTGTATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATAtttttt                                                                                                     <  5:372179/149‑1 (MQ=255)
                                                  cTGATTTACGGTCTTGGTGAAGTAATCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTc                                                                                                    <  6:72441/149‑1 (MQ=255)
                                                   tGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttctttttc                                                                                                                                     >  3:258979/1‑115 (MQ=255)
                                                   tGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttctttttc                                                                                                                                     <  4:258979/115‑1 (MQ=255)
                                                     aTTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTa                                                                                                 <  1:224733/149‑1 (MQ=255)
                                                      tttACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGAAAAAATCAATGATTCCCTGAATGCACGCTTAGTCGGGCATTTTCTTTTTCTAAACTTCCTGCTTTTAATGCCGATATTTTTTCTTa                                                                                                 >  7:249973/1‑148 (MQ=255)
                                                             tCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCt                                                                                          >  8:240490/1‑148 (MQ=255)
                                                               ttGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTca                                                                                        >  3:83736/1‑148 (MQ=255)
                                                                   tGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAg                                                                                   <  7:240490/149‑1 (MQ=255)
                                                                    gAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAgg                                                                                  >  8:334/1‑149 (MQ=255)
                                                                     aaGTATTCAACATCTCTTCTGGTGTTGGTATAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCAGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGGCGATATTTTTTCTTATCTACCTAACAAAGGt                                                                                 >  6:94426/1‑149 (MQ=255)
                                                                       gTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGt                                                                                 >  8:312611/1‑147 (MQ=255)
                                                                             aaCATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTAt                                                                                                <  1:439601/126‑1 (MQ=255)
                                                                             aaCATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTAt                                                                                                >  2:439601/1‑126 (MQ=255)
                                                                             aaCATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAAt                                                                          >  8:50621/1‑148 (MQ=255)
                                                                             aaCATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAAt                                                                          >  3:48532/1‑148 (MQ=255)
                                                                             aaCATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAAt                                                                          >  2:202952/1‑148 (MQ=255)
                                                                             aaCATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAAt                                                                          >  8:288488/1‑148 (MQ=255)
                                                                              aCATCTATACTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATAGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTATTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATa                                                                         >  8:335073/1‑148 (MQ=255)
                                                                                aTCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTCCTTATCTACCTCACAAAGGTTAGCAATAAct                                                                      >  5:124491/1‑149 (MQ=255)
                                                                                aTCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGACCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTCCTTTTCCTGCCGATATTTTTTCTTATCTACCTAACAAAGCTTAGCAATAAct                                                                      >  7:275427/1‑149 (MQ=255)
                                                                                   ttttCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTAGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAActgct                                                                   >  1:122265/1‑149 (MQ=255)
                                                                                     ttCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAActgctg                                                                  <  7:50621/148‑1 (MQ=255)
                                                                                     ttCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTCCTTATCTACCTCACAAAGGTTAGCAATAActg                                                                     <  6:124491/145‑1 (MQ=255)
                                                                                      tCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTgg                                                                 <  1:252341/148‑1 (MQ=255)
                                                                                        tGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCCGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGa                                                               >  4:79263/1‑148 (MQ=255)
                                                                                          gtgtTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGaaaa                                                            >  4:316226/1‑149 (MQ=255)
                                                                                            gtTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAAtt                                                          >  1:29845/1‑149 (MQ=255)
                                                                                             ttGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTc                                                         >  5:185789/1‑149 (MQ=255)
                                                                                                gTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTcc                                                        >  3:212050/1‑147 (MQ=255)
                                                                                                gTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGa                                                      <  1:306763/149‑1 (MQ=255)
                                                                                                gTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGa                                                      <  6:180003/149‑1 (MQ=255)
                                                                                                gTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGa                                                      <  1:202952/149‑1 (MQ=255)
                                                                                                gTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCCGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCg                                                       <  3:79263/148‑1 (MQ=255)
                                                                                                 tAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCACATTCGTGGTTTTCGTGCCGATATTTTTTGTTATCTCCCTCACAAAGTTTAGCATTAAGTGCTGGGAAAATTCCGAg                                                     >  2:22045/1‑149 (MQ=255)
                                                                                                 tAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAg                                                     >  4:234836/1‑149 (MQ=255)
                                                                                                           cAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCg                                                                                                              >  8:504562/1‑82 (MQ=255)
                                                                                                           cAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCg                                                                                                              <  7:504562/82‑1 (MQ=255)
                                                                                                           cAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTta                                           <  7:288488/149‑1 (MQ=255)
                                                                                                            ccgtgccaCAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTtat                                          <  5:197519/142‑1 (MQ=255)
                                                                                                             gCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTGTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTtata                                         >  2:411532/1‑149 (MQ=255)
                                                                                                                 ggACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGt                                     <  7:312611/149‑1 (MQ=255)
                                                                                                                  gACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAAct                                                                      <  8:362427/115‑1 (MQ=255)
                                                                                                                  gACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAAct                                                                      >  7:362427/1‑115 (MQ=255)
                                                                                                                  gACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTc                                    <  6:209531/149‑1 (MQ=255)
                                                                                                                     aaaaTCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAg                                                     >  2:500195/1‑129 (MQ=255)
                                                                                                                     aaaaTCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAg                                                     <  1:500195/129‑1 (MQ=255)
                                                                                                                     aaaaTCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGcc                                 >  1:312786/1‑149 (MQ=255)
                                                                                                                       aaTCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTa                               >  1:217989/1‑149 (MQ=255)
                                                                                                                          cATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTc                                    <  7:4266/141‑1 (MQ=255)
                                                                                                                          cATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTc                                    >  8:4266/1‑141 (MQ=255)
                                                                                                                            ttGATTCCCTGAATGCCCGCTTATTCGGGCATTTTCTATTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCATATCAACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATaa                          <  5:94426/149‑1 (MQ=255)
                                                                                                                                 tCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTa                     >  3:291618/1‑149 (MQ=255)
                                                                                                                                   ccTGAATGCTCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTAc                    >  3:394072/1‑148 (MQ=255)
                                                                                                                                   ccTGAATGCTCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTAc                    <  4:394072/148‑1 (MQ=255)
                                                                                                                                    cTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACAt                  <  6:431691/149‑1 (MQ=255)
                                                                                                                                       aaTGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACAtt                 >  2:424652/1‑147 (MQ=255)
                                                                                                                                       aaTGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACAtt                 >  2:451236/1‑147 (MQ=255)
                                                                                                                                       aaTGCCCGACTAAGCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTgg                                                                 <  4:28064/99‑1 (MQ=255)
                                                                                                                                       aaTGCCCGACTAAGCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTgg                                                                 >  3:28064/1‑99 (MQ=255)
                                                                                                                                       aaTGCCCGACTAAGCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTa                               >  7:389738/1‑133 (MQ=255)
                                                                                                                                       aaTGCCCGACTAAGCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTa                               <  8:389738/133‑1 (MQ=255)
                                                                                                                                           cccGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAgggg           >  7:279737/1‑149 (MQ=255)
                                                                                                                                           cccGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAgggg           <  5:269878/149‑1 (MQ=255)
                                                                                                                                           cccGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAgggg           <  3:316226/149‑1 (MQ=255)
                                                                                                                                            ccGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAg                                                                                   >  6:98862/1‑76 (MQ=255)
                                                                                                                                            ccGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGt          >  2:129064/1‑149 (MQ=255)
                                                                                                                                            ccGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACATCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCg                                                       >  7:196707/1‑104 (MQ=255)
                                                                                                                                            acGCTTAGTCGTGCATTTTCTTTTTCTGAACTTCCTGCTTTTCCTGCAGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCg                                                       <  8:196707/103‑1 (MQ=255)
                                                                                                                                             cGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGtt         <  1:129064/149‑1 (MQ=255)
                                                                                                                                               cTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGGTATATTGTCGCCTACATAACGTTACATTAAGGGGTTAc       <  4:83736/149‑1 (MQ=255)
                                                                                                                                                ttAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAg                                                                                   <  1:437406/72‑1 (MQ=255)
                                                                                                                                                ttAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAg                                                                                   >  2:437406/1‑72 (MQ=255)
                                                                                                                                                   gTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTACCAAt   <  2:312786/149‑1 (MQ=255)
                                                                                                                                                    tCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTACCAATg  >  3:442790/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGGTATTCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTACCAATG  >  NC_000913/794795‑795091

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