Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 870,892 C→A 7.1% T494T (ACC→ACA) ‡ gsiB → glutathione ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913870,8920CA 7.1% 295.7 / 8.9 84T494T (ACC→ACA) ‡gsiBglutathione ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (39/39);  new base A (3/3);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTT  >  NC_000913/870767‑871038
                                                                                                                             |                                                                                                                                                  
aaaaCTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGTCCGGTTTTTgg                                                                                                                                           >  2:737282/1‑135 (MQ=255)
aaaaCTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTgg                                                                                                                                           <  1:737282/135‑1 (MQ=255)
    cTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCgg                                                                                                                          <  3:624469/148‑1 (MQ=255)
     tAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTgg                                                                                                                                           >  1:207019/1‑130 (MQ=255)
     tAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTgg                                                                                                                                           <  2:207019/130‑1 (MQ=255)
     tAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGc                                                                                                                         >  1:230146/1‑148 (MQ=255)
       aTGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCtt                                                                                                                       >  4:278380/1‑148 (MQ=255)
           tCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGtt                                                                                                                                                <  4:923035/119‑1 (MQ=255)
           tCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGtt                                                                                                                                                >  3:923035/1‑119 (MQ=255)
           tCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGtt                                                                                                                                                <  2:6444/119‑1 (MQ=255)
           tCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGtt                                                                                                                                                >  1:6444/1‑119 (MQ=255)
             cggcggAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCtt                                                                                                                 <  2:230146/148‑1 (MQ=255)
                                      ggcggcGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGtttt                                                                                                                                              >  2:67334/1‑94 (MQ=255)
                                      ggcggcGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGtttt                                                                                                                                              <  1:67334/94‑1 (MQ=255)
                                       gcggcgCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCt                                                                                                                  <  4:666420/121‑1 (MQ=255)
                                       gcggcgCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCt                                                                                                                  >  3:666420/1‑121 (MQ=255)
                                           cgcAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCa                                                                                   >  2:843059/1‑148 (MQ=255)
                                                 atatCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAg                                                                                  >  2:103213/1‑143 (MQ=255)
                                                 atatCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAg                                                                                  <  1:103213/143‑1 (MQ=255)
                                                 atatCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGt                                                                             <  3:278380/148‑1 (MQ=255)
                                                  tatCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTg                                                                            <  2:239232/148‑1 (MQ=255)
                                                     catcTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGaa                                                                         <  1:843059/148‑1 (MQ=255)
                                                       tcTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATg                                                                       >  2:501068/1‑148 (MQ=255)
                                                           gCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGa                                                                                                                                          >  1:998398/1‑77 (MQ=255)
                                                           gCAAGAATCGCAGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACCGTAAAAACCTGACCGGTTTTTGGa                                                                                                                                          <  2:998398/77‑1 (MQ=255)
                                                            cAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTgg                                                                                                                                           <  1:1080407/75‑1 (MQ=255)
                                                            cAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTgg                                                                                                                                           >  2:1080407/1‑75 (MQ=255)
                                                            cAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAggg                                                                                <  2:348468/134‑1 (MQ=255)
                                                            cAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAggg                                                                                >  1:348468/1‑134 (MQ=255)
                                                              aGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGttttt                                                                                                                                             >  3:660180/1‑71 (MQ=255)
                                                              aGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGttttt                                                                                                                                             <  4:660180/71‑1 (MQ=255)
                                                              aGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGAttt                                                                                                  >  4:1128818/1‑114 (MQ=255)
                                                              aGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGAttt                                                                                                  <  3:1128818/114‑1 (MQ=255)
                                                              aGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTgg                                                                           <  4:798068/137‑1 (MQ=255)
                                                              aGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTgg                                                                           >  3:798068/1‑137 (MQ=255)
                                                                aaTCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGa                                                                          >  1:678435/1‑136 (MQ=255)
                                                                aaTCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGa                                                                          <  2:678435/136‑1 (MQ=255)
                                                                 aTCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAAt                                                                 <  2:1097278/144‑1 (MQ=255)
                                                                 aTCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAAt                                                                 >  1:1097278/1‑144 (MQ=255)
                                                                  tCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAg                                                                                          >  1:658154/1‑118 (MQ=255)
                                                                  tCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAg                                                                                          <  2:658154/118‑1 (MQ=255)
                                                                     ccGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTAt                                                         <  1:501068/148‑1 (MQ=255)
                                                                         ggATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCaaa                                                     >  4:561725/1‑148 (MQ=255)
                                                                           aTCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGcc                                                                                                                                   >  1:56989/1‑68 (MQ=255)
                                                                           aTCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGcc                                                                                                                                   <  2:56989/68‑1 (MQ=255)
                                                                                 cTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGt                                                            >  3:471455/1‑133 (MQ=255)
                                                                                 cTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGt                                                            <  4:471455/133‑1 (MQ=255)
                                                                                    gtggtAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTggg                                          <  3:561725/148‑1 (MQ=255)
                                                                                      ggtAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGt                                                                             >  2:364238/1‑111 (MQ=255)
                                                                                      ggtAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGt                                                                             <  1:364238/111‑1 (MQ=255)
                                                                                        tAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAAt                                                                                             >  1:683215/1‑93 (MQ=255)
                                                                                        tAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAAt                                                                                             <  2:683215/93‑1 (MQ=255)
                                                                                         aGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGa                                     >  3:77957/1‑148 (MQ=255)
                                                                                          gAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATg                                                                                                                                     >  4:365194/1‑51 (MQ=255)
                                                                                          gAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATg                                                                                                                                     <  3:365194/51‑1 (MQ=255)
                                                                                          gAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAg                                                                                                            >  2:262199/1‑76 (MQ=255)
                                                                                          gAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAg                                                                                                            <  1:262199/76‑1 (MQ=255)
                                                                                           aaaaaCTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGAtt                                   >  3:554323/1‑148 (MQ=255)
                                                                                              aaCTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCtt                                                                                                                 <  3:1041711/67‑1 (MQ=255)
                                                                                              aaCTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCtt                                                                                                                 >  4:1041711/1‑67 (MQ=255)
                                                                                              aaCTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATg                                                                       >  2:816411/1‑109 (MQ=255)
                                                                                              aaCTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATg                                                                       <  1:816411/109‑1 (MQ=255)
                                                                                                  ggTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGc                            >  1:135800/1‑148 (MQ=255)
                                                                                                    tgtCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCa                                                                                               >  1:102887/1‑79 (MQ=255)
                                                                                                    tgtCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCa                                                                                               <  2:102887/79‑1 (MQ=255)
                                                                                                     gtCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATg                                                                       <  3:1141260/102‑1 (MQ=255)
                                                                                                     gtCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATg                                                                       >  4:1141260/1‑102 (MQ=255)
                                                                                                     gtCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGt                         <  4:554323/148‑1 (MQ=255)
                                                                                                              atccaaaaaCCGGTCAGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGa                                                                          <  3:853533/86‑1 (MQ=255)
                                                                                                              atccaaaaaCCGGTCAGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGa                                                                          >  4:853533/5‑90 (MQ=255)
                                                                                                                gTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACg                                                             >  3:190405/1‑101 (MQ=255)
                                                                                                                gTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACg                                                             <  4:190405/101‑1 (MQ=255)
                                                                                                                gTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCg              <  4:77957/148‑1 (MQ=255)
                                                                                                                  aaaaaCCGGTCAGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTa                                                                   <  4:1135331/93‑1 (MQ=255)
                                                                                                                  aaaaaCCGGTCAGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTa                                                                   >  3:1135331/1‑93 (MQ=255)
                                                                                                                  aaaaaCCGGTCAGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCg                                <  1:230356/128‑1 (MQ=255)
                                                                                                                  aaaaaCCGGTCAGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCg                                >  2:230356/1‑128 (MQ=255)
                                                                                                                    aaaCCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTa                                                               <  2:1006304/95‑1 (MQ=255)
                                                                                                                    aaaCCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTa                                                               >  1:1006304/1‑95 (MQ=255)
                                                                                                                      aCCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTCGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTg        >  1:603689/1‑148 (MQ=255)
                                                                                                                          gACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAg                                                                              >  3:432955/1‑74 (MQ=255)
                                                                                                                          gACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAg                                                                              <  4:432955/74‑1 (MQ=255)
                                                                                                                           aCCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGAcg                                                                                                         <  4:218247/46‑1 (MQ=255)
                                                                                                                           aCCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGAcg                                                                                                         >  3:218247/1‑46 (MQ=255)
                                                                                                                            ccGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATaa                                                                                           >  4:241380/1‑59 (MQ=255)
                                                                                                                            ccGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATaa                                                                                           <  3:241380/59‑1 (MQ=255)
                                                                                                                            ccGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCt         >  1:90505/1‑141 (MQ=255)
                                                                                                                            ccGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCt         <  2:90505/141‑1 (MQ=255)
                                                                                                                            ccGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTCGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGttt  <  2:603689/148‑1 (MQ=255)
                                                                                                                             cGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCtt                                                >  4:6151/1‑101 (MQ=255)
                                                                                                                             cGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCtt                                                <  3:6151/101‑1 (MQ=255)
                                                                                                                             |                                                                                                                                                  
AAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTGGTGTTT  >  NC_000913/870767‑871038

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