Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,788,507 T→G 10.1% Y25D (TAT→GAT)  aroH → 3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,788,5070TG10.1% 150.2 / 9.0 60Y25D (TAT→GAT) aroH3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible
Reads supporting (aligned to +/- strand):  ref base T (22/31);  new base G (4/2);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

CAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGTACG  >  NC_000913/1788359‑1788648
                                                                                                                                                    |                                                                                                                                             
cagcGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTa                                                                                                                                              <  2:16556/150‑1 (MQ=255)
      ttGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTGCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGt                                                                                                                                        <  2:354081/150‑1 (MQ=255)
         gTTTATCGTGATGCTCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGTACTGACGAACACCGTACAGCGCGTATTGAGTGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAAc                                                                                                                                     <  2:161935/150‑1 (MQ=255)
            tATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGcc                                                                                                                                  >  2:161136/1‑150 (MQ=255)
            tATCGGGCTGCGCAGCACGTCCCGGCAGTCCTGCCGGAGAAGCACCAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTCACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGcc                                                                                                                                  <  1:103646/150‑1 (MQ=255)
             aTCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCt                                                                                                                                 >  2:172924/1‑150 (MQ=255)
             aTCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCt                                                                                                                                 >  1:148520/1‑150 (MQ=255)
              tCGTGATGTGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGTTACGGGATACCGTAgcgcgag                                                                                                                                <  2:1060252/150‑4 (MQ=255)
                 ggATGCGCATAACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGTTACGGGATACCGTAgcgcgagttc                                                                                                                             <  2:200112/149‑7 (MQ=255)
                       gcATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCa                                                                                                                       <  2:348212/150‑1 (MQ=255)
                        cATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGcccc                                                                                                                      >  1:675909/1‑148 (MQ=255)
                         aTCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCAcc                                                                                                                     <  1:521015/150‑1 (MQ=255)
                                cccGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCa                                                                                                              <  2:1115782/150‑1 (MQ=255)
                                          cTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGGATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCTCCCCCCATGTCCCCGCCGCccg                                                                                                    >  2:1170841/1‑150 (MQ=255)
                                                  gAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGaa                                                                                            >  1:479126/1‑150 (MQ=255)
                                                    aGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAAtt                                                                                          >  1:755663/1‑150 (MQ=255)
                                                    aGCAACAAATTTCTGAGACTTGAAATGAACAGAACTTACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAAtt                                                                                          <  2:845125/150‑1 (MQ=255)
                                                      caacaaATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCAGCAGAATTGa                                                                                        >  1:181565/1‑150 (MQ=255)
                                                       aacaaATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGaa                                                                                       <  2:1010826/150‑1 (MQ=255)
                                                        acaaATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGaaa                                                                                      >  1:585302/1‑150 (MQ=255)
                                                                 tGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGa                                                                             >  1:1275392/1‑150 (MQ=255)
                                                                     aCTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATgg                                                                         <  1:1044876/150‑1 (MQ=255)
                                                                          tAATGAACAGAACGGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAg                                                                    <  1:809925/150‑1 (MQ=255)
                                                                               aaCAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAg                                                               <  1:572722/150‑1 (MQ=255)
                                                                                aCAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGc                                                              <  2:14701/150‑1 (MQ=255)
                                                                                  aGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGa                                                            <  1:265672/150‑1 (MQ=255)
                                                                                        gACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTg                                                      >  1:727368/1‑150 (MQ=255)
                                                                                         aCGACCTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCGGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTTAAAAAATACTGAATGGTGAAGATAAGCGACTGTTgg                                                     <  1:509494/150‑1 (MQ=255)
                                                                                          cGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGt                                                    <  1:161136/150‑1 (MQ=255)
                                                                                            aaCTCCGTACTGCGCGTAGTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTAGCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCa                                                  <  2:146121/150‑1 (MQ=255)
                                                                                             aCTCCGACCTGCGCTGATTGAGAGCCTGGTACCGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTCAAGATAAGCGACTGTTGGTCat                                                 <  2:1135991/150‑1 (MQ=255)
                                                                                                  gTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGTTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCTCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTg                                            <  1:72702/150‑1 (MQ=255)
                                                                                                   tACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTgg                                           <  1:493473/150‑1 (MQ=255)
                                                                                                   tACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTgg                                           <  2:524726/150‑1 (MQ=255)
                                                                                                    aCTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACACCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGc                                          <  1:1040885/150‑1 (MQ=255)
                                                                                                       gcgcGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCCCCCATGTCCCCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGcccc                                       >  1:713566/1‑150 (MQ=255)
                                                                                                           gTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGcc                                         <  2:267987/144‑1 (MQ=255)
                                                                                                           gTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGcc                                         >  1:267987/1‑144 (MQ=255)
                                                                                                            tATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAACTACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCACTGCTc                                  >  1:892086/1‑150 (MQ=255)
                                                                                                            ccaagagaGCCCTGGAACGCCCTCCCAACTCTCGCTACGGAATCCCGTAACGCCGGGCGTCGCCCCCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTCAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTc                                  <  1:70707/146‑1 (MQ=255)
                                                                                                             aTTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACAGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCg                                 <  2:516591/150‑1 (MQ=255)
                                                                                                                gagaGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATc                              >  2:600613/1‑150 (MQ=255)
                                                                                                                  gagCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCa                            >  1:288782/1‑150 (MQ=255)
                                                                                                                    gCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACg                          >  2:423256/1‑150 (MQ=255)
                                                                                                                    gCCGGGTAACGCCCGCCGAACTCGGGCTACGGTATCCCGTAACGCCTGGCGGCGCCACCCATGTCACCGACTCCCGCCGCATAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACg                          <  2:969497/150‑1 (MQ=255)
                                                                                                                    cccTGCTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGGGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACg                          <  2:1001393/149‑1 (MQ=255)
                                                                                                                        ggTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGAtct                      >  1:613915/1‑150 (MQ=255)
                                                                                                                         ttAACGCCCGCCGAACTCGCGCTACGTTATCCCGAAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAATTACTCAATGGTTAATATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGAtctc                     <  1:1125235/149‑1 (MQ=255)
                                                                                                                             cgcccgccGAACTCGCGCTACGGTATCCCGTAACGCCTGGCTTCGCCACCCATGTCACCGACTCCCGCCGCATAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCg                 <  2:956480/150‑1 (MQ=255)
                                                                                                                             cgcccgccGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCg                 >  2:1262479/1‑150 (MQ=255)
                                                                                                                                 cgccgAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGc             >  1:756126/1‑150 (MQ=255)
                                                                                                                                 cgccgAACTCGCGCTACGGTATCCCGTAACGCCTGGAGTCGCCACCCATGTCACAGACTCCCGCCGCAGAAATGAAAAAATACTGAATCGTGAAGATAAGCGACTGGTGCACATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGc             >  1:195123/1‑150 (MQ=255)
                                                                                                                                     gAACTCGCGCTACGGTATCCCATAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATg         >  2:158220/1‑150 (MQ=255)
                                                                                                                                      caCTGGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACGCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCGCGATCCACGATCTCACCGCTGCAATgg        <  2:224255/149‑1 (MQ=255)
                                                                                                                                      aaCTCGCGCTACGGTATCCCGTAACGCCGGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATgg        <  1:1209686/150‑1 (MQ=255)
                                                                                                                                       ccaggcgtTACGGGATACCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATACGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGa       >  1:1004438/5‑150 (MQ=255)
                                                                                                                                        caggcgtTACGGGATACCGTAGCGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAg      >  1:333126/4‑150 (MQ=255)
                                                                                                                                        cTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAg      <  2:11491/150‑1 (MQ=255)
                                                                                                                                           gcgtTACGGGATACCGTAGCGCCTGGCGTCGCCACCAATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGAGAACATAAGAGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGTAc   >  2:42162/1‑150 (MQ=255)
                                                                                                                                            cgttACGGGATACCGTAGCGCCTGGCGTCGCCACCCATGTCACCGCCTCCCGCCGCAGAAGTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGTACg  >  2:798934/4‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
CAGCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGTACG  >  NC_000913/1788359‑1788648

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