Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,326,887 A→T 7.1% intergenic (‑153/‑248) yjdM ← / → crfC conserved protein YjdM/clamp‑binding sister replication fork colocalization protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,326,8870AT 7.1% 487.6 / 16.9 154intergenic (‑153/‑248)yjdM/crfCconserved protein YjdM/clamp‑binding sister replication fork colocalization protein
Reads supporting (aligned to +/- strand):  ref base A (78/65);  new base T (7/4);  total (85/69)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 159 total aligned reads displayed.

ATTTCCTCAAAAGTCAGCAGGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAACACACTAAATAT  >  NC_000913/4326739‑4327036
                                                                                                                                                    |                                                                                                                                                     
 tttCCTAAAAAGTCAGCAGGGCAAAAGCGACAAAAAACGGCAGATAGCCGAAAAAGGCCCCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGaaaa                                                                                                                                                     >  2:1066177/1‑150 (MQ=255)
  ttCCTCAAAAGTCAGCAGGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAaat                                                                                                                                                    <  1:727231/150‑1 (MQ=255)
     cTCAAAAGTCAGCAGGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAaataat                                                                                                                                                 >  2:1368330/1‑150 (MQ=255)
         aaaGTCAGCAGGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTTCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCttt                                                                                                                                             >  1:39668/1‑150 (MQ=255)
         aaaGTCAGCAGGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCttt                                                                                                                                             <  2:1409220/150‑1 (MQ=255)
            gTCAGCAGGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTg                                                                                                                                          <  2:840474/150‑1 (MQ=255)
              cagcagGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGAGGACGATCGTTTTCATTTGAAAGTAAACGCAAAAAGATTATTTTCTTTTTGcg                                                                                                                                        <  2:988075/150‑2 (MQ=255)
              cagcagGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGATTATTTTCTTTTTGcg                                                                                                                                        <  1:1529098/150‑2 (MQ=255)
              aagcagGGCAAAATCGCCAAAAAACGGAAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGTGGTACTTCTGATTTTAGACGACGATCTTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGcc                                                                                                                                        <  2:1103805/149‑1 (MQ=255)
                gcagGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGATTATTTTCTTTTTGcgtt                                                                                                                                      <  1:493355/150‑4 (MQ=255)
                gcagGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCctc                                                                                                                                      >  1:1432701/1‑150 (MQ=255)
                 cagGGCAAAAGCGCCAAAAATCGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGTGGTACTGCTTATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGATTATTTTCTTTTTGcgttt                                                                                                                                     <  2:1192392/150‑5 (MQ=255)
                     gCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCtattt                                                                                                                                 >  1:614918/1‑150 (MQ=255)
                         aaGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATttatta                                                                                                                             <  2:1126237/150‑1 (MQ=255)
                             gcCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTa                                                                                                                         >  1:430744/1‑150 (MQ=255)
                              ccAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTaa                                                                                                                        <  2:1317489/150‑1 (MQ=255)
                              ccAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTAATTATTATTTaa                                                                                                                        <  2:57112/150‑1 (MQ=255)
                               cAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAg                                                                                                                       >  1:426219/1‑150 (MQ=255)
                                     aCGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATg                                                                                                                 >  1:130110/1‑150 (MQ=255)
                                      cGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGt                                                                                                                >  1:55133/1‑150 (MQ=255)
                                               gCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGctct                                                                                                       >  1:123557/1‑150 (MQ=255)
                                               gCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGctct                                                                                                       >  2:878915/1‑150 (MQ=255)
                                               gCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGctct                                                                                                       >  1:890747/1‑150 (MQ=255)
                                               gCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGatct                                                                                                       >  2:1200141/1‑150 (MQ=255)
                                                  gAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACt                                                                                                    <  2:123557/150‑1 (MQ=255)
                                                        ggCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCt                                                                                              >  2:1070769/1‑150 (MQ=255)
                                                         gCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCtt                                                                                             >  1:592476/1‑150 (MQ=255)
                                                         gCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCtt                                                                                             >  1:938707/1‑150 (MQ=255)
                                                          ccGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCttt                                                                                            >  2:1390692/1‑150 (MQ=255)
                                                          ccGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCttt                                                                                            >  2:491274/1‑150 (MQ=255)
                                                             cAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCtttttt                                                                                         <  1:11381/150‑1 (MQ=255)
                                                              aaTTATACACAAAATCCCTGAGTTTGTCGGGTGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCttttttt                                                                                        <  2:736125/150‑1 (MQ=255)
                                                                ttATACACAAAATCCCTGAGTTTGTCGGGGGTACTTCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTaa                                                                                      <  1:790915/150‑1 (MQ=255)
                                                                 tataCACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAAt                                                                                     <  2:283165/150‑1 (MQ=255)
                                                                   aacacAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTc                                                                                   <  2:919244/149‑1 (MQ=255)
                                                                     cacaAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTcac                                                                                 >  1:444137/1‑150 (MQ=255)
                                                                     cacaAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTcaa                                                                                 >  1:946494/1‑149 (MQ=255)
                                                                        aaaaTCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATg                                                                              <  2:125122/150‑1 (MQ=255)
                                                                           aTCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGata                                                                           >  2:1436619/1‑150 (MQ=255)
                                                                              ccTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATAttt                                                                        >  2:1483639/1‑150 (MQ=255)
                                                                              ccTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATAttt                                                                        >  2:9314/1‑150 (MQ=255)
                                                                               cTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATAtttt                                                                       >  1:944935/1‑150 (MQ=255)
                                                                                 gAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTa                                                                     >  1:1045202/1‑150 (MQ=255)
                                                                                   gTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTca                                                                                  <  1:1220663/135‑1 (MQ=255)
                                                                                   gTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTca                                                                                  >  2:1220663/1‑135 (MQ=255)
                                                                                   gTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATc                                                                   >  2:1150670/1‑150 (MQ=255)
                                                                                       gTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTc                                                               <  1:1368330/150‑1 (MQ=255)
                                                                                        tCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCt                                                              >  1:195282/1‑150 (MQ=255)
                                                                                           ggggTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGACTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGaa                                                           >  1:1357490/1‑149 (MQ=255)
                                                                                              gTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGtt                                                        <  2:890747/150‑1 (MQ=255)
                                                                                                aCTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGtttt                                                      <  1:1072886/150‑1 (MQ=255)
                                                                                                aCTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGtttt                                                      <  1:389924/150‑1 (MQ=255)
                                                                                                 ctgctgATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTa                                                     <  2:944935/150‑1 (MQ=255)
                                                                                                  tgctgATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTaa                                                    >  1:469999/1‑150 (MQ=255)
                                                                                                   gctgATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTaaa                                                   >  2:340106/1‑150 (MQ=255)
                                                                                                   gctgATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTaaa                                                   >  1:177387/1‑150 (MQ=255)
                                                                                                   gctgATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGACTCTATTTATTATTTAAGATAATTTTAATTCTCTAATTGTTCTTTTTTTAATTAACATGATATTTTTATCTTTCTGACGTTTTaaa                                                   >  2:501644/1‑150 (MQ=255)
                                                                                                        ttttAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTAt                                              >  2:1429855/1‑150 (MQ=255)
                                                                                                        ttttAGACGACGATCGCTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTAt                                              >  2:1011837/1‑150 (MQ=255)
                                                                                                         tttAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATc                                             <  2:402347/150‑1 (MQ=255)
                                                                                                             tacgacgATCGTTCTCATTTGAAATGAAACGCAAAAAGAAAATAATCTTTTTTCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCagaa                                         <  1:1514174/149‑1 (MQ=255)
                                                                                                             gacgacGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAAt                                                                                     <  1:342010/106‑1 (MQ=255)
                                                                                                             gacgacGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAAt                                                                                     >  2:342010/1‑106 (MQ=255)
                                                                                                             gacgacGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCagaa                                         >  1:235162/1‑150 (MQ=255)
                                                                                                                 acgaTCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACg                                     <  2:1503539/150‑1 (MQ=255)
                                                                                                                  cgaTCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGt                                    >  1:1091397/1‑150 (MQ=255)
                                                                                                                  cgaTCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGt                                    <  2:430744/150‑1 (MQ=255)
                                                                                                                  cgaTCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGt                                    <  2:1535793/150‑1 (MQ=255)
                                                                                                                   gaTCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGtt                                   <  2:873067/150‑1 (MQ=255)
                                                                                                                    aTCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGtta                                  <  2:177387/150‑1 (MQ=255)
                                                                                                                    aTCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGtta                                  >  1:394061/1‑150 (MQ=255)
                                                                                                                      cGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGttatt                                <  2:953328/150‑1 (MQ=255)
                                                                                                                      cGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGttatt                                <  2:592476/150‑1 (MQ=255)
                                                                                                                        ttttCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCCGACGTTTTAAACGTATCAGAAGACGttattta                              >  2:1001193/1‑150 (MQ=255)
                                                                                                                           tCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTtatttatta                           >  1:1464477/1‑149 (MQ=255)
                                                                                                                           tCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTc                           <  2:1224417/150‑1 (MQ=255)
                                                                                                                           tCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTc                           <  2:809615/150‑1 (MQ=255)
                                                                                                                           tCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTc                           <  1:998245/150‑1 (MQ=255)
                                                                                                                             aTTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTg                         >  2:439800/1‑150 (MQ=255)
                                                                                                                             aTTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTg                         <  2:1058398/150‑1 (MQ=255)
                                                                                                                              tttGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGa                        >  1:1260877/1‑150 (MQ=255)
                                                                                                                              tttGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGa                        >  1:690194/1‑150 (MQ=255)
                                                                                                                                tGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATa                      <  1:9314/150‑1 (MQ=255)
                                                                                                                                 gAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATaa                     >  2:644439/1‑150 (MQ=255)
                                                                                                                                 gAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATaa                     >  1:859389/1‑150 (MQ=255)
                                                                                                                                    aGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATaaaaa                  >  1:1509752/1‑150 (MQ=255)
                                                                                                                                    aGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATaaaaa                  >  2:1019152/1‑150 (MQ=255)
                                                                                                                                       aaaCGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAaaca               <  2:39668/150‑1 (MQ=255)
                                                                                                                                        gaggCAAAAAGATTATTTTCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTGTTTTTTTAATTCACATGATTTTTTTATCTTTATGACGTTTTAAACGTAGCAGAAGACGTTATTTATTGTGATAAAAaacaa              >  2:387722/4‑150 (MQ=255)
                                                                                                                                          ggCAAAAAGATTATTTTCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAacaaca            >  1:256754/2‑150 (MQ=255)
                                                                                                                                          cGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAacaaca            >  1:97330/1‑150 (MQ=255)
                                                                                                                                           gCAAAAAGATTATTTTCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAacac           >  2:586314/1‑150 (MQ=255)
                                                                                                                                           gCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAacac           >  1:549815/1‑150 (MQ=255)
                                                                                                                                             aaaaaGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAacacac         >  2:22815/1‑150 (MQ=255)
                                                                                                                                              aaaaGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAACACACt        >  1:575367/1‑150 (MQ=255)
                                                                                                                                               aaaGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTAATTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGAAGTTTTAAACGTATCAGAAGAAGTTATTTATTCTGATAAAAAAAAACACACTa       >  2:305332/1‑150 (MQ=255)
                                                                                                                                                  gAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAACACACTAAat    <  1:340106/150‑1 (MQ=255)
                                                                                                                                                  gAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAACACACTAAat    <  1:1070769/150‑1 (MQ=255)
                                                                                                                                                   aaaaTAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAACACACTAAata   >  1:1021438/1‑150 (MQ=255)
                                                                                                                                                    aaaTAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAACACACTAAatat  <  2:614918/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATTTCCTCAAAAGTCAGCAGGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAACACACTAAATAT  >  NC_000913/4326739‑4327036

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: