Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 870,888 T→G 6.9% L493R (CTG→CGG)  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,8880TG 6.9% 314.4 / 11.4 87L493R (CTG→CGG) gsiBglutathione ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (41/40);  new base G (3/3);  total (44/43)
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.89e-01

ATGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTG  >  NC_000913/870744‑871032
                                                                                                                                                |                                                                                                                                                
aTGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGAc                                                                                                                                               >  1:239232/1‑148 (MQ=255)
   aCTTCCTGGCTCCGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGACACACAGTAAAAACCTg                                                                                                                                                 >  4:942757/1‑143 (MQ=255)
   aCTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTg                                                                                                                                                 <  3:942757/143‑1 (MQ=255)
                       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                                                                                                                                          >  3:923035/1‑119 (MQ=255)
                                  tCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGtt                                                                                                                                          <  2:6444/119‑1 (MQ=255)
                                  tCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGtt                                                                                                                                          <  4:923035/119‑1 (MQ=255)
                                  tCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGtt                                                                                                                                          >  1:6444/1‑119 (MQ=255)
                                    cggcggAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCtt                                                                                                           <  2:230146/148‑1 (MQ=255)
                                              gACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGa                                                                                                                                                >  4:175623/1‑101 (MQ=255)
                                              gACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGa                                                                                                                                                <  3:175623/101‑1 (MQ=255)
                                                             ggcggcGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGtttt                                                                                                                                        >  2:67334/1‑94 (MQ=255)
                                                             ggcggcGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGtttt                                                                                                                                        <  1:67334/94‑1 (MQ=255)
                                                              gcggcgCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCt                                                                                                            >  3:666420/1‑121 (MQ=255)
                                                              gcggcgCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCt                                                                                                            <  4:666420/121‑1 (MQ=255)
                                                                  cgcAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCa                                                                             >  2:843059/1‑148 (MQ=255)
                                                                        atatCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAg                                                                            <  1:103213/143‑1 (MQ=255)
                                                                        atatCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAg                                                                            >  2:103213/1‑143 (MQ=255)
                                                                        atatCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGt                                                                       <  3:278380/148‑1 (MQ=255)
                                                                         tatCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTg                                                                      <  2:239232/148‑1 (MQ=255)
                                                                            catcTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGaa                                                                   <  1:843059/148‑1 (MQ=255)
                                                                              tcTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTg                                                                                                                                                 <  1:1100068/68‑1 (MQ=255)
                                                                              tcTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTg                                                                                                                                                 >  2:1100068/1‑68 (MQ=255)
                                                                              tcTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATg                                                                 >  2:501068/1‑148 (MQ=255)
                                                                                  gCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGa                                                                                                                                    >  1:998398/1‑77 (MQ=255)
                                                                                  gCAAGAATCGCAGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACCGTAAAAACCTGACCGGTTTTTGGa                                                                                                                                    <  2:998398/77‑1 (MQ=255)
                                                                                   cAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTgg                                                                                                                                     >  2:1080407/1‑75 (MQ=255)
                                                                                   cAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTgg                                                                                                                                     <  1:1080407/75‑1 (MQ=255)
                                                                                   cAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAggg                                                                          <  2:348468/134‑1 (MQ=255)
                                                                                   cAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAggg                                                                          >  1:348468/1‑134 (MQ=255)
                                                                                     aGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGttttt                                                                                                                                       <  4:660180/71‑1 (MQ=255)
                                                                                     aGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGttttt                                                                                                                                       >  3:660180/1‑71 (MQ=255)
                                                                                     aGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGAttt                                                                                            <  3:1128818/114‑1 (MQ=255)
                                                                                     aGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGAttt                                                                                            >  4:1128818/1‑114 (MQ=255)
                                                                                     aGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTgg                                                                     <  4:798068/137‑1 (MQ=255)
                                                                                     aGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTgg                                                                     >  3:798068/1‑137 (MQ=255)
                                                                                       aaTCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGa                                                                    <  2:678435/136‑1 (MQ=255)
                                                                                       aaTCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGa                                                                    >  1:678435/1‑136 (MQ=255)
                                                                                        aTCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAAt                                                           >  1:1097278/1‑144 (MQ=255)
                                                                                        aTCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAAt                                                           <  2:1097278/144‑1 (MQ=255)
                                                                                         tCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAg                                                                                    <  2:658154/118‑1 (MQ=255)
                                                                                         tCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAg                                                                                    >  1:658154/1‑118 (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                                                                                                                               <  3:365194/51‑1 (MQ=255)
                                                                                                                 gAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATg                                                                                                                               >  4:365194/1‑51 (MQ=255)
                                                                                                                 gAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAg                                                                                                      <  1:262199/76‑1 (MQ=255)
                                                                                                                 gAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAg                                                                                                      >  2:262199/1‑76 (MQ=255)
                                                                                                                  aaaaaCTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGAtt                             >  3:554323/1‑148 (MQ=255)
                                                                                                                     aaCTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCtt                                                                                                           <  3:1041711/67‑1 (MQ=255)
                                                                                                                     aaCTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCtt                                                                                                           >  4:1041711/1‑67 (MQ=255)
                                                                                                                     aaCTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATg                                                                 <  1:816411/109‑1 (MQ=255)
                                                                                                                     aaCTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATg                                                                 >  2:816411/1‑109 (MQ=255)
                                                                                                                         ggTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGc                      >  1:135800/1‑148 (MQ=255)
                                                                                                                           tgtCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCa                                                                                         >  1:102887/1‑79 (MQ=255)
                                                                                                                           tgtCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCa                                                                                         <  2:102887/79‑1 (MQ=255)
                                                                                                                            gtCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATg                                                                 >  4:1141260/1‑102 (MQ=255)
                                                                                                                            gtCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATg                                                                 <  3:1141260/102‑1 (MQ=255)
                                                                                                                            gtCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGt                   <  4:554323/148‑1 (MQ=255)
                                                                                                                                     atccaaaaaCCGGTCAGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGa                                                                    <  3:853533/86‑1 (MQ=255)
                                                                                                                                     atccaaaaaCCGGTCAGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGa                                                                    >  4:853533/5‑90 (MQ=255)
                                                                                                                                       gTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACg                                                       <  4:190405/101‑1 (MQ=255)
                                                                                                                                       gTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACg                                                       >  3:190405/1‑101 (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                                                         >  1:1006304/1‑95 (MQ=255)
                                                                                                                                           aaaCCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTa                                                         <  2:1006304/95‑1 (MQ=255)
                                                                                                                                             aCCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTCGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTg  >  1:603689/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
ATGACTTCCTGGCTCAGGCACTGAAAACTAATGATCCGGCGGAAAAGACCCGCTTATATAAGGCGGCGCAGGATATCATCTGGCAAGAATCGCCGTGGATCCCGCTGGTGGTAGAAAAACTGGTGTCGGCACACAGTAAAAACCTGACCGGTTTTTGGATCATGCCAGACACCGGCTTCAGCTTTGAAGACGCGGATTTGCAATAAGCAACGCAGGGAGTGGAATGCTTAATTACGTTATCAAACGCTTACTGGGGTTGATTCCGACGCTGTTTATCGTCTCGGTGCTG  >  NC_000913/870744‑871032

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