Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 794,939 T→C 5.2% intergenic (+17/‑150) 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,9390TC 5.2% 298.6 / 3.7 95intergenic (+17/‑150)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/43);  new base C (2/3);  total (49/46)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

AGTCTGGTATTCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTAC  >  NC_000913/794791‑795086
                                                                                                                                                    |                                                                                                                                                   
aggctcgTATTCGCCGTAATCATGGTACCTGTCGTGATGGCCAGCATCCTGGGGCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGCATGCCCGCt                                                                                                                                                     <  6:141304/143‑1 (MQ=255)
    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                                                                                                                                 <  7:269089/148‑1 (MQ=255)
                     aTGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCAttttcttttt                                                                                                                                 <  6:201130/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                                                                                           >  2:439601/1‑126 (MQ=255)
                                                                                 aaCATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTAt                                                                                           <  1:439601/126‑1 (MQ=255)
                                                                                 aaCATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAAt                                                                     >  8:288488/1‑148 (MQ=255)
                                                                                 aaCATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAAt                                                                     >  2:202952/1‑148 (MQ=255)
                                                                                 aaCATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAAt                                                                     >  3:48532/1‑148 (MQ=255)
                                                                                 aaCATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAAt                                                                     >  8:50621/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:202952/149‑1 (MQ=255)
                                                                                                    gTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGa                                                 <  6:180003/149‑1 (MQ=255)
                                                                                                    gTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGa                                                 <  1:306763/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                                                                                                         <  7:504562/82‑1 (MQ=255)
                                                                                                               cAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCg                                                                                                         >  8:504562/1‑82 (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                               >  8:4266/1‑141 (MQ=255)
                                                                                                                              cATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTc                               <  7:4266/141‑1 (MQ=255)
                                                                                                                                ttGATTCCCTGAATGCCCGCTTATTCGGGCATTTTCTATTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCATATCAACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATaa                     <  5:94426/149‑1 (MQ=255)
                                                                                                                                     tCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTa                >  3:291618/1‑149 (MQ=255)
                                                                                                                                       ccTGAATGCTCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTAc               <  4:394072/148‑1 (MQ=255)
                                                                                                                                       ccTGAATGCTCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTAc               >  3:394072/1‑148 (MQ=255)
                                                                                                                                        cTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACAt             <  6:431691/149‑1 (MQ=255)
                                                                                                                                           aaTGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACAtt            >  2:424652/1‑147 (MQ=255)
                                                                                                                                           aaTGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACAtt            >  2:451236/1‑147 (MQ=255)
                                                                                                                                           aaTGCCCGACTAAGCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTgg                                                            >  3:28064/1‑99 (MQ=255)
                                                                                                                                           aaTGCCCGACTAAGCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTgg                                                            <  4:28064/99‑1 (MQ=255)
                                                                                                                                           aaTGCCCGACTAAGCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTa                          <  8:389738/133‑1 (MQ=255)
                                                                                                                                           aaTGCCCGACTAAGCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTa                          >  7:389738/1‑133 (MQ=255)
                                                                                                                                               cccGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAgggg      >  7:279737/1‑149 (MQ=255)
                                                                                                                                               cccGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAgggg      <  3:316226/149‑1 (MQ=255)
                                                                                                                                               cccGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAgggg      <  5:269878/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)
                                                                                                                                                    |                                                                                                                                                   
AGTCTGGTATTCGCCGTAATCATGGTACCTGTCGTGATGGCCATCATCCTGGGTCTGATTTACGGTCTTGGTGAAGTATTCAACATCTTTTCTGGTGTTGGTAAAAAAGACCAGCCCGGACAAAATCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGTTAGTCGTTATATTGTCGCCTACATAACGTTACATTAAGGGGTTAC  >  NC_000913/794791‑795086

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: