Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,531,418 A→G 5.1% intergenic (‑185/‑45) zipA ← / → cysZ cell division protein ZipA/sulfate:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,531,4180AG 5.1% 269.8 / 2.6 80intergenic (‑185/‑45)zipA/cysZcell division protein ZipA/sulfate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (39/35);  new base G (2/2);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GGCTGGTCCAGAAACCATGTACCAGTAAAGCGATTATGGCGATCGCGCCAACAATGATTAATATCAGACGCAAATCCTGCATCATTATATTCTCTGTTGTTCTAACACCTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGT  >  NC_000913/2531149‑2531681
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                       
ggCTGGTCCAGAAACCATGTACCAGTAAAGCGATTATGGCGATCGCGCCAACAATGATTAATATCAGACGCAAATCCTGCATCATTATATTCTCTGTTGTTCTAACACCTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAATTGTCCATTGCGCAATTACCCGGTAAAATACg                                                                                                                                                                                                                                                                       >  1:629966/1‑272 (MQ=255)
        cAGAAACCATGTACCAGTAAAGCGATTATGGCGATCGCGCCAACAATGATTAATATCAGACGCAAATCCTGCATCATTATATTCTCTGTTGTTCTAACACCTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAAttt                                                                                                                                                                                                                                                               >  1:376809/1‑272 (MQ=255)
                                      gCGATCGCGCCAACAATGATTAATATCAGACGCAAATCCTGCATCATTATATTCTCTGTTGTTCTAACACCTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTg                                                                                                                                                                                                                                                          >  1:362634/1‑247 (MQ=255)
                                       cGATCGCGCCAACAATGATTAATATCAGACGCAAATCCTGCATCATTATATTCTCTGTTGTTCTAACACCTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTgg                                                                                                                                                                                                                                                         <  2:362633/247‑1 (MQ=255)
                                       cGATCCGGCAAACATTGATTATTATCTGACGAAAATCCTGCATCATTATATTCTCTGTTGTTCTAACACCTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTg                                                                                                                                                                                                                                <  2:224019/272‑1 (MQ=255)
                                                caacaaTGATTAATATCAGACGCAAATCCTGCATCATTATATTCTCTGTTGTTCTAACACCTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCTAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTGTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGtt                                                                                                                                                                                                                       <  1:865931/272‑1 (MQ=255)
                                                     aattttAATATCATACGCAAATCCTGCATCATTATTTTCTCCGTTGTTCTAACACCTTGCCACAACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTtcatc                                                                                                                                                                                                                  <  2:376808/268‑1 (MQ=255)
                                                      tGATTAATATCAGACGCAAATCCTGCATCATTATATTCTCTGTTGTTCTAACACCTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTtcatc                                                                                                                                                                                                                  <  1:166740/271‑1 (MQ=255)
                                                      tGATTAATATCAGACGCAAATCCTGCATCATTATATTCTCTGTTGTTCTAACACCTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTAATTATggactgtatc                                                                                                                                                                                                                 >  1:686337/1‑264 (MQ=255)
                                                                                                           ccTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATGCGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGtt                                                                                                                                                                                        >  2:466209/1‑244 (MQ=255)
                                                                                                            cTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATGCGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGttt                                                                                                                                                                                       <  1:466210/244‑1 (MQ=255)
                                                                                                               gCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTtcat                                                                                                                                                                                                                   >  1:231458/1‑213 (MQ=255)
                                                                                                                ccaccaCGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTtcatc                                                                                                                                                                                                                  <  2:231458/213‑1 (MQ=255)
                                                                                                                       ggCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATg                                                                                                                                                                                                                          >  1:374799/1‑198 (MQ=255)
                                                                                                                        gCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATgg                                                                                                                                                                                                                         <  2:374798/198‑1 (MQ=255)
                                                                                                                         cAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCCAGGCTGGAAGCTCGTCTCGCAACCTg                                                                                                                                              >  1:42590/1‑272 (MQ=255)
                                                                                                                                           gagTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTg                                                                                                                                                                                                                                <  2:686898/172‑1 (MQ=255)
                                                                                                                                           gagTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTg                                                                                                                                                                                                                                >  1:686899/1‑172 (MQ=255)
                                                                                                                                           gagTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCCCGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCCCGTCTCGCCACCTGGGAGTCGGTGTTTCGttt                                                                                                                            >  2:952536/1‑271 (MQ=255)
                                                                                                                                           gagTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTa                                                                                                                            >  1:788449/1‑272 (MQ=255)
                                                                                                                                                   gCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATc                                                                                                                                                                                                         >  1:1039297/1‑187 (MQ=255)
                                                                                                                                                   gCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAg                                                                                                                                                                                             >  2:970919/1‑199 (MQ=255)
                                                                                                                                                   gCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCACGGCTGGAAGCTCGTCTCGCACCCTGGGATTCGTCGTTTCGTTATTTTACCg                                                                                                                    >  1:1092307/1‑272 (MQ=255)
                                                                                                                                                   gCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCGCGCAACCTGGGATTCGGCGTTCCGTTATTTTACCg                                                                                                                    >  1:491411/1‑272 (MQ=255)
                                                                                                                                                   gCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCCCGTCCCGCAACCTGGGATTCGGCGTTTCGTTATTTTTccc                                                                                                                    >  2:377562/1‑271 (MQ=255)
                                                                                                                                                    ccGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCt                                                                                                                                                                                                        <  2:1039296/187‑1 (MQ=255)
                                                                                                                                                    ccGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGc                                                                                                                                                                                            <  1:970920/199‑1 (MQ=255)
                                                                                                                                                    ccGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCGAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTGCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGtt                                                                                                                                                                                                                       <  1:598809/172‑1 (MQ=255)
                                                                                                                                                         tACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCt                                                                                                                                                                                                                                                           >  1:639632/1‑131 (MQ=255)
                                                                                                                                                         tACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCAc                                                                                                                                                                                                 >  1:892938/1‑189 (MQ=255)
                                                                                                                                                          aCCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTg                                                                                                                                                                                                                                                          <  2:639631/131‑1 (MQ=255)
                                                                                                                                                          aCCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACg                                                                                                                                                                                                <  2:892937/189‑1 (MQ=255)
                                                                                                                                                                         gcaCTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCCGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCCGCCTTTCTGTTTTTTAACTCCGCCGGTGAAAATttttttta                                                                                              >  2:1084479/1‑270 (MQ=255)
                                                                                                                                                                                  aCTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTtcat                                                                                                                                                                                                                   >  2:446416/1‑146 (MQ=255)
                                                                                                                                                                                   cTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTtcatc                                                                                                                                                                                                                  <  1:446417/146‑1 (MQ=255)
                                                                                                                                                                                         cgccaaaaaGTTAGATAAAATCGTGCTTTTTCTTGTGTTTTCCACCATATCCTACCTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCtg                                                                              <  2:42590/266‑1 (MQ=255)
                                                                                                                                                                                          aGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTc                                                                                                                                                                                                              >  1:320954/1‑143 (MQ=255)
                                                                                                                                                                                           gCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTca                                                                                                                                                                                                             <  2:320953/143‑1 (MQ=255)
                                                                                                                                                                                            cacaAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGt                                                                                                                                                                                                                        >  2:1105959/1‑131 (MQ=255)
                                                                                                                                                                                            cacaAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTAtttt                                                                                                                                                                             >  2:602779/1‑174 (MQ=255)
                                                                                                                                                                                            cacaAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCtggtg                                                                           <  1:1084480/272‑1 (MQ=255)
                                                                                                                                                                                             acaAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGtt                                                                                                                                                                                                                       <  1:1105960/131‑1 (MQ=255)
                                                                                                                                                                                             acaAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTg                                                                                                                                                                            <  1:602780/174‑1 (MQ=255)
                                                                                                                                                                                                         gAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCACGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTAcaca                                                               >  1:55809/1‑271 (MQ=255)
                                                                                                                                                                                                         gAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGt                                                                                                                                                         >  1:430625/1‑181 (MQ=255)
                                                                                                                                                                                                         gAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTAcaca                                                               >  1:410875/1‑271 (MQ=255)
                                                                                                                                                                                                         gAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTCCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTAcaca                                                               >  1:834711/1‑271 (MQ=255)
                                                                                                                                                                                                         gAAATAGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGACTGGAAGCTCGTATCGCAACCTGAGATACGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGAATTATGGTGGATCTTTAcaca                                                               >  1:750136/1‑271 (MQ=255)
                                                                                                                                                                                                          aaaTCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGt                                                                                                                                                                                                                        >  2:212083/1‑117 (MQ=255)
                                                                                                                                                                                                          aaaTCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGtc                                                                                                                                                        <  2:430624/181‑1 (MQ=255)
                                                                                                                                                                                                           aaTTGTGCTTTTTTCTGTTGTTTCGAACCTATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATTCGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTAcaca                                                               <  1:377563/269‑1 (MQ=255)
                                                                                                                                                                                                            tttgttttttttCCTGTTTTTTCTAACTTATCCTAACTGTCCATGGCGAAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTc                                                           <  2:410874/265‑1 (MQ=255)
                                                                                                                                                                                                                tGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTtcat                                                                                                                                                                                                                   >  1:751876/1‑116 (MQ=255)
                                                                                                                                                                                                                 gCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTtcatc                                                                                                                                                                                                                  <  2:751875/116‑1 (MQ=255)
                                                                                                                                                                                                                    ttttGCTGTTTTTTCAAACATTTCATAATTTCCATTGGCGCAATTACCCGGTAAAATACGCCGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTg                                                   <  2:55809/272‑1 (MQ=255)
                                                                                                                                                                                                                                  cGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTtcat                                                                                                                                                                                                                   >  2:915241/1‑98 (MQ=255)
                                                                                                                                                                                                                                   gAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTtcatc                                                                                                                                                                                                                  <  1:915242/98‑1 (MQ=255)
                                                                                                                                                                                                                                                  gTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGa                     >  1:56221/1‑272 (MQ=255)
                                                                                                                                                                                                                                                  gTCAATTGCGACAGTAACCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGGGCTCATCCTGATTATGGTTCCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGa                     <  2:834710/272‑1 (MQ=255)
                                                                                                                                                                                                                                                   tCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTAAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATcc                                              <  1:1129192/246‑1 (MQ=255)
                                                                                                                                                                                                                                                    ccATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTtcat                                                                                                                                                                                                                   >  2:1114296/1‑80 (MQ=255)
                                                                                                                                                                                                                                                     cATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTtcatc                                                                                                                                                                                                                  <  1:1114297/80‑1 (MQ=255)
                                                                                                                                                                                                                                                     cATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGTTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTcc                        <  1:1057989/266‑1 (MQ=255)
                                                                                                                                                                                                                                                      aTTGCGCAATTACCCGGTAAAATGCGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGCTTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCGGGTGGCTCTTTACCCAGCTCGATGTTTGGCACCGGACACCCAAGAGTTACgggtctggcggg                 >  2:1098405/1‑261 (MQ=255)
                                                                                                                                                                                                                                                      aTTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCaa                                                                                                         >  2:759224/1‑184 (MQ=255)
                                                                                                                                                                                                                                                      aTTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATgggg                                                                                         >  2:582173/1‑200 (MQ=255)
                                                                                                                                                                                                                                                       ttGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAat                                                                                                        <  1:759225/184‑1 (MQ=255)
                                                                                                                                                                                                                                                       ttGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATggggg                                                                                        <  1:582174/200‑1 (MQ=255)
                                                                                                                                                                                                                                                        tttcgcAATTACCCGGTAAAATGCGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTAGTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTcc                        <  1:1098406/261‑1 (MQ=255)
                                                                                                                                                                                                                                                        tGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGCTGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCTCGCACCCCATGCGTTTCGTTCCGGACgggtc               >  2:236333/1‑270 (MQ=255)
                                                                                                                                                                                                                                                             cAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGt                                                                                                                                                         >  1:520904/1‑130 (MQ=255)
                                                                                                                                                                                                                                                              aaTTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGtc                                                                                                                                                        <  2:520903/130‑1 (MQ=255)
                                                                                                                                                                                                                                                              aaTTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACCCAGCTCGATGTCTgg                                                  >  2:879716/1‑232 (MQ=255)
                                                                                                                                                                                                                                                               aTTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGa                                                 <  1:879717/232‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTAtttt                                                                                                                                                                             >  1:475400/1‑102 (MQ=255)
                                                                                                                                                                                                                                                                      tttAAAAACCGCAGAATTTTCCGGGGTTTGGTCAAAAGGGCCTATTCTGGATTATGGTTTAATCATTAACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGt  <  2:56221/270‑1 (MQ=255)
                                                                                                                                                                                                                                                                      ggTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTg                                                                                                                                                                            <  2:475399/102‑1 (MQ=255)
                                                                                                                                                                                                                                                                      ggTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGt  <  1:236333/272‑1 (MQ=255)
                                                                                                                                                                                                                                                                           aaTACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTAcc                                                                                                                     >  2:316150/1‑152 (MQ=255)
                                                                                                                                                                                                                                                                            aTACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCg                                                                                                                    <  1:316151/152‑1 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                       
GGCTGGTCCAGAAACCATGTACCAGTAAAGCGATTATGGCGATCGCGCCAACAATGATTAATATCAGACGCAAATCCTGCATCATTATATTCTCTGTTGTTCTAACACCTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGT  >  NC_000913/2531149‑2531681

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