Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,326,892 A→T 11.2% intergenic (‑158/‑243) yjdM ← / → crfC zinc‑ribbon family protein/clamp‑binding sister replication fork colocalization protein, dynamin‑related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,326,8920AT11.2% 221.7 / 15.3 89intergenic (‑158/‑243)yjdM/crfCzinc‑ribbon family protein/clamp‑binding sister replication fork colocalization protein, dynamin‑related
Reads supporting (aligned to +/- strand):  ref base A (36/43);  new base T (6/4);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

CCTCAAAAGTCAGCAGGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAACACACTAAATATAAAA  >  NC_000913/4326743‑4327040
                                                                                                                                                     |                                                                                                                                                    
ccTCAAAAGTCAGCAGGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATGTTAGCCGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAaataa                                                                                                                                                      >  2:781377/1‑150 (MQ=255)
 cTCAAAAGTCAGCAGGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAaataat                                                                                                                                                     >  1:1295712/1‑150 (MQ=255)
 cTCAAAAGTCAGCAGGGCAAAAGCGCCAAAAAACAGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAaataat                                                                                                                                                     >  1:752857/1‑150 (MQ=255)
  tCAAAAGTCAGCAGGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATc                                                                                                                                                    <  2:437148/150‑1 (MQ=255)
     aaaGTCAGCAGGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGGTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCttt                                                                                                                                                 >  1:622154/1‑150 (MQ=255)
      aaGGCAGCAGGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGACGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTAATGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCtttt                                                                                                                                                <  1:476090/150‑1 (MQ=255)
        gTCAGCAGGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGCCGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTg                                                                                                                                              >  2:125047/1‑150 (MQ=255)
        gTCAGCAGGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTg                                                                                                                                              <  1:467991/150‑1 (MQ=255)
         tCAGCAGGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGATTATTTTCTTTTTGc                                                                                                                                             <  1:837943/150‑1 (MQ=255)
          cagcagGGCAAAAGCGCCAAAAACCGGCAGATTGCCTAATAAGGCCGCAAATATAGACAATATGCCTGAGTTTGTCGGTTGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAATCGCAAAAAGATTATTTTCTTTTTGcg                                                                                                                                            <  1:52413/150‑2 (MQ=255)
          cagcagGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGATTATTTTCTTTTTGcg                                                                                                                                            <  1:1256838/150‑2 (MQ=255)
           agcagGTCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCct                                                                                                                                           <  1:267833/150‑1 (MQ=255)
           agcagGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTCCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGcat                                                                                                                                           >  2:1025776/1‑148 (MQ=255)
            gcagGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGATTATTTTCTTTTTGcgtt                                                                                                                                          <  2:651541/150‑4 (MQ=255)
               gggCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCtat                                                                                                                                       >  2:348123/1‑150 (MQ=255)
                     aaGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATttatta                                                                                                                                 <  1:1241435/150‑1 (MQ=255)
                      aGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTtattat                                                                                                                                <  2:1066271/150‑1 (MQ=255)
                       gcgcCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTattatt                                                                                                                               <  1:366438/150‑1 (MQ=255)
                       gcgcCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTattatt                                                                                                                               <  2:1117974/150‑1 (MQ=255)
                        cgcCAAAAAACGGCAGATTGACGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTAttt                                                                                                                              >  1:308939/1‑150 (MQ=255)
                         gcCAAATAACGGAAGATTGCCGAAAAAGGCCGCAATTATACACAAATTCCCTGAGTTTTTCGGGGGTACTTCTTATTTTAGACGACGATCGTTTTCTTTTGAAAGGAATCGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTa                                                                                                                             <  1:665918/150‑1 (MQ=255)
                          ccAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTaa                                                                                                                            >  1:934848/1‑150 (MQ=255)
                          ccAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCATCTATTTATTATTTaa                                                                                                                            >  1:963910/1‑150 (MQ=255)
                            aaaaaaCGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGa                                                                                                                          >  1:637186/1‑150 (MQ=255)
                                aaCGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAAt                                                                                                                      <  2:377910/150‑1 (MQ=255)
                                  cGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATATTGt                                                                                                                    >  2:182608/1‑150 (MQ=255)
                                        aTTGCCGAAACAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAACGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGc                                                                                                              <  2:508443/150‑1 (MQ=255)
                                            ccGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTa                                                                                                          >  1:138798/1‑150 (MQ=255)
                                            ccGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTa                                                                                                          >  1:1206610/1‑150 (MQ=255)
                                            ccGAAAAAGGCCGCAAATATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTATTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTAATTATTATTTAAGATAATGTTAATGCTCTa                                                                                                          >  1:753949/1‑150 (MQ=255)
                                             cGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGACTCTATTTATTATTTAAGATAATGTTAATGCTCTAc                                                                                                         >  1:228210/1‑150 (MQ=255)
                                                  aaGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGtt                                                                                                    >  1:491150/1‑150 (MQ=255)
                                                   aGGCCTCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTc                                                                                                   <  2:1070233/150‑1 (MQ=255)
                                                    ggCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGTGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTcn                                                                                                  <  2:161968/150‑2 (MQ=255)
                                                    ggCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCt                                                                                                  >  2:806665/1‑150 (MQ=255)
                                                    ggCCGCAATTAAACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCt                                                                                                  <  2:639177/150‑1 (MQ=255)
                                                     gCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCtt                                                                                                 >  2:975689/1‑150 (MQ=255)
                                                     gCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTATATTTATTATTTAAGATAATGTTAATGATCTACTTGTTCtt                                                                                                 >  2:1007110/1‑150 (MQ=255)
                                                        gCAATTATACACAGAATACCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCttttt                                                                                              <  1:1113166/150‑1 (MQ=255)
                                                          aaTTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCttttttt                                                                                            <  2:1229259/150‑1 (MQ=255)
                                                           aTTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTa                                                                                           <  2:308939/150‑1 (MQ=255)
                                                             tataCACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAAt                                                                                         >  1:1028311/1‑150 (MQ=255)
                                                                   caAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATAGTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTATATTTAATTCACAt                                                                                   >  2:464245/1‑150 (MQ=255)
                                                                    aaaaTCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAATAAAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATg                                                                                  <  2:463909/150‑1 (MQ=255)
                                                                          ccTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATAttt                                                                            >  1:957834/1‑150 (MQ=255)
                                                                          ccTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATAttt                                                                            <  2:251642/150‑1 (MQ=255)
                                                                              aGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTAt                                                                        <  2:138798/150‑1 (MQ=255)
                                                                               gTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTAta                                                                       >  2:1289603/1‑149 (MQ=255)
                                                                                     cGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTg                                                                 <  2:640069/150‑1 (MQ=255)
                                                                                       ggggTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGAc                                                               <  2:484875/150‑1 (MQ=255)
                                                                                           tACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGttt                                                           <  2:653256/150‑1 (MQ=255)
                                                                                             ctgctgATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTa                                                         >  1:372702/1‑150 (MQ=255)
                                                                                               gctgATTTTAGACGACGATCGTTTTCAGTTGAAAGGAAACGCAAAAAGAAAATAATCTTATTGACTCTATTTATTATATAAGATAATGTAAATGCTCTACTTGTTCTTTTTTTAATTCAAATGATATTTTTATCTTTCTGACGTTTTaaa                                                       >  1:86660/1‑150 (MQ=255)
                                                                                                 tgATTTTAGACGACTATATTTTTAATTTGAAAGGAAACGCATAAAGATAATAATATTTTTGCCTCTATTTATTATTTAAGATTATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACg                                                     <  1:1025776/150‑1 (MQ=255)
                                                                                                       tAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCag                                               >  1:343457/1‑150 (MQ=255)
                                                                                                       tAGACGACGATCGTTTACATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCat                                               >  1:306770/1‑149 (MQ=255)
                                                                                                         gacgacGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCag                                               >  1:853813/1‑148 (MQ=255)
                                                                                                         gacgacGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCag                                               <  2:853813/148‑1 (MQ=255)
                                                                                                           cgacgaTCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTTCTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTATGACGTTTTAAACGTATCagaaga                                           >  2:1252608/1‑150 (MQ=255)
                                                                                                             acgaTCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACg                                         <  2:901209/150‑1 (MQ=255)
                                                                                                              cgaTCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGt                                        <  2:637186/150‑1 (MQ=255)
                                                                                                               gaTCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATATTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGtt                                       <  1:746810/150‑1 (MQ=255)
                                                                                                                       tCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTc                               <  1:781377/150‑1 (MQ=255)
                                                                                                                        cATTTTAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCt                              <  1:426686/150‑1 (MQ=255)
                                                                                                                        cATTTGAAAGGAAACTCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAATATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCt                              <  1:242946/150‑1 (MQ=255)
                                                                                                                        cATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCt                              <  1:1103659/150‑1 (MQ=255)
                                                                                                                            tGAAAGTAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATa                          <  1:512686/150‑1 (MQ=255)
                                                                                                                            tGAAAGGAAACGCAAAAATAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATa                          <  2:459606/150‑1 (MQ=255)
                                                                                                                             taaaGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATaa                         <  2:361220/149‑1 (MQ=255)
                                                                                                                              aaaGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATaaa                        >  1:684160/1‑150 (MQ=255)
                                                                                                                               aaGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATaaaa                       >  2:379122/1‑150 (MQ=255)
                                                                                                                                aGGATACGCAAAAAGAAAATAATCTATTTGCCTGTATTTATTATTTAAGATAATGTTAATGCTCTTCTTGTTCTTTTTTTAATTCACATGATATTTTTTTCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTTATaaaat                      <  1:1104679/150‑2 (MQ=255)
                                                                                                                                aGGAAACTCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTATGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATaaaaa                      <  1:4922/150‑1 (MQ=255)
                                                                                                                                 ggAAACGCTAAAAGAAAATAATCTTTTAGCCTCTATTTATTATTTAAGATAATGTTTATGCTCTACTTGTTCTTTTTTTAATTCACATTATATTTTTATCAATCTGACGTTTTAAACGTATCAGAAGACGCTATTTTTTCTGaaaaaaaa                     <  1:1134846/150‑1 (MQ=255)
                                                                                                                                   agaggCAAAAAGATTATTTTCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTATTTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAaaca                   >  1:1010513/5‑150 (MQ=255)
                                                                                                                                   aaaCGCAAAAAGAAAATAATCTTTTTGGCTCTATTTATTATTTAAGATAATGTTAATGCTCTTCTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAaaca                   <  2:820208/150‑1 (MQ=255)
                                                                                                                                   aaaCGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTATTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAaaca                   >  2:45188/1‑150 (MQ=255)
                                                                                                                                     aggCAAAAAGATTATTTTCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAaacaac                 >  1:390872/3‑150 (MQ=255)
                                                                                                                                      ggCAAAAAGATTATTTTCTTTTTGCCTCTATTTTTTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTAACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAacaaca                >  2:397333/2‑150 (MQ=255)
                                                                                                                                      ggCAAAAAGATTATTTTCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAacaaca                >  2:1088392/2‑150 (MQ=255)
                                                                                                                                      cGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAacaaca                >  2:700635/1‑150 (MQ=255)
                                                                                                                                       tcAAAAAGAAAATAATCTTTTTTCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAacac               <  1:125047/149‑1 (MQ=255)
                                                                                                                                       gCAAAAAGATTATTTTCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAacac               >  2:656479/1‑150 (MQ=255)
                                                                                                                                        cAAAAAGATTATTTTCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAacaca              >  2:1312315/1‑150 (MQ=255)
                                                                                                                                        cAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAATCGTATCAGAAGACGTTATTTTTTCTGATAAAAAACAacaca              >  1:49059/1‑150 (MQ=255)
                                                                                                                                        cAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTAACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAacaca              >  1:443382/1‑150 (MQ=255)
                                                                                                                                         aaaaaGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAacacac             <  1:806665/150‑1 (MQ=255)
                                                                                                                                         aaaaaGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAacacac             <  2:1237205/150‑1 (MQ=255)
                                                                                                                                           aaaGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAACACACTa           <  2:934848/150‑1 (MQ=255)
                                                                                                                                           aaaGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTGCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACTTTTTAAACGTATCAGATGACGTTATATATTCTGATAAAAAACAACACACTa           >  2:160711/1‑150 (MQ=255)
                                                                                                                                            aaGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAACACAATaa          >  2:99788/1‑150 (MQ=255)
                                                                                                                                              gAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAACACACTAAat        >  2:926265/1‑150 (MQ=255)
                                                                                                                                                 aataatCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAACAACAACACCCTAAAtata     >  2:852299/1‑150 (MQ=255)
                                                                                                                                                  ataatCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAACACACTAAATATaa    <  2:737778/150‑1 (MQ=255)
                                                                                                                                                  ataatCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAACACACTAAATATaa    >  1:10394/1‑150 (MQ=255)
                                                                                                                                                    aatCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAACACACTAAATATaaaa  <  2:661587/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCTCAAAAGTCAGCAGGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAACACACTAAATATAAAA  >  NC_000913/4326743‑4327040

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