Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,635,114 A→C L66R (CTA→CGA)  rob ← right oriC‑binding transcriptional activator, AraC family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,635,1140AC96.7% 179.9 / ‑5.9 62L66R (CTA→CGA) robright oriC‑binding transcriptional activator, AraC family
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base C (35/24);  minor base G (0/2);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAG  >  NC_000913/4634966‑4635259
                                                                                                                                                    |                                                                                                                                                 
gcgAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGcgtcgt                                                                                                                                                 <  1:435544/151‑1 (MQ=255)
 cgAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTTAAGTGAAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGGAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTg                                                                                                                                                >  1:261869/1‑151 (MQ=255)
 cgAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTg                                                                                                                                                >  2:740280/1‑151 (MQ=255)
      aCCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGAcc                                                                                                                                           <  2:789259/151‑1 (MQ=255)
        cAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGAccgc                                                                                                                                         >  2:855645/1‑151 (MQ=255)
          aaGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGAccgccg                                                                                                                                       >  2:290964/1‑151 (MQ=255)
                 tCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCg                                                                                                                                 >  1:856002/1‑150 (MQ=255)
                  ccATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGAc                                                                                                                               <  2:343096/151‑1 (MQ=255)
                   cATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGa                                                                                                                                >  1:802565/1‑149 (MQ=255)
                    aTTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACaa                                                                                                                             <  2:942805/151‑1 (MQ=255)
                     ttCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGAAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACaa                                                                                                                             <  2:327371/150‑1 (MQ=255)
                       cAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAAGGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAAAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAacg                                                                                                                          <  2:591448/151‑1 (MQ=255)
                        aGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGGGCGACCGCCGATTTCGACAAacg                                                                                                                          >  2:167251/1‑150 (MQ=255)
                        aGGAGAACGGCGGTAAAGTGCAGGAGTCGGGGCAAACGGCTTCTTGAATGCGGGGGTAAATGTCTGTTGAGAGTCGAAGCGTTATTGCGGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAacga                                                                                                                         <  1:767052/151‑1 (MQ=255)
                          gagaACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAacgacg                                                                                                                       <  1:887078/151‑1 (MQ=255)
                            gaACGGCGGTAAAGTGCAGGGGGCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAacg                                                                                                                          >  1:362684/1‑146 (MQ=255)
                                      aaaGTGCAGGAGTCTGGCCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAGGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAAtata                                                                                                           <  2:64610/151‑1 (MQ=255)
                                            cAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGcc                                                                                                     >  1:615486/1‑151 (MQ=255)
                                             aGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGcc                                                                                                     >  1:58083/1‑150 (MQ=255)
                                             aGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGcc                                                                                                     >  2:718488/1‑150 (MQ=255)
                                                aGTCTGGGCAAACTGTTTCTTGAAGGCGGGGGTAAAGTTCTGTTGGGAGTGGAGGCGGTATTGCAGCGCGATGTCCAGAACCGGCCGCGAAGTCAGGCTGCTTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATa                                                                                                 <  1:39659/151‑1 (MQ=255)
                                                  tCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGc                                                                                               >  2:194341/1‑151 (MQ=255)
                                                    tGGGCAAACTGCTTTTTGAATGCGCGGGTAATTGTCTGTTGAGAGTCGAAGCGTTATTGCAGCGCTATGTCCAGAATCGGCCGCGCATTCAGGCGTCGTGCCACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAg                                                                                                <  1:408552/148‑1 (MQ=255)
                                                     gggCAAACTGCTTCTTGAAGGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATGGCAGCGCGATGTCCAAAACCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATg                                                                                            <  1:559055/151‑1 (MQ=255)
                                                         aaaCTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGcc                                                                                         >  1:174980/1‑150 (MQ=255)
                                                            cTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGGGGGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGt                                                                                      <  2:859115/150‑1 (MQ=255)
                                                             tGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTg                                                                                     >  2:867/1‑150 (MQ=255)
                                                              gCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGAc                                                                                   >  1:465210/1‑151 (MQ=255)
                                                              gCTTCTGAAATGCGCGGGAAAATTTCTGTGGAGAGCGGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCGGCCAGGGTTCGTGCTACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGa                                                                                    <  1:30321/150‑1 (MQ=255)
                                                                 tcttAAATGCGCGGGTAAATGTCTGTTGAGAGTCAAAGCGTTATTGCAGCGCGATGTCCAGATTCGGACGCGCAGTCAGGCGTGGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGa                                                                                    <  2:427429/147‑1 (MQ=255)
                                                                  cttGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCt                                                                               >  2:91657/1‑151 (MQ=255)
                                                                     gAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCttt                                                                             >  1:781070/1‑150 (MQ=255)
                                                                      aaTGCGCGGGTAAATGTCTGTTGAGAGTCGAGGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTaa                                                                           <  2:242369/151‑1 (MQ=255)
                                                                      aaTGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCttt                                                                             <  2:796075/149‑1 (MQ=255)
                                                                           gcgGGTAAATGTCTGTTGAGAGTCGAAGCGGTATGGCAGCGCGATGTCCAGAACCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACAt                                                                       <  1:39007/150‑1 (MQ=255)
                                                                              ggTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATtctc                                                                   >  2:190593/1‑151 (MQ=255)
                                                                                tAAAGGTTTGTTGAGAGTGGAGGCGGTTTTGCAGCGGGATCTCCAGAATCGGCCGGGCAGTAGGGCGTAGTGCGGCCGCCGATGTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATc                                                                                <  2:89400/136‑1 (MQ=255)
                                                                                 aaaTGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTg                                                                 <  2:199952/150‑1 (MQ=255)
                                                                                     gTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGt                                                            >  1:531524/1‑151 (MQ=255)
                                                                                      tCTGTTGAGAGTCGAGGCGGTATTGCAGCGCGATGTCAAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGt                                                            <  2:649794/150‑1 (MQ=255)
                                                                                        tGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGcc                                                         >  2:34678/1‑151 (MQ=255)
                                                                                           tGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACt                                                      >  2:615086/1‑151 (MQ=255)
                                                                                            gagaGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACt                                                      <  1:292125/150‑1 (MQ=255)
                                                                                             agagTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTg                                                    >  2:228508/1‑151 (MQ=255)
                                                                                             agagTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTg                                                    >  1:771282/1‑151 (MQ=255)
                                                                                                    aaGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAAc                                             >  2:243364/1‑151 (MQ=255)
                                                                                                        ggTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTg                                          <  2:14660/150‑1 (MQ=255)
                                                                                                              gCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTcgc                                   >  1:882798/1‑151 (MQ=255)
                                                                                                                  cgcgATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCt                               >  1:112267/1‑151 (MQ=255)
                                                                                                                     gATGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTAc                             >  2:612418/1‑150 (MQ=255)
                                                                                                                       tGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACAtt                          >  1:35848/1‑151 (MQ=255)
                                                                                                                       tGTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACAt                           <  2:39884/150‑1 (MQ=255)
                                                                                                                        gTCCAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTg                         >  1:934096/1‑151 (MQ=255)
                                                                                                                          caagaaggggAGGCGCAGTCGGGCGTGGGTCAACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTg                         <  1:913181/141‑1 (MQ=255)
                                                                                                                           cAGAATCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCg                      >  2:853161/1‑151 (MQ=255)
                                                                                                                              aaTCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGc                   <  2:849366/151‑1 (MQ=255)
                                                                                                                              aaTCGGACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGc                   >  1:735253/1‑151 (MQ=255)
                                                                                                                                  ggAGGCGCAGGCTGGCGTCGGGCCCGCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGAc                <  1:356306/150‑1 (MQ=255)
                                                                                                                                  ggACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGAc                >  1:229418/1‑150 (MQ=255)
                                                                                                                                  ggACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGAc                >  1:5112/1‑150 (MQ=255)
                                                                                                                                   gACGCGCAGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGAc                >  2:675435/1‑149 (MQ=255)
                                                                                                                                          aGTCAGGCGTCGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTg        >  2:567952/1‑150 (MQ=255)
                                                                                                                                               ggCGTCGGGCGCCCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAg  <  1:798359/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAG  >  NC_000913/4634966‑4635259

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑