Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,640,501 C→A 5.3% D29Y (GAT→TAT)  essQ ← Qin prophage; putative S lysis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,640,5010CA 5.3% 332.6 / 6.7 114D29Y (GAT→TAT) essQQin prophage; putative S lysis protein
Reads supporting (aligned to +/- strand):  ref base C (48/60);  new base A (3/3);  total (51/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-01

Only 100 of 116 total aligned reads displayed.

TAGTCTTGAGTCATTGGATTACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCATATCCTCACCCCGATGTACACGGATGGTGCAATATGTTTGAAAAGATCGGAGTCTACGGGGTAG  >  NC_000913/1640352‑1640648
                                                                                                                                                     |                                                                                                                                                   
tAGTCTTGAGTCATTGGATTACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATc                                                                                                                                                     >  1:150242/1‑150 (MQ=255)
 aGTCTTGAGTCATTGGATTACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCg                                                                                                                                                    <  2:1422690/150‑1 (MQ=255)
 aGTCTTGAGTCATTGGATTACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCg                                                                                                                                                    <  1:1454327/150‑1 (MQ=255)
 aGTCTTGAGTCATTGGATTACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCg                                                                                                                                                    >  1:1253783/1‑150 (MQ=255)
          tCATTGGATTACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGc                                                                                                                                           <  2:1001120/150‑1 (MQ=255)
          tCATTGGATTACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGc                                                                                                                                           >  1:108736/1‑150 (MQ=255)
           cATTGGATTACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCCCTGTGACGGAGTTACTTTATCGAGTAACTGCa                                                                                                                                          >  1:1456752/1‑150 (MQ=255)
           cATTGGATTACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCa                                                                                                                                          >  2:1034120/1‑150 (MQ=255)
           cATTGGATTACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACAACGATTGCAGCCCACTGTGACGGAGTTACTTTTTCGAGTAACTGCa                                                                                                                                          >  1:457486/1‑150 (MQ=255)
            aTTGGATTACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCaa                                                                                                                                         >  1:1091824/1‑150 (MQ=255)
            aTTGGATTACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGCAGTTACTCGATAAAGTAACtccgt                                                                                                                                         <  1:1204891/150‑5 (MQ=18)
             ttGGATTACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTCGATAAAGTAACtccgtc                                                                                                                                        <  1:1560581/150‑6 (MQ=18)
             ttGGATTACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGCAGTTACTCGATAAAGTAACtccgtc                                                                                                                                        <  2:347809/150‑6 (MQ=18)
               ggATTACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTTATCTTGATATAAAGTTTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACCCCGATTGCAGCCCACTGTGACGGAGTTACATTTTCGAGGATCTGCAATGc                                                                                                                                      >  2:158264/1‑150 (MQ=255)
               ggATTACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGAATGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGc                                                                                                                                      >  1:988894/1‑150 (MQ=255)
                   tACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCTGGCCAAAAACCAGGCTACCCAGCACTCCGATTGCAGCCCACTGTTACGGAGTTACTTTATCGAGTAACTGCAATGCGCAg                                                                                                                                  <  1:312277/150‑1 (MQ=255)
                   tACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAg                                                                                                                                  <  1:1113492/150‑1 (MQ=255)
                    aCTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGa                                                                                                                                 >  1:154176/1‑150 (MQ=255)
                     ctctcctctcGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAACAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTGACTTTATCGAGGTACTGCAATGCCCAGaa                                                                                                                                >  2:103718/1‑150 (MQ=255)
                     ctctcctctcGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGaa                                                                                                                                >  2:1552802/1‑150 (MQ=255)
                     ctctcctctcGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGaa                                                                                                                                >  2:1224870/1‑150 (MQ=255)
                     cgctcctctcGCAGCCTTTCGCCTGTGGTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCCGCAGGCCAAAAACCAGGCGACCGAGCACCCCGATTGCAGGGCCCTGTGACGGAGTAACTTTATCGAGTAACTGCACTGCCCAGac                                                                                                                                >  2:378978/3‑149 (MQ=18)
                      tctcctctcGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGaaa                                                                                                                               <  2:1546725/150‑1 (MQ=255)
                      tctcctctcGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGaaa                                                                                                                               >  1:1238366/1‑150 (MQ=255)
                       ctcctctcGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAAc                                                                                                                              >  1:177381/1‑150 (MQ=255)
                         cctctcGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCa                                                                                                                            >  2:1147722/1‑150 (MQ=255)
                         cctctcGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCa                                                                                                                            >  2:368272/1‑150 (MQ=255)
                          ctctcGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAg                                                                                                                           <  2:621793/150‑1 (MQ=255)
                           tctcGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGc                                                                                                                          <  2:721625/150‑1 (MQ=255)
                            ctcGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCa                                                                                                                         <  1:951828/150‑1 (MQ=255)
                                  gCCTTACGCCTGTCTTCTTTAATGTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGGAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTAccc                                                                                                                   <  2:903679/150‑1 (MQ=255)
                                     ttACGCCTGTCTTCTTTAATCTTGTAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACccgcc                                                                                                                <  2:1228329/150‑1 (MQ=255)
                                         gCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCCGAAACCAGCATTACCCGCCGAtg                                                                                                            >  2:1653722/1‑150 (MQ=255)
                                           cTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGAtgtg                                                                                                          <  1:368272/150‑1 (MQ=255)
                                           cTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGAtgtg                                                                                                          <  1:1147722/150‑1 (MQ=255)
                                           cTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGAtgtg                                                                                                          <  2:64577/150‑1 (MQ=255)
                                             gTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGcc                                                                                                        <  1:875777/150‑1 (MQ=255)
                                             gTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGcc                                                                                                        >  2:698007/1‑150 (MQ=255)
                                              tcttctTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCa                                                                                                       <  2:657524/150‑1 (MQ=255)
                                               cttcttTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCCCTGTGACGGAGTTACTTTATCGAGTAACTGCTATGCCCAGAAACCTGCATTACCCGCCGATGTGCCAt                                                                                                      >  2:250976/1‑150 (MQ=255)
                                                 tcttTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAg                                                                                                    >  1:1702192/1‑150 (MQ=255)
                                                    ttAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCa                                                                                                 <  2:182027/150‑1 (MQ=255)
                                                      aaTCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAac                                                                                               <  2:154176/150‑1 (MQ=255)
                                                         catGAAATACAGAGTTGTCAGATACGTCAGCAGGCCAAACCCCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGGTACTTAATCGAGTAACTGCAAAGCCCAGAAACCAGCATTACCCGCCGATGTGCCAAAGGCAACAcc                                                                                            <  1:607648/148‑1 (MQ=255)
                                                         cTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACAcc                                                                                            <  2:150242/150‑1 (MQ=255)
                                                         cTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACAcc                                                                                            >  2:1716060/1‑150 (MQ=255)
                                                            gacaTAAAGATTTGTCAGATACGTCAGCAGGCCAACAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTAATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCtgt                                                                                         <  1:1629089/147‑1 (MQ=255)
                                                             aaaTAAAGATTTGTCAGATACGTCAGCAGGCCAAAACCCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTAATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCtgtt                                                                                        <  1:91634/150‑1 (MQ=255)
                                                               aTAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCtgttgt                                                                                      >  1:1555500/1‑150 (MQ=255)
                                                                tAAAGATTAGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTgttgtt                                                                                     <  2:1411142/150‑1 (MQ=255)
                                                                 aaaGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTa                                                                                    >  2:573934/1‑150 (MQ=255)
                                                                 aaaGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTa                                                                                    <  1:1376387/150‑1 (MQ=255)
                                                                   aGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCCGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAAc                                                                                  <  1:847672/150‑1 (MQ=255)
                                                                    taTTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACt                                                                                 <  2:787016/149‑1 (MQ=255)
                                                                    gATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACt                                                                                 >  2:229214/1‑150 (MQ=255)
                                                                    gATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACt                                                                                 <  1:496116/150‑1 (MQ=255)
                                                                    gATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACt                                                                                 >  2:790933/1‑150 (MQ=255)
                                                                    gATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACt                                                                                 <  2:1733570/150‑1 (MQ=255)
                                                                        tGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATc                                                                             >  1:673645/1‑150 (MQ=255)
                                                                         gTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATcc                                                                            <  2:1671089/150‑1 (MQ=255)
                                                                           cAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCACCACCTGTTGTTAACTTATACAt                                                                          >  2:1048799/1‑150 (MQ=255)
                                                                           cAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCAt                                                                          <  2:607384/150‑1 (MQ=255)
                                                                            aGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCAtt                                                                         <  1:935394/150‑1 (MQ=255)
                                                                             gATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTg                                                                        <  2:1253783/150‑1 (MQ=255)
                                                                             acgccGTCAGCAGGCCAAAAACCAGGATACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTg                                                                        <  1:825496/146‑1 (MQ=255)
                                                                                 cGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGAttt                                                                    <  1:847640/150‑1 (MQ=255)
                                                                                    cagcagGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCTATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCat                                                                 >  2:871869/1‑150 (MQ=255)
                                                                                    cagcagGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCat                                                                 <  1:492712/150‑1 (MQ=255)
                                                                                     agcagGCCAAAAACCAGTCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTATCTGCAATTCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCata                                                                <  1:516963/150‑1 (MQ=255)
                                                                                     agcagGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCata                                                                <  2:1236335/150‑1 (MQ=255)
                                                                                     agcagGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCata                                                                <  1:1552802/150‑1 (MQ=255)
                                                                                     agcagGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCata                                                                <  1:492270/150‑1 (MQ=255)
                                                                                     agcagGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCata                                                                <  1:1224870/150‑1 (MQ=255)
                                                                                      gcagGCGAAAAACAAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCatat                                                               <  2:1247701/150‑1 (MQ=255)
                                                                                      gcagGCCAAAAACCAGGCTACCCAGCACCCCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGACACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCatat                                                               >  1:710113/1‑150 (MQ=255)
                                                                                       cagGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATTGGCAACACCTGTTGTTAACTTATCCATTGATTTCATATc                                                              >  2:670180/1‑150 (MQ=255)
                                                                                           ccAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCATATCCTCa                                                          >  1:542820/1‑150 (MQ=255)
                                                                                           ccAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTGACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTTTCCATTGATTTCATATCCTCa                                                          >  2:1059114/1‑150 (MQ=255)
                                                                                            cAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCATATCCTCAc                                                         >  1:1602716/1‑150 (MQ=255)
                                                                                                  ccAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCACCACCTGTTGTTAACTTATCCATTGATTTCATATCCTCACCCCGAt                                                   >  2:127215/1‑150 (MQ=255)
                                                                                                  ccAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTAATTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTGCATATCCTCACCCCGAt                                                   >  1:73277/1‑150 (MQ=255)
                                                                                                    aGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCATATCCTCACCCCGATGt                                                 <  2:108736/150‑1 (MQ=255)
                                                                                                         aCCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCATATCCTCACCCCGATGTACACg                                            <  1:698007/150‑1 (MQ=255)
                                                                                                         aCCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCATATCCTCACCCCGATGTACACg                                            <  1:959791/150‑1 (MQ=255)
                                                                                                          cccAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTACCTTATCCATTGATTTCATATCCTCACCCCGATGTACACgg                                           >  2:1017235/1‑150 (MQ=255)
                                                                                                               cacaCCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCATATCCTCACCCCGATGTACACGGATGGt                                      <  2:1091824/150‑1 (MQ=255)
                                                                                                                 cacCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCTGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCATATCCTCACCCCGATGTACACGGATGGTGc                                    <  2:287543/150‑1 (MQ=255)
                                                                                                                          cAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCATATCCTCACCCCGATGTACACGGATGGTGCAATATGttt                           >  2:1018287/1‑150 (MQ=255)
                                                                                                                           aGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCATATCCTCACCCCGATGTACACGGATGGTGCAATATGTTTg                          <  2:1650776/150‑1 (MQ=255)
                                                                                                                           aGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCATATCCTCACCCCGATGTACACGGATGGTGCAATATGTTTg                          <  1:1048728/150‑1 (MQ=255)
                                                                                                                            tcccACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCATATCCTCACCCCGATGTACACGGATGGTGCAATATGTTTGa                         >  1:1069577/2‑150 (MQ=255)
                                                                                                                               cACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCATATCCTCACCCCGATGTACACGGATGGTGCAATATGTTTGaaaa                      <  1:669572/150‑1 (MQ=255)
                                                                                                                                   gtgACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCATATCCTCACCCCGATGTACACGGATGGTGCAATATGTTTGAAAAGATc                  <  1:1034120/150‑1 (MQ=255)
                                                                                                                                    tgACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCATATCCTCACCCCGATGTACACGGATGGTGCAATATGTTTGAAAAGATCg                 <  1:546121/150‑1 (MQ=255)
                                                                                                                                     gACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATCACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCATATCCTCACCCCGATGTACACGGATGGTGCAATATGTTTGAAAAGATCgg                <  2:177381/150‑1 (MQ=255)
                                                                                                                                       tgcaGTTACTCGATAAAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCATATCCTCACCCCGATGTACACGGATGGTGCAATATGTTTGAAAAGATCGGAg              >  2:1561486/4‑150 (MQ=255)
                                                                                                                                        gcaGTTACTCGATAAAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCATATCCTCACCCCGATGTACACGGATGGTGCAATATGTTTGAAAAGATCGGAGt             >  1:1218600/3‑150 (MQ=255)
                                                                                                                                        gcaGTTACTCGATAAAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCATATCCTCACCCCGATGTACACGGATGGTGCAATATGTTTGAAAAGATCGGAGt             >  1:730817/3‑150 (MQ=255)
                                                                                                                                                 ttATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCATATCCTCACCCCGATGTACACGGATGGTGCAATATGTTTGAAAAGATCGGAGTCTACGGGGt    >  1:1704277/1‑150 (MQ=255)
                                                                                                                                                   aTCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCATATCCTCACCCCCATGTACACGGATGGTGCAATATGTTTGAAAAGATCGGAGTCTACGGGGTAg  <  2:1324457/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TAGTCTTGAGTCATTGGATTACTCTCCTCTCGCAGCCTTACGCCTGTCTTCTTTAATCTTGAAATAAAGATTTGTCAGATACGTCAGCAGGCCAAAAACCAGGCTACCCAGCACACCGATTGCAGCCCACTGTGACGGAGTTACTTTATCGAGTAACTGCAATGCCCAGAAACCAGCATTACCCGCCGATGTGCCATAGGCAACACCTGTTGTTAACTTATCCATTGATTTCATATCCTCACCCCGATGTACACGGATGGTGCAATATGTTTGAAAAGATCGGAGTCTACGGGGTAG  >  NC_000913/1640352‑1640648

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