Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,596,878:1 +G 7.5% intergenic (‑141/+27) ypfJ ← / ← purC putative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,596,8781.G 7.5% 348.0 / 26.4 93intergenic (‑141/+27)ypfJ/purCputative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase
Reads supporting (aligned to +/- strand):  ref base . (39/47);  new base G (4/3);  total (43/50)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

CCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAGGCGGC  >  NC_000913/2596729‑2597026
                                                                                                                                                      |                                                                                                                                                    
ccAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTTCCGTGCCAGCt                                                                                                                                                       >  1:411853/1‑150 (MQ=255)
ccAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCt                                                                                                                                                       >  1:897428/1‑150 (MQ=255)
ccAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCt                                                                                                                                                       <  1:1544308/150‑1 (MQ=255)
    cGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCAc                                                                                                                                                  >  2:959536/1‑150 (MQ=255)
     gCATAGATACCTCAACAATATATTCATTATAGAAATCATCGTAGGTGGTGCAGAGGAAGATTACCACAGGAAAGCGGGTTGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACg                                                                                                                                                 <  1:1094836/150‑1 (MQ=255)
      cATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACgg                                                                                                                                                <  2:205632/150‑1 (MQ=255)
       aTAGATAACTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGc                                                                                                                                               <  2:283187/150‑1 (MQ=255)
         aGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCaa                                                                                                                                             <  1:1212036/150‑1 (MQ=255)
           aTACCTCAACAATATATTCATTATAGAAATGATCGAACGTGGGTCAGAGGAAGATTACCACAGCAAAGCGGGATGAAAAGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCAAga                                                                                                                                           <  1:1051880/150‑1 (MQ=255)
              ccTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTTCCCTGCCAGCT‑GCACGGCAAgatgt                                                                                                                                        >  1:610988/1‑149 (MQ=255)
              ccTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgatga                                                                                                                                        >  2:58279/1‑150 (MQ=255)
              ccTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgatga                                                                                                                                        >  1:1367345/1‑150 (MQ=255)
                   acaaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTTCCGTGCCAGCT‑GAACGGCAAGATGATTaaa                                                                                                                                   >  1:394205/1‑149 (MQ=255)
                     aaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAg                                                                                                                                 <  2:558432/150‑1 (MQ=255)
                        atatTCATTATAGAAATGATAGTAGGTGGTTCAGAGGAAGATTAACACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATAGTCAGAGAGCGGTTGTCTTTCCGTGCAAGCT‑GCAGGGAAAGATGATGAACAGaaa                                                                                                                              >  2:84569/1‑150 (MQ=255)
                         tatTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑g                                                                                                                                                     <  2:288498/126‑1 (MQ=255)
                         tatTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑g                                                                                                                                                     >  1:288498/1‑126 (MQ=255)
                          attaattaTAGAAATGATCTTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGTAATGAGCGTTAAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGaaaaa                                                                                                                            <  1:1392892/150‑1 (MQ=255)
                          atTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGaaaaa                                                                                                                            >  2:915846/1‑150 (MQ=255)
                            tCATTATAGAAATGATCTTAGGTGGTTCAGAGGAAGATTACCACAGGAAATCGTGCTGAAATGGCGGAGTTATGAGCGTTCAATGCATAATATCCGAGAGCGGTTGTCTTGCCGGGCCAGCT‑GCACGGCAAGATGATGAACCGAAAAATc                                                                                                                          <  2:1419708/150‑1 (MQ=255)
                             cATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCAGAGAGCGGTTGTCTTGCCGTGCCCGCT‑GCACGGCAAGATGATGAACAGAAAAATCa                                                                                                                         >  2:25908/1‑150 (MQ=255)
                              aTTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCAAga                                                                                                                                           >  2:966265/1‑131 (MQ=255)
                              aTTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCAAga                                                                                                                                           <  1:966265/131‑1 (MQ=255)
                                tataGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTc                                                                                                                      <  1:233744/150‑1 (MQ=255)
                                  taGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCa                                                                                                                    >  1:1621791/1‑150 (MQ=255)
                                   aGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAg                                                                                                                   <  1:948244/150‑1 (MQ=255)
                                   aGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAg                                                                                                                   <  2:1621791/150‑1 (MQ=255)
                                      aaTGATCGTCGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGGTCAATGCATATATTCCGAGAGCGGGTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTg                                                                                                                <  1:1278701/150‑1 (MQ=255)
                                      aaTGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTg                                                                                                                <  1:1334178/150‑1 (MQ=255)
                                       aTGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGt                                                                                                               >  1:668595/1‑150 (MQ=255)
                                        tGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTa                                                                                                              <  2:668595/150‑1 (MQ=255)
                                         gATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGACAAGATGATGAACAGAAAAATCAGTCCAGCTGTac                                                                                                             >  1:411802/1‑150 (MQ=255)
                                            cGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATTCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACAcc                                                                                                          <  1:98052/150‑1 (MQ=255)
                                              tAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCa                                                                                                        <  2:795933/150‑1 (MQ=255)
                                                gggggtTAAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGCAAAATCAGTCACGCGGTTCACCCacc                                                                                                      >  2:738247/4‑148 (MQ=255)
                                                   ggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTTCCGTGCCAGCT‑GCCCGGCAATATGATTAACAGAAAAATCAGTCCCGCTGTACACCCAggcgg                                                                                                   >  1:798499/1‑150 (MQ=255)
                                                   ggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGAACACCCAggcgg                                                                                                   <  2:1441789/150‑1 (MQ=255)
                                                     ttCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCAAGGCAAGATGATGAACAGAAAAATCAGTCAAGCTGTACACCCAGgcggcg                                                                                                 >  2:88362/1‑150 (MQ=255)
                                                       cAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcggg                                                                                               <  1:1290103/150‑1 (MQ=255)
                                                        agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTTCCGTGCCAGCT‑GCACGGCAAGATGATGAACACAAAAAACAGTCCAGCTGTACACCCAGGCggcgggc                                                                                              >  2:516127/1‑150 (MQ=255)
                                                        agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcgggc                                                                                              <  2:236761/150‑1 (MQ=255)
                                                         gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcgggcg                                                                                             >  2:738129/1‑150 (MQ=255)
                                                         gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcgggcg                                                                                             >  2:1228967/1‑150 (MQ=255)
                                                         gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACCGCAAGATGATGAACAGAAAAATCAGTCCCGCTGGACACCCCGGCggcgggcg                                                                                             >  1:1433218/1‑150 (MQ=255)
                                                            gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCACCT‑GCACCGCAAGATGATGATCAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACg                                                                                          >  2:498011/1‑150 (MQ=255)
                                                             aaGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                                                                         <  1:986904/150‑1 (MQ=255)
                                                                aTTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGGTCAATGCACATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCtt                                                                                      <  2:1518099/150‑1 (MQ=255)
                                                                    ccACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                                                                                  >  1:212213/1‑150 (MQ=255)
                                                                       caGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGc                                                                               <  2:580428/150‑1 (MQ=255)
                                                                         ggAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GGACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCtt                                                                             >  1:1547017/1‑150 (MQ=255)
                                                                         ggAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCtt                                                                             >  1:988065/1‑150 (MQ=255)
                                                                           aaaGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGCTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCg                                                                           <  1:1347883/150‑1 (MQ=255)
                                                                                gggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAg                                                                      <  1:32499/150‑1 (MQ=255)
                                                                                  gATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGc                                                                    <  2:897428/150‑1 (MQ=255)
                                                                                     gAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAc                                                                 >  2:207019/1‑150 (MQ=255)
                                                                                      aaaTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCATCT‑GCAAGGCCAGATGATGAACAGAAAAATCAGTCAAGCTGTACACCCACGCGGCTGGCGACGGCTTCATAGGCGAAGATCAGTCCAcc                                                                >  1:326099/1‑150 (MQ=255)
                                                                                             ggaggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCAAGCT‑ACACTGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGctc                                                         <  2:268172/150‑1 (MQ=255)
                                                                                                ggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTgg                                                      <  2:518397/150‑1 (MQ=255)
                                                                                                ggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTgg                                                      <  2:1367345/150‑1 (MQ=255)
                                                                                                 caaacaGCGCTCAATGCATATATTCCGAGAGCGGCTGCCTTGCCGTGCCAGCG‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTGCATAGGCTTCGATCAGGCCACCGAGGCTCTGGc                                                     <  1:84569/145‑1 (MQ=255)
                                                                                                  aaTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCg                                                    <  1:402631/150‑1 (MQ=255)
                                                                                                  aaTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCg                                                    <  1:396480/150‑1 (MQ=255)
                                                                                                        cGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCCCGGCAAGATGATGAACAGAAAAATCAGACCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGAACAGGCCACCCCGGCTCTGGCGGAAAcc                                              >  1:1626928/1‑149 (MQ=255)
                                                                                                            cAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAAACAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTGATAGGCTTAGATCAGGCCACCGACGCTCGGGCGGAAACGGTCt                                          >  1:704566/1‑150 (MQ=255)
                                                                                                             aaTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGTCAAGATGATGAACAGATAAATCAGTCCAGCTGTACACCCTGGCGGCTGGCGACGGCTTCATAGGCGTAGATCAGGCCACCGACGCTCTGCCGGAATCGGTCtt                                         >  1:265461/1‑150 (MQ=255)
                                                                                                             aaTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCtt                                         >  2:67223/1‑150 (MQ=255)
                                                                                                             aaTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCtt                                         >  2:580999/1‑150 (MQ=255)
                                                                                                             aaTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCtt                                         <  2:514793/150‑1 (MQ=255)
                                                                                                             aaTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCtt                                         <  1:1032347/150‑1 (MQ=255)
                                                                                                             aaTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCtt                                         <  1:1296709/150‑1 (MQ=255)
                                                                                                                 cATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTc                                     <  1:1261271/150‑1 (MQ=255)
                                                                                                                 cATATATTCCGAGAGCGGTTGTCATTCCGTGCCAGCA‑GCACGGCAAGATGATTAACACAAAAATCAGTACAGCTCTAAACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCCCCGAGGCACAGGCGGATACGGTCTTTGtt                                     >  2:1400539/1‑149 (MQ=255)
                                                                                                                  atatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCTAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTcc                                    <  1:959536/150‑1 (MQ=255)
                                                                                                                        tCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTtc                              <  2:20374/150‑1 (MQ=255)
                                                                                                                         ccGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTtct                             >  1:1110107/1‑150 (MQ=255)
                                                                                                                          cGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTtctc                            <  1:1228967/150‑1 (MQ=255)
                                                                                                                          cGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTtctc                            <  1:738129/150‑1 (MQ=255)
                                                                                                                              cgCGGGTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGc                        <  2:1433218/149‑1 (MQ=255)
                                                                                                                               gCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCg                       <  2:798499/150‑1 (MQ=255)
                                                                                                                                 ggTTGTCTTGCCGTGCCAGCT‑GCTCGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTCCCCCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCGGGCGGAAACGGTCTTTGTCCATTTTCTGCAGCGtt                     >  1:890548/1‑150 (MQ=255)
                                                                                                                                 ggTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAACCGGTCTTTGTCCATTTTCTCCCGCGtt                     >  2:112414/1‑150 (MQ=255)
                                                                                                                                 ggTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCCTTTTCTCCAGCggt                     >  1:131596/1‑148 (MQ=255)
                                                                                                                                 ggTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGtt                     >  2:1023028/1‑150 (MQ=255)
                                                                                                                                    tGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCt                  <  2:394205/150‑1 (MQ=255)
                                                                                                                                     gTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCGCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCtt                 <  1:34159/150‑1 (MQ=255)
                                                                                                                                      tCTTGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCttt                >  1:1139056/1‑150 (MQ=255)
                                                                                                                                       cTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTg               >  1:1033022/1‑150 (MQ=255)
                                                                                                                                        ttGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGt              <  2:988065/150‑1 (MQ=255)
                                                                                                                                        ttGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGt              >  2:1322746/1‑150 (MQ=255)
                                                                                                                                        ttGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGt              <  1:1322746/150‑1 (MQ=255)
                                                                                                                                         tGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAAGCCCCCGCGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCCGCGTTTCTTTGTc             >  2:1564393/1‑150 (MQ=255)
                                                                                                                                          gCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTcc            >  2:449166/1‑150 (MQ=255)
                                                                                                                                          gCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTcc            <  1:983744/150‑1 (MQ=255)
                                                                                                                                          gCCGTGCCAGCT‑GCACGACAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTcc            <  2:411802/150‑1 (MQ=255)
                                                                                                                                           ccGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTccc           >  2:1564177/1‑150 (MQ=255)
                                                                                                                                           ccGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTccc           <  1:1430124/150‑1 (MQ=255)
                                                                                                                                               gCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCACCGTTTCTTTGTCCCACAg       >  1:300745/1‑150 (MQ=255)
                                                                                                                                                    cT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAggcggc  <  1:1387489/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAGGCGGC  >  NC_000913/2596729‑2597026

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