Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,592,990 A→C 8.2% intergenic (+28/+82) ypfN → / ← ypfH UPF0370 protein YpfN/esterase YpfH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,592,9900AC 8.2% 210.8 / 11.0 73intergenic (+28/+82)ypfN/ypfHUPF0370 protein YpfN/esterase YpfH
Reads supporting (aligned to +/- strand):  ref base A (27/40);  new base C (3/3);  total (30/43)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

CAAGCGCGTCGACCATAAGAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTT  >  NC_000913/2592842‑2593139
                                                                                                                                                    |                                                                                                                                                     
cAAGCGCGTCGACCATAAGAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAg                                                                                                                                                      <  1:981738/150‑1 (MQ=255)
cAAGCGCGTCGACCATAAGAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAg                                                                                                                                                      <  2:338487/150‑1 (MQ=255)
cAAGCGCGTCGACCATAAGAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAg                                                                                                                                                      >  1:1358020/1‑150 (MQ=255)
   gcgcGTCGACCATAAGAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCg                                                                                                                                                   <  2:1004042/150‑1 (MQ=255)
     gcgTCGACCATAAGAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGtt                                                                                                                                                 <  1:619250/150‑1 (MQ=255)
      ggTCGACCATAAGAAATTTCTCGCCTACAAGCCGGAGCTTGCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTa                                                                                                                                                <  1:597606/149‑1 (MQ=255)
      cgTCGACCATAAGAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTa                                                                                                                                                <  2:1041908/150‑1 (MQ=255)
       gTCGACCATAAGAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATGCGGCGTTaa                                                                                                                                               <  1:1379396/150‑1 (MQ=255)
       gTCGACCATAAGAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTaa                                                                                                                                               >  2:752490/1‑150 (MQ=255)
         cGACCATAAGAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAAcg                                                                                                                                             <  1:230240/150‑1 (MQ=255)
                aaGAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATGCGGCGTTAAcgccttatc                                                                                                                                      <  1:1314748/150‑6 (MQ=255)
                  gAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCg                                                                                                                                    >  1:763764/1‑150 (MQ=255)
                   aaaTTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCgg                                                                                                                                   >  2:1338290/1‑150 (MQ=255)
                    aaTTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                                                  <  1:1001721/150‑1 (MQ=255)
                    aaTTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                                                  <  2:450959/150‑1 (MQ=255)
                     aTTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCa                                                                                                                                 >  1:574296/1‑150 (MQ=255)
                        tctcGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATg                                                                                                                              <  2:1110657/150‑1 (MQ=255)
                         ctcGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATgg                                                                                                                             <  2:242905/150‑1 (MQ=255)
                             ccAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGaa                                                                                                                         >  1:1277761/1‑150 (MQ=255)
                              caacaaGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAAc                                                                                                                        <  2:594045/150‑1 (MQ=255)
                              caacaaGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAAc                                                                                                                        <  2:515614/150‑1 (MQ=255)
                                 caaGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCgg                                                                                                                                   >  2:593282/1‑136 (MQ=255)
                                 caaGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCgg                                                                                                                                   <  1:593282/136‑1 (MQ=255)
                                 caaGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGAt                                                                                                                     <  2:433187/150‑1 (MQ=255)
                                   aGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGc                                                                                                                   >  1:1245217/1‑150 (MQ=255)
                                   aGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGc                                                                                                                   <  2:787777/150‑1 (MQ=255)
                                     ccGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAAGGCCGCATCGGGCAATGGTGAAGGATGCCt                                                                                                                 <  2:1083027/150‑1 (MQ=255)
                                     ccGGAGATTCCCCAGCAACGTGATTTCAACGATAAGTGCGACGATGACGACGACTGGCCGAAAACGGATCAAAGGAAGAAGTAAGCTGCGACCGCAACCGTAACCAGCATAAGGCGTTAACGCAGCAGCAGGCAATGGTGAACGATACCt                                                                                                                 >  2:1248966/1‑150 (MQ=255)
                                          gCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGc                                                                                                            <  1:1072289/150‑1 (MQ=255)
                                          gCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGc                                                                                                            <  1:1037126/150‑1 (MQ=255)
                                             tCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGAc                                                                                                         >  1:189049/1‑150 (MQ=255)
                                              cccccGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACg                                                                                                        <  2:574296/150‑1 (MQ=255)
                                                      tCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCAcgcg                                                                                                <  2:29776/150‑1 (MQ=255)
                                                      tCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCAcgcg                                                                                                <  1:1338290/150‑1 (MQ=255)
                                                        gTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTc                                                                                              <  1:1026839/150‑1 (MQ=255)
                                                         tGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCt                                                                                             <  1:760944/150‑1 (MQ=255)
                                                         tGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCt                                                                                             <  1:746333/150‑1 (MQ=255)
                                                                   gATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATAAGACGCGTGAGCGTCGCATCAGGcat                                                                                   <  2:1358020/150‑3 (MQ=37)
                                                                       aGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTa                                                                               >  1:692936/1‑150 (MQ=255)
                                                                          gggACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATc                                                                            <  2:1245217/150‑1 (MQ=255)
                                                                           ggACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCt                                                                           <  1:752490/150‑1 (MQ=255)
                                                                           ggACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCt                                                                           >  1:1081279/1‑150 (MQ=255)
                                                                                  gacgacgacTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCtttt                                                                    >  2:1192329/1‑150 (MQ=255)
                                                                                  gacgacgacTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCtttt                                                                    >  1:721300/1‑150 (MQ=255)
                                                                                    cgacgacTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGTCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTAc                                                                  <  2:1324475/150‑1 (MQ=255)
                                                                                    cgacgacTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTAc                                                                  <  1:226254/150‑1 (MQ=255)
                                                                                     gacgacTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTCGGCCGGCTAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTAca                                                                 <  2:1395665/150‑1 (MQ=255)
                                                                                          cTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCAtt                                                            >  2:1014673/1‑150 (MQ=255)
                                                                                          cTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCAtt                                                            <  2:710274/150‑1 (MQ=255)
                                                                                            ggCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTc                                                          >  1:924070/1‑150 (MQ=255)
                                                                                            ggCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTc                                                          >  2:346553/1‑150 (MQ=255)
                                                                                              ccGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCaa                                                        <  2:692936/150‑1 (MQ=255)
                                                                                                   aaaGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATcaca                                                   <  1:153312/150‑1 (MQ=255)
                                                                                                        aTCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCAcatcatc                                              <  2:80286/150‑1 (MQ=255)
                                                                                                         tCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCAcatcatca                                             <  2:721300/150‑1 (MQ=255)
                                                                                                           aaCCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACatcatcatc                                           >  2:1273582/1‑150 (MQ=255)
                                                                                                               gaagaaGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCa                                       >  1:459174/1‑150 (MQ=255)
                                                                                                                 agaagTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAgg                                     <  1:1156224/150‑1 (MQ=255)
                                                                                                                         gCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGc                             <  2:238180/150‑1 (MQ=255)
                                                                                                                             aGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACt                         >  2:1127863/1‑150 (MQ=255)
                                                                                                                                gCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAg                      <  2:1066362/150‑1 (MQ=255)
                                                                                                                                gCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAg                      >  2:159368/1‑150 (MQ=255)
                                                                                                                                gCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAg                      <  1:148417/150‑1 (MQ=255)
                                                                                                                                 cTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAgc                     >  2:250956/1‑150 (MQ=255)
                                                                                                                                 cTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAgc                     >  2:900624/1‑150 (MQ=255)
                                                                                                                                      accattgCCGGATGCGGCGTTAACGCCTTATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTtc                <  2:60433/144‑1 (MQ=39)
                                                                                                                                         attgCCGGATGCGGCGTTAACGCCTTATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTtcgtc             >  1:1070655/4‑150 (MQ=255)
                                                                                                                                          aGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTtcgtcg            >  2:384105/1‑150 (MQ=255)
                                                                                                                                           tgacGGATGCGGCGTTAACGCCTTATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGa           >  2:805893/4‑150 (MQ=39)
                                                                                                                                           ggCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGa           >  2:606089/1‑150 (MQ=255)
                                                                                                                                           ggCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGa           >  2:1221113/1‑150 (MQ=255)
                                                                                                                                           ggCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGa           >  1:1236957/1‑150 (MQ=255)
                                                                                                                                            gCCGGATGCGGCGTTAACGCCTTATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGaa          >  1:443278/1‑150 (MQ=255)
                                                                                                                                                    aGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGtt  <  1:1221113/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAAGCGCGTCGACCATAAGAAATTTCTCGCCAACAAGCCGGAGCTTCCCCCGCATCGTGATTTCAACGATAAGTGGGACGATGACGACGACTGGCCGAAAAAGGATCAACCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAATCACATCATCATCGCCAGGCTTACCGCCACTTAGCGCTTCGTCGAAATAATGTT  >  NC_000913/2592842‑2593139

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