Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,448,424 T→A 9.0% intergenic (‑179/‑23) yzfA ← / → chpS putative YjfA/ChpS antitoxin of the ChpB‑ChpS toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,448,4240TA 9.0% 321.5 / 14.8 109intergenic (‑179/‑23)yzfA/chpSputative YjfA/ChpS antitoxin of the ChpB‑ChpS toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base T (43/56);  new base A (5/5);  total (48/61)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 111 total aligned reads displayed.

CGAGACGACTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATCAACTGATTC  >  NC_000913/4448275‑4448573
                                                                                                                                                     |                                                                                                                                                     
cGAGACGACTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACtt                                                                                                                                                       >  1:1355940/1‑150 (MQ=255)
   gacgacTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGt                                                                                                                                                    <  2:483692/150‑1 (MQ=255)
    ccgacTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGTGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGta                                                                                                                                                   <  2:1124578/149‑1 (MQ=255)
    acgacTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGta                                                                                                                                                   <  1:9346/150‑1 (MQ=255)
      gacTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGtata                                                                                                                                                 <  1:879238/150‑1 (MQ=255)
          cGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGta                                                                                                                                                   >  2:397857/1‑144 (MQ=255)
           gCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTa                                                                                                                                            >  1:243590/1‑150 (MQ=255)
           gCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACAAAGTATAACtca                                                                                                                                            <  1:1125065/150‑3 (MQ=255)
            cATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACAAAGTATAACtcac                                                                                                                                           <  2:1345252/150‑4 (MQ=255)
            cATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACAAAGTATAACtcac                                                                                                                                           <  1:1443610/150‑4 (MQ=255)
             aTTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAag                                                                                                                                          >  1:1340009/1‑150 (MQ=255)
             aTTCTGATTTGTAATTCCCTCACCTTTTCCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACTCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGCGTAAAACGTTTCCTGTAGCACCGTGAGTTATACAAAGTATAACtcacg                                                                                                                                          <  2:407539/150‑5 (MQ=255)
              ttCTGTTTTGTAATTCACTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAagg                                                                                                                                         <  2:1414413/150‑1 (MQ=255)
                  gTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGt                                                                                                                                     >  1:1450970/1‑150 (MQ=255)
                    tttGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGc                                                                                                                                   >  1:315074/1‑150 (MQ=255)
                    tttGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGc                                                                                                                                   <  1:291630/150‑1 (MQ=255)
                        tAATTCCCTCACCTTTTGCTTTTCTCTCCGAGACGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATAATTTGTATAACTTAAGGAGGTGCAGAt                                                                                                                               >  2:533667/1‑150 (MQ=255)
                         aaTTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATg                                                                                                                              <  2:1519635/150‑1 (MQ=255)
                           ttCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTTCAGATGCg                                                                                                                            <  2:30386/150‑1 (MQ=255)
                             cccTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTa                                                                                                                          >  1:106342/1‑150 (MQ=255)
                              ccTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTAt                                                                                                                         >  1:803785/1‑150 (MQ=255)
                                   ccTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGc                                                                                                                             >  2:103432/1‑141 (MQ=255)
                                   ccTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGc                                                                                                                             <  1:103432/141‑1 (MQ=255)
                                   ccTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCa                                                                                                                    <  2:1138132/150‑1 (MQ=255)
                                     ttttGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATa                                                                                                                  >  2:291709/1‑150 (MQ=255)
                                     ttttGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATa                                                                                                                  >  1:465572/1‑150 (MQ=255)
                                      tttGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATaa                                                                                                                 >  1:182877/1‑150 (MQ=255)
                                         gCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATaaaaa                                                                                                              >  1:174462/1‑150 (MQ=255)
                                         gCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCATGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATaaaaa                                                                                                              <  1:642359/150‑1 (MQ=255)
                                              tctctcCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACGGTGTGTTATACTTTGTATATCTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATgg                                                                                                         <  2:287652/150‑1 (MQ=255)
                                               ctctcCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATggg                                                                                                        <  1:915829/150‑1 (MQ=255)
                                                 ctcCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTTAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATggggg                                                                                                      <  2:1447285/150‑1 (MQ=255)
                                                 ctcCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATggggg                                                                                                      <  2:283259/150‑1 (MQ=255)
                                                 ctcCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATggggg                                                                                                      >  1:983695/1‑150 (MQ=255)
                                                 ctcCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTACTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATggggg                                                                                                      >  2:64587/1‑150 (MQ=255)
                                                   ccGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGaa                                                                                                    <  2:1382801/150‑1 (MQ=255)
                                                          gCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCa                                                                                             <  1:1301456/150‑1 (MQ=255)
                                                            tttCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAgg                                                                                           <  2:1538682/150‑1 (MQ=255)
                                                              tCCATATCTATTAACGCATAAAAAACTCTGGTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTa                                                                                         <  1:1208272/150‑1 (MQ=255)
                                                              tCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGtg                                                                                         <  1:1073239/150‑2 (MQ=255)
                                                              tCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTa                                                                                         <  2:931058/150‑1 (MQ=255)
                                                                cATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATg                                                                                       >  1:832979/1‑150 (MQ=255)
                                                                 atatCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATgg                                                                                      >  2:398471/1‑150 (MQ=255)
                                                                        ttAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTcc                                                                               >  1:1180962/1‑150 (MQ=255)
                                                                          aaCGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCa                                                                             >  1:311847/1‑150 (MQ=255)
                                                                          aaCGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCa                                                                             >  1:306351/1‑150 (MQ=255)
                                                                              cATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAatat                                                                         >  1:1121026/1‑150 (MQ=255)
                                                                               aaaaaaaaATCTGAAGGCAAACACAAAAGCGCAGGGGTAAAAACTTTCCTCTAGAACACAGAAATATAAATTGTATAACTTAAGTAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGAAGGAATGGTCATTCCCAATATc                                                                        >  2:21141/3‑150 (MQ=255)
                                                                               aTAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGCATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATc                                                                        >  1:404513/1‑150 (MQ=255)
                                                                                tAAAAAACTCTGCTGGCATTCATAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCg                                                                       <  2:200686/150‑1 (MQ=255)
                                                                                tAAAAAACTCTGCTGGCATTCACAAATGCGCATGGGTAAAACGTTTCCTGTACCACCGTGAGTTATACTTTGTATAACTTAAGGATGTGCATATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCg                                                                       <  1:714058/150‑1 (MQ=255)
                                                                                tAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCg                                                                       >  1:1505542/1‑150 (MQ=255)
                                                                                  aaaaaCTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTa                                                                     <  1:349440/150‑1 (MQ=255)
                                                                                  aaaaaCTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTa                                                                     >  2:508227/1‑150 (MQ=255)
                                                                                    aaaCTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAAt                                                                   <  1:794408/150‑1 (MQ=255)
                                                                                     aaCTCTGCTGGCATTCACAAATGCGCATGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATg                                                                  <  1:533667/150‑1 (MQ=255)
                                                                                     aaCTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATg                                                                  <  1:508227/150‑1 (MQ=255)
                                                                                      aCTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGa                                                                 <  2:383644/150‑1 (MQ=255)
                                                                                       ctctGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGaa                                                                >  2:1147872/1‑150 (MQ=255)
                                                                                         ctgctgGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAg                                                              <  2:315074/150‑1 (MQ=255)
                                                                                           gctgGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGaa                                                            >  2:1032178/1‑150 (MQ=255)
                                                                                           gctgGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGaa                                                            >  1:1425186/1‑150 (MQ=255)
                                                                                            ctgGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGaac                                                           <  1:291709/150‑1 (MQ=255)
                                                                                               gCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGaactta                                                        <  2:1340892/150‑1 (MQ=255)
                                                                                                 aTTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGaacttaac                                                      <  2:1450970/150‑1 (MQ=255)
                                                                                                   tCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGaacttaactt                                                    <  2:263064/150‑1 (MQ=255)
                                                                                                    cacaAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCAt                                                                                                                   >  1:541573/1‑86 (MQ=255)
                                                                                                    cacaAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCAt                                                                                                                   <  2:541573/86‑1 (MQ=255)
                                                                                                      caAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACa                                                 <  2:465572/150‑1 (MQ=255)
                                                                                                      caAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACa                                                 >  2:969825/1‑150 (MQ=255)
                                                                                                      caAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACa                                                 <  2:1505542/150‑1 (MQ=255)
                                                                                                       aaaTGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAg                                                <  1:620226/150‑1 (MQ=255)
                                                                                                            cgcAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCggg                                           >  2:1143201/1‑150 (MQ=255)
                                                                                                                ggggTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCag                                       >  1:43038/1‑150 (MQ=255)
                                                                                                                ggggTAAAACGATTCCTGAAGCACCGTGAGTTATACTTTGAATAACTTAAGGAGGTGCAGATGCGTATGACCATAAAAAGATGCGGGAACAGTGCAGGTATGGTCATTCCGAATATCGTAATGAAAGAACTTAACTTCCAGCCGGGGCag                                       >  2:464441/1‑150 (MQ=255)
                                                                                                                  ggTAAACCGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTAGTAGCATAAAAAGTTGGTGTAACAGTGCTGGGATGGTCATTGCCAATATCTTTCTGAAAGAAGTTAACTTACAGCCgggggtgag                                     <  2:1029548/150‑6 (MQ=255)
                                                                                                                    tAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCg                                   >  1:552384/1‑150 (MQ=255)
                                                                                                                     aaaaCGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGt                                  >  2:1188381/1‑150 (MQ=255)
                                                                                                                       aaCGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTgg                                >  1:1150508/1‑150 (MQ=255)
                                                                                                                       aaCGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTgg                                <  2:243590/150‑1 (MQ=255)
                                                                                                                        aCGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGa                               <  2:103012/150‑1 (MQ=255)
                                                                                                                         cGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAg                              <  2:803785/150‑1 (MQ=255)
                                                                                                                          gTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAgg                             <  2:1340009/150‑1 (MQ=255)
                                                                                                                             tCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGaact                                                          <  2:1103904/118‑1 (MQ=255)
                                                                                                                             tCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGaact                                                          >  1:1103904/1‑118 (MQ=255)
                                                                                                                              ccTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCa                         >  1:97396/1‑150 (MQ=255)
                                                                                                                               cTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCagag                                     >  2:9033/1‑137 (MQ=255)
                                                                                                                                      tccttaaGTTATACAAAGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGca                 >  2:1391377/7‑150 (MQ=255)
                                                                                                                                      aCCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGca                 <  2:571968/150‑1 (MQ=255)
                                                                                                                                        cttaaGTTATACAAAGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTgg                                <  2:1236331/129‑1 (MQ=255)
                                                                                                                                        cttaaGTTATACAAAGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTgg                                >  1:1236331/5‑133 (MQ=255)
                                                                                                                                         ttaaGTTATACAAAGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGcaaca              >  1:421474/4‑150 (MQ=255)
                                                                                                                                         gTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGcaaca              >  1:807230/1‑150 (MQ=255)
                                                                                                                                           gAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAAt            <  2:9690/150‑1 (MQ=255)
                                                                                                                                            aGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATc           <  1:1032178/150‑1 (MQ=255)
                                                                                                                                            aGTTATACAAAGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATc           >  1:1447615/1‑150 (MQ=255)
                                                                                                                                             gTTATACAAAGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATCa          >  1:405414/1‑150 (MQ=255)
                                                                                                                                               tataCTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATCAAc        <  2:1425186/150‑1 (MQ=255)
                                                                                                                                                ataCTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATCAACt       >  2:853818/1‑150 (MQ=255)
                                                                                                                                                     ttGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATCAACTGATTc  >  1:880408/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGAGACGACTCGCATTCTGTTTTGTAATTCCCTCACCTTTTGCTTTTCTCTCCGAGCCGCTTTCCATATCTATTAACGCATAAAAAACTCTGCTGGCATTCACAAATGCGCAGGGGTAAAACGTTTCCTGTAGCACCGTGAGTTATACTTTGTATAACTTAAGGAGGTGCAGATGCGTATTACCATAAAAAGATGGGGGAACAGTGCAGGTATGGTCATTCCCAATATCGTAATGAAAGAACTTAACTTACAGCCGGGGCAGAGCGTGGAGGCGCAAGTGAGCAACAATCAACTGATTC  >  NC_000913/4448275‑4448573

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: