Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,326,887 A→T 11.7% intergenic (‑153/‑248) 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,8870AT11.7% 218.4 / 17.2 92intergenic (‑153/‑248)yjdM/crfCzinc‑ribbon family protein/clamp‑binding sister replication fork colocalization protein, dynamin‑related
Reads supporting (aligned to +/- strand):  ref base A (39/42);  new base T (6/5);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 7.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.20e-01

GATTTCCTCAAAAGTCAGCAGGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAACACACTAAAT  >  NC_000913/4326738‑4327034
                                                                                                                                                     |                                                                                                                                                   
gATTTCCTCAAAAGTCAGCAGGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGaa                                                                                                                                                     >  2:596976/1‑150 (MQ=255)
 aTTACCTCAACACGCTGCAGGGCTAAAGAGACTCAAAACGTCAGATTTCCGAAAAAGGCCGCAATTATAAACAAAATCCCTGACTTTGTGGGGGGTCCTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAATGaat                                                                                                                                                    <  2:238571/150‑2 (MQ=255)
  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                                                                      <  2:251642/150‑1 (MQ=255)
                                                                               ccTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATAttt                                                                      >  1:957834/1‑150 (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)
                                                                                                                                                     |                                                                                                                                                   
GATTTCCTCAAAAGTCAGCAGGGCAAAAGCGCCAAAAAACGGCAGATTGCCGAAAAAGGCCGCAATTATACACAAAATCCCTGAGTTTGTCGGGGGTACTGCTGATTTTAGACGACGATCGTTTTCATTTGAAAGGAAACGCAAAAAGAAAATAATCTTTTTGCCTCTATTTATTATTTAAGATAATGTTAATGCTCTACTTGTTCTTTTTTTAATTCACATGATATTTTTATCTTTCTGACGTTTTAAACGTATCAGAAGACGTTATTTATTCTGATAAAAAACAACACACTAAAT  >  NC_000913/4326738‑4327034

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