Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,326,888 A→T 12.1% intergenic (‑154/‑247) 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,8880AT12.1% 218.2 / 17.3 89intergenic (‑154/‑247)yjdM/crfCzinc‑ribbon family protein/clamp‑binding sister replication fork colocalization protein, dynamin‑related
Reads supporting (aligned to +/- strand):  ref base A (36/42);  new base T (6/5);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

TTTCCTCAAAAGTCAGCAGGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAACACACTAAATATA  >  NC_000913/4326740‑4327037
                                                                                                                                                    |                                                                                                                                                     
tttCCTCAAAAGTCAGCAGGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAACTCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGaaaa                                                                                                                                                      <  1:756375/150‑1 (MQ=255)
   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:1206610/1‑150 (MQ=255)
                                               ccGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTa                                                                                                       >  1:138798/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                                            <  2:853813/148‑1 (MQ=255)
                                                                                                            gacgacGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCag                                            >  1:853813/1‑148 (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)
                                                                                                                                                    |                                                                                                                                                     
TTTCCTCAAAAGTCAGCAGGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAACACACTAAATATA  >  NC_000913/4326740‑4327037

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: ‑