Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,640,231 C→A 72.0% E26* (GAA→TAA)  arcA ← DNA‑binding transcriptional dual regulator ArcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,640,2310CA72.0% 84.2 / 31.7 62E26* (GAA→TAA) arcADNA‑binding transcriptional dual regulator ArcA
Reads supporting (aligned to +/- strand):  ref base C (8/10);  new base A (22/21);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 7.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

GCCTGCTCGCGCAGTTCACGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGCTACCTAAATTGCCAACTAAATCGAAACAGGAAGTACAAAAGTCCCTGACCTGCCTGATGCA  >  NC_000913/4640089‑4640373
                                                                                                                                               |                                                                                                                                              
gCCTGCTCGCGCAGTTCACGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAACACATCATAGCCTTc                                                                                                                                                >  2:592217/1‑143 (MQ=255)
 ccTGCTCGCGCAGTTCACGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTcc                                                                                                                                               >  2:46604/1‑143 (MQ=255)
   tGCTCGCGCAGTTCACGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTCCGc                                                                                                                                             <  1:573139/143‑1 (MQ=255)
    gCTCGCGCAGTTCACGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGGCATATTCAGAGAGGATCTGATGCATTTCCGCGCCCTCCGTCGCTTCGCAACCATCATAGCCTTGCGCt                                                                                                                                            >  1:226146/1‑143 (MQ=255)
       cgcgCAGTTCACGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCCATCGCTTCGAAAACATCATAACCTTcacctccg                                                                                                                                         >  2:1107068/1‑136 (MQ=255)
         cgcAGTTCACGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTACGCTTCGaa                                                                                                                                       <  1:398093/143‑1 (MQ=255)
          gcAGTTCACGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTACGCTTCGaaa                                                                                                                                      <  1:195585/143‑1 (MQ=255)
             gTTCACGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTACAAAACCTCATAGCCGCACGCCCCGAAAATa                                                                                                                                   >  2:382002/1‑143 (MQ=255)
              ttCACGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTACGCTTCGAAAATAc                                                                                                                                  >  2:562862/1‑143 (MQ=255)
                cACGCGCTGACAGAAGACCGTGCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTACGCTTCGaatttcgaa                                                                                                                                <  2:161472/143‑8 (MQ=255)
                cACGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTACGCTTCGaatttcgaa                                                                                                                                <  2:699266/143‑8 (MQ=255)
                 aCGCGCTACCAGAAGACCGTTCTCACCCGGCAGATTGATATCCATGATCACCATGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTACGCTTCGAAAATACttt                                                                                                                               <  2:434487/143‑1 (MQ=255)
                 aCGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACttt                                                                                                                               >  2:405747/1‑143 (MQ=255)
                 aCGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACttt                                                                                                                               >  2:413775/1‑143 (MQ=255)
                  cgcgCTAACAGAAGACCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGCTCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATACCCTTCCGCCTCGAAAATACtttt                                                                                                                              >  2:466499/1‑143 (MQ=255)
                         aCAGAAGACCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCCTTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTACGCTTCGAAAATACTTTTCAACgtg                                                                                                                       >  2:551148/1‑143 (MQ=255)
                         aCAGAAGACCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTACGCTTCGAAAATACTTTTCAACgtg                                                                                                                       <  1:932891/143‑1 (MQ=255)
                           agaagaCCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTACGCTTCGAAAATACTTTTCAACGTGtt                                                                                                                     <  1:334061/143‑1 (MQ=255)
                           agaagaCCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCATCCTCTTCGAAAACACTTGTCTACGTGtt                                                                                                                     >  2:265162/1‑143 (MQ=255)
                            gaagaCCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGACGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTCCGCTTCGCAACTACTTTTCAACGTGTTg                                                                                                                    >  1:870869/1‑143 (MQ=255)
                             aagaCCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGc                                                                                                                   <  1:963826/143‑1 (MQ=255)
                             aagaCCGCGCTGACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGc                                                                                                                   <  1:218625/143‑1 (MQ=255)
                                aCCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCCACGGTTTGCggg                                                                                                                >  1:977962/1‑141 (MQ=255)
                                   gTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTACGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTa                                                                                                             >  1:669100/1‑143 (MQ=255)
                                       ttgcccGGCAGATTGATAGCCATGAGGAGCGGGGTGAGGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGGTTCGAAAACATCATAGCCTTACGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCaa                                                                                                         <  1:898368/140‑1 (MQ=255)
                                           ccGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGGCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCg                                                                                                     >  1:819937/1‑143 (MQ=255)
                                            cGGCAGATTGATATCCATGATCGCCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTACGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGt                                                                                                    <  2:762137/143‑1 (MQ=255)
                                              gCAGATTGATAGCCATGATCACCAGGTTGGTGTCATATTCAGAGAGGCTCTGATGCAGTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTGACGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCt                                                                                                  <  2:531649/143‑1 (MQ=255)
                                                  aTTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTGCCCACTCGTCTTCaa                                                                                              >  1:176252/1‑143 (MQ=255)
                                                    ggataggcATGATCACCGGGTTGATGTCGTGTTCAGGGGGGATCTGGTGCATTTCCGCGCCATCTGTCGCTGCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACg                                                                                            <  2:344503/136‑1 (MQ=255)
                                                      atatCCATGATCACCAGGTTGATGTCATATTCAGAGCGGATCGGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTACGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGAt                                                                                          <  1:624214/143‑1 (MQ=255)
                                                         tCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTACGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTGGTCTTCACCGATAAg                                                                                       >  2:950249/1‑143 (MQ=255)
                                                           cATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGaa                                                                                     <  2:471612/143‑1 (MQ=255)
                                                            aTGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTACGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAAt                                                                                    <  2:661218/143‑1 (MQ=255)
                                                              gATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTACGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTTCCATCTTGTCTTCAACGATAAGAAtgt                                                                                  <  2:396104/143‑1 (MQ=255)
                                                                 cACCAGGTTGATGTCATATTCAGAGAGGATCGGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTACGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCg                                                                               <  1:298240/143‑1 (MQ=255)
                                                                   ccAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTACGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCggg                                                                             <  1:880967/143‑1 (MQ=255)
                                                                         tGATGTCATATTCAGAGAGGATCTGATGCAGTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTGCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGc                                                                       <  1:892231/143‑1 (MQ=255)
                                                                             gTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTACGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGt                                                                   >  1:624700/1‑143 (MQ=255)
                                                                             gTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTACGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGt                                                                   >  2:285414/1‑143 (MQ=255)
                                                                               cATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGACAACATCATAGCCTTACGCTTCGAAAATACTTTTCCACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGACTGTGCGGGCTCTGCATGttt                                                                 >  2:15096/1‑143 (MQ=255)
                                                                                atatTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTACGCTTCGAAAATACTTTTCAACGTGTTGTTGGGGACCAACTTGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTg                                                                >  1:754763/1‑143 (MQ=255)
                                                                                 tatTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTACGCTTCGATATTTCTTTTCTACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGc                                                               <  1:36441/143‑1 (MQ=255)
                                                                                          agGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAACCATCATAGCCTTACGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCACCGATAAGACTGTGCGGGGTCCGCATGTTTGCTACCTACAt                                                      >  2:151266/1‑143 (MQ=255)
                                                                                                ttttGCATTTCCGCCTCATCTGTCTCTTCTCGAGTCATATATCGTTCCGCTTCGAAATTACTTTTCAACGTGTTGCGTGTTACCATCTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGCTACCTAAATTGCCaa                                                <  2:778117/140‑1 (MQ=255)
                                                                                                tGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTACGCTTCGAAAATACTTTTCAACGTGATGCGTGTTACCAACTCGTCTTCAACGATAACAATGTGCGGAGTCTGCATGTATGCTACCTAAAATGCCaa                                                >  2:291647/1‑143 (MQ=255)
                                                                                                tGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTACGCTTCGAAAATACTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGCTACCTAAATTGCCAAc                                               >  1:228279/1‑143 (MQ=255)
                                                                                                   tGCATTTCCGAGCCATCTGTCGCTTCGAAGACGTCATAGCCTTCCTCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGCTACCTAAATTGCCAACTa                                             <  2:66365/143‑1 (MQ=255)
                                                                                                    gCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTACGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGCTACCTAAATTGCCAACTaa                                            >  1:476385/1‑143 (MQ=255)
                                                                                                         tCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTACGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGCTACCTAAATTGCCAACTAAATCGa                                       >  1:749193/1‑143 (MQ=255)
                                                                                                           cgcgCCATCTGTCGCTTCGAAAACATCATAGCCTTACGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGCTACCTAAATTGCCAACTAAATCGaaa                                     >  2:316632/1‑143 (MQ=255)
                                                                                                                    tGTCGCTTCGAAAACATCATAGCCTTACGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGCTACCTACATTGCCAACTAAATCGAAAAAGGAAGTa                            >  1:704919/1‑143 (MQ=255)
                                                                                                                     gTCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGCTACCTAAATTGCCAACTAAATCGAAACAGGAAGTAc                           <  1:643978/143‑1 (MQ=255)
                                                                                                                      tCGCTTCGAAAACATCATAGCCTTACGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGCTACCTAAATTGCCAACTAAATCGAAACAGGAAGTACa                          >  2:287105/1‑143 (MQ=255)
                                                                                                                         cTTCGAAAACATCATAGCCTTACGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCCACGATCAGAATTTGCGGGGTCTGCATGTTTGCTACCTAAATTGCCAACTAAATCGAAACAGGAAGTACaaaa                       <  1:416007/143‑1 (MQ=255)
                                                                                                                          ttCGAAAACATCATAGCCTTAGGCTGCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGCTACCTAAATTGACAACTAGATCGAAACAGGAAGTACAAAAg                      >  2:291963/1‑143 (MQ=255)
                                                                                                                          ttCGAAAACATCATAGCCTTACGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAAGGTGCGGGGTCTGCATGTTTGCTACCTAAATTGCCAACTAAATCGAAACAGGAAGTACAAAAg                      >  1:529115/1‑143 (MQ=255)
                                                                                                                           tCGAAAACATCATAGCCTTACGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGCTACCTAAATTGCCAACTAAATCGAAACAGGAAGTCCAAAAGt                     >  1:463817/1‑143 (MQ=255)
                                                                                                                                  ccaccatAGCGTTACGCTGCGAAAATACTTTTCAGCGTGTTGCGTGTTACCACCTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGCTACCTAAACTGCCAACTAAATCGAAACAGGAAGAACAAATGTCCCTGa               <  2:973046/139‑1 (MQ=255)
                                                                                                                                  aCATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGCTACCTAAATTGCCAACTAAATCGAAACAGGAAGTACAAAAGTCCCTGa               <  1:972341/143‑1 (MQ=255)
                                                                                                                                   catcatTGCCTTACGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGCTACCTAAATTGCCAACTAAATCGAAACAGGAAGTACAAAAGTCCCTGAc              <  2:754763/143‑1 (MQ=255)
                                                                                                                                   catcatAGCCTTACGCTTCGAAAATACTTTTCAACGTGTTACGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGCTACCTAAATTGCCAACTCAATCGCAACAGGAAGTACAAAACGCCCTGAc              >  2:903675/1‑143 (MQ=255)
                                                                                                                                       atAGCCTTACGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGCTACCTAAATTGCCAACTAAATCGAAACAGGAAGTACAAAAGTCCCTGAcctgc          <  1:382002/143‑1 (MQ=255)
                                                                                                                                         aGCCTTACGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGCTTGCTACCTAAATTGCCAACTAAATAGAAACAGGAAGTACAAAAGCCCCTGAcctgcca        >  2:89045/1‑142 (MQ=255)
                                                                                                                                           ccTTACGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGCTACCTAAATTGCCAACTAAATCGAAACAGGAAGTACAAAAGTCCCTGACCTGCCTGa      <  1:667700/143‑1 (MQ=255)
                                                                                                                                               ccGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGCTACCTAAATTGCCAACTAAATCGAAACAGGAAGTACAAAAGTCCCTGACCTGCCTGatgcc  >  1:695699/1‑142 (MQ=255)
                                                                                                                                               |                                                                                                                                              
GCCTGCTCGCGCAGTTCACGCGCTAACAGAAGACCGTTCTTACCCGGCAGATTGATATCCATGATCACCAGGTTGATGTCATATTCAGAGAGGATCTGATGCATTTCCGCGCCATCTGTCGCTTCGAAAACATCATAGCCTTCCGCTTCGAAAATACTTTTCAACGTGTTGCGTGTTACCAACTCGTCTTCAACGATAAGAATGTGCGGGGTCTGCATGTTTGCTACCTAAATTGCCAACTAAATCGAAACAGGAAGTACAAAAGTCCCTGACCTGCCTGATGCA  >  NC_000913/4640089‑4640373

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