Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 928,458 C→T 10.2% L246L (TTG→TTA cydC ← glutathione/L‑cysteine ABC exporter subunit CydC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913928,4580CT10.2% 280.6 / 15.8 98L246L (TTG→TTAcydCglutathione/L‑cysteine ABC exporter subunit CydC
Reads supporting (aligned to +/- strand):  ref base C (45/43);  new base T (6/4);  total (51/47)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

GGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAGTTGTTGGCGATACTG  >  NC_000913/928309‑928604
                                                                                                                                                     |                                                                                                                                                  
ggTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGAc                                                                                                                                                    <  1:310737/150‑1 (MQ=255)
ggTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGAc                                                                                                                                                    <  1:1042870/150‑1 (MQ=255)
 gTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCGTTATCGCTTGCGACa                                                                                                                                                   >  1:1371806/1‑150 (MQ=255)
   gCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAAt                                                                                                                                                 <  1:929569/150‑1 (MQ=255)
     cAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGc                                                                                                                                               >  1:1143878/1‑150 (MQ=255)
     cAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGc                                                                                                                                               >  1:1149197/1‑150 (MQ=255)
     cAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGc                                                                                                                                               <  1:1364709/150‑1 (MQ=255)
       gTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTGTCGCAAGCGATAATGctg                                                                                                                                             <  2:1443434/150‑3 (MQ=255)
        tGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTGTCGCAAGCGATAATGctgc                                                                                                                                            <  1:1298812/150‑4 (MQ=255)
        tGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACCCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTGTCGCAAGCGATAATGctgc                                                                                                                                            <  1:1261497/150‑4 (MQ=255)
               aaCGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCa                                                                                                                                     >  2:1015722/1‑150 (MQ=255)
               aaCGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCa                                                                                                                                     >  1:1447459/1‑150 (MQ=255)
                     gCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGa                                                                                                                               <  1:1522912/150‑1 (MQ=255)
                      cTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGAc                                                                                                                              >  1:621068/1‑150 (MQ=255)
                      cTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGAc                                                                                                                              >  1:72153/1‑150 (MQ=255)
                       tAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACg                                                                                                                             >  1:862594/1‑150 (MQ=255)
                        aaCGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACgg                                                                                                                            <  2:91351/150‑1 (MQ=255)
                        aaCGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACgg                                                                                                                            >  1:1275554/1‑150 (MQ=255)
                            cgcAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGtt                                                                                                                        <  2:190775/150‑1 (MQ=255)
                            cgcAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGtt                                                                                                                        >  1:513814/1‑150 (MQ=255)
                                     cAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCa                                                                                                               <  1:1129037/150‑1 (MQ=255)
                                        aCAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGcc                                                                                                            >  1:646764/1‑150 (MQ=255)
                                              cAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGaa                                                                                                      <  2:900067/150‑1 (MQ=255)
                                               aaTTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAAt                                                                                                     >  2:373703/1‑150 (MQ=255)
                                                aTTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAAtt                                                                                                    <  1:1243402/150‑1 (MQ=255)
                                                  tAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTc                                                                                                  <  1:1015722/150‑1 (MQ=255)
                                                  tAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTc                                                                                                  <  2:1447459/150‑1 (MQ=255)
                                                    aCGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCGTTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTg                                                                                                <  2:1371806/150‑1 (MQ=255)
                                                     cgcgCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTCCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGt                                                                                               >  1:1284215/1‑150 (MQ=255)
                                                      gcgcCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTa                                                                                              <  2:1185904/150‑1 (MQ=255)
                                                       cgcCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTAt                                                                                             >  1:6350/1‑150 (MQ=255)
                                                          cGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATtct                                                                                          <  2:1275554/150‑1 (MQ=255)
                                                           gggTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGc                                                                                                                           <  2:131294/116‑1 (MQ=255)
                                                           gggTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGc                                                                                                                           >  1:131294/1‑116 (MQ=255)
                                                            ggTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTctct                                                                                        <  1:1505948/150‑1 (MQ=255)
                                                            ggTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTctct                                                                                        <  1:1217466/150‑1 (MQ=255)
                                                                gAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGTGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGtt                                                                                    >  1:453980/1‑150 (MQ=255)
                                                                 aGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTtg                                                                                   <  1:396907/150‑1 (MQ=255)
                                                                 aGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTtg                                                                                   <  2:130526/150‑1 (MQ=255)
                                                                 aGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTCCCAGCCATTGAATTTCTGTATTCTCTAGTtg                                                                                   <  1:524730/150‑1 (MQ=255)
                                                                     ttGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTtgcgtg                                                                               >  1:1415709/1‑150 (MQ=255)
                                                                        ccgccAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGa                                                                            >  2:820820/1‑150 (MQ=255)
                                                                                 ccgccAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCg                                                                   <  1:912901/150‑1 (MQ=255)
                                                                                   gccAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCt                                                                 >  1:1402361/1‑150 (MQ=255)
                                                                                   gccAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCt                                                                 >  1:1459774/1‑150 (MQ=255)
                                                                                     cAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCATATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTgg                                                               <  2:511522/150‑1 (MQ=255)
                                                                                      aGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGTGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGc                                                              <  2:453980/150‑1 (MQ=255)
                                                                                       gACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCa                                                             >  1:1056897/1‑150 (MQ=255)
                                                                                       gACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCa                                                             <  2:6350/150‑1 (MQ=255)
                                                                                       gACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGATTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGTTCGCTGTCa                                                             <  2:1499237/150‑1 (MQ=255)
                                                                                          gCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCa                                                          >  1:1095052/1‑150 (MQ=255)
                                                                                          gCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCa                                                          >  2:1035156/1‑150 (MQ=255)
                                                                                               ccACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAAt                                                     >  2:14576/1‑150 (MQ=255)
                                                                                               ccACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAAt                                                     <  2:1262351/150‑1 (MQ=255)
                                                                                               ccACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCATGCGACAAAGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAAt                                                     >  2:779047/1‑150 (MQ=255)
                                                                                               ccACAGCATAAGGATCACCGCTAACGCGCCAATCAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAAt                                                     >  2:335831/1‑150 (MQ=255)
                                                                                                cacaGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGCTCGCTGGCACCAAAAATg                                                    >  2:879270/1‑150 (MQ=255)
                                                                                                   aGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCTTTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGtc                                                 <  1:1075060/150‑1 (MQ=255)
                                                                                                    gCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGtca                                                <  1:694242/150‑1 (MQ=255)
                                                                                                    gCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGtca                                                >  1:462019/1‑150 (MQ=255)
                                                                                                      aTCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGTtcagc                                              <  2:646764/150‑1 (MQ=255)
                                                                                                            aTCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTcagctcagct                                        <  2:513814/150‑1 (MQ=255)
                                                                                                              cACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTg                                      <  2:462019/150‑1 (MQ=255)
                                                                                                                ccGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGcc                                    >  1:1455502/1‑150 (MQ=255)
                                                                                                                    tAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTg                                <  2:6850/150‑1 (MQ=255)
                                                                                                                         cgcCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGcc                           >  1:996552/1‑150 (MQ=255)
                                                                                                                         cgcCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGcc                           >  1:96985/1‑150 (MQ=255)
                                                                                                                            cAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAgg                        <  2:1095052/150‑1 (MQ=255)
                                                                                                                             aaTGACCAGGATTATCGCTTGCGTCCATGCGGTCATTTCAGAATGACGGCGGTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGGCAGCTCAGCTTGCCCTTGCAGCCAGGc                       <  1:1409281/150‑1 (MQ=255)
                                                                                                                              aTGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGcc                      >  2:1299462/1‑150 (MQ=255)
                                                                                                                               tGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGccn                     <  1:73669/150‑2 (MQ=255)
                                                                                                                                gAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCa                              <  2:697220/140‑1 (MQ=255)
                                                                                                                                gAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCa                              >  1:697220/1‑140 (MQ=255)
                                                                                                                                 agcagcATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTc                   >  2:1014026/1‑150 (MQ=255)
                                                                                                                                   ccgcaTTGTCGCAAGCGATAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAg                 >  2:47486/3‑150 (MQ=255)
                                                                                                                                   ccgcaTTGTCGCAAGCGATAATGCGGTCAGTTCAGATTGAAGGAGTTGCGATTCCAGCAATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAACAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGGCGTCAg                 >  2:968275/3‑150 (MQ=255)
                                                                                                                                   cagcaTTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAg                 <  1:335831/150‑1 (MQ=255)
                                                                                                                                   cagcaTTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAg                 <  1:100634/150‑1 (MQ=255)
                                                                                                                                   cagcaTTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAg                 <  2:773138/150‑1 (MQ=255)
                                                                                                                                    cgcaTTGTCGCAAGCGATAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAgt                >  2:369050/2‑150 (MQ=255)
                                                                                                                                    agcaTTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAgt                >  2:30023/1‑150 (MQ=255)
                                                                                                                                    agcaTTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAgt                >  1:619220/1‑150 (MQ=255)
                                                                                                                                    agcaTTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAgt                <  2:1402361/150‑1 (MQ=255)
                                                                                                                                    agcaTTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAgt                >  1:1442452/1‑150 (MQ=255)
                                                                                                                                     gcaTTGTCGCAAGCGATAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTc                   <  2:1482320/146‑1 (MQ=255)
                                                                                                                                     gcaTTGTCGCAAGCGATAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTc                   >  1:1482320/1‑146 (MQ=255)
                                                                                                                                     gcaTTGTCGCAAGCGATAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAgtt               >  2:1387/1‑150 (MQ=255)
                                                                                                                                     gcaTTGTCGCAAGCGATAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAgtt               >  1:151875/1‑150 (MQ=255)
                                                                                                                                     gcaTTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAgtt               >  2:817207/1‑150 (MQ=255)
                                                                                                                                      caTTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAgttg              <  1:503177/150‑1 (MQ=255)
                                                                                                                                      caTTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAgttg              <  1:1035156/150‑1 (MQ=255)
                                                                                                                                        ttATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCGTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAgtggtt            >  2:651030/1‑150 (MQ=255)
                                                                                                                                           tCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAGTTGTTGGc         >  1:933778/1‑150 (MQ=255)
                                                                                                                                           tCGCTTGCGACAATGCGGTCAGTCCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAGTTGTTGGc         <  2:1284215/150‑1 (MQ=255)
                                                                                                                                             gCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAGTTGTTGGCGa       >  1:269028/1‑150 (MQ=255)
                                                                                                                                             gCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAGTTGTTGGCGa       >  1:1271149/1‑150 (MQ=255)
                                                                                                                                                tGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGtcagctc                                            >  1:926521/1‑110 (MQ=255)
                                                                                                                                                tGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGtcagctc                                            <  2:926521/110‑1 (MQ=255)
                                                                                                                                                 gCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAGTTGTTGGCGATACt   >  1:1060626/1‑150 (MQ=255)
                                                                                                                                                  cGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAGTTGTTGGCGATACTg  >  2:1188932/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAGTTGTTGGCGATACTG  >  NC_000913/928309‑928604

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