Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,951,448 C→A 7.0% P338P (CCG→CCT)  recD ← exonuclease V (RecBCD complex), alpha chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,951,4480CA 7.0% 194.5 / 3.5 57P338P (CCG→CCT) recDexonuclease V (RecBCD complex), alpha chain
Reads supporting (aligned to +/- strand):  ref base C (24/29);  new base A (2/2);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CTGAAAAACGGTTTTCACTGCCGTTTTATCACCACGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGC  >  NC_000913/2951301‑2951594
                                                                                                                                                   |                                                                                                                                                  
ctgAAAAACGGTTTTCACTGCCGTTTTATCACCACGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCg                                                                                                                                                   <  4:157160/149‑1 (MQ=255)
   aaaaaCGGTTTTCACTGCCGTTTTATCACCACGCTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGaa                                                                                                                                                <  4:386998/149‑1 (MQ=255)
      aaCGGTTTTCACTGCCGTTTTATCACCACGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGt                                                                                                                                             >  1:413549/1‑149 (MQ=255)
      aaCGGTTTTCACTGCCGTTTTATCACCACGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGt                                                                                                                                             >  4:52930/1‑149 (MQ=255)
      aaCGGTTTTCACTGCCGTTTTATCACCACGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCAGGAACGt                                                                                                                                             <  2:326568/149‑1 (MQ=255)
       aCGGTTTTCACTGCCGTTTTATCACCACGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTg                                                                                                                                            <  2:429159/149‑1 (MQ=255)
        cGGTTTTCACTGCCGTTTTATCACCACGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGa                                                                                                                                           <  1:432720/149‑1 (MQ=255)
        cGGTTTTCACTGCCGTTTTATCACCACGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGa                                                                                                                                           <  2:392908/149‑1 (MQ=255)
                       ttttATCACCACGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAggcgg                                                                                                                            >  1:266764/1‑149 (MQ=255)
                       ttttATCACCACGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAggcgg                                                                                                                            >  4:31590/1‑149 (MQ=255)
                        tttATCACCACGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAggcggc                                                                                                                           <  2:413549/149‑1 (MQ=255)
                              accacGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCt                                                                                                                     >  2:434301/1‑149 (MQ=255)
                                 acGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGcc                                                                                                                  <  2:266764/149‑1 (MQ=255)
                                  cGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCt                                                                                                                 <  1:438624/149‑1 (MQ=255)
                                   ggTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTg                                                                                                                <  1:543862/149‑1 (MQ=255)
                                    gTTGATCGCCGCAGCTCACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGATCGTGAGTCCCCGTCAGGCAGCTAAGCCGCCTga                                                                                                               <  3:371127/149‑2 (MQ=255)
                                    gTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTgg                                                                                                               >  4:584620/1‑149 (MQ=255)
                                    gTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCAGCTAAGCCGCCTga                                                                                                               <  3:296186/149‑2 (MQ=255)
                                     ttGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGc                                                                                                              <  3:116350/149‑1 (MQ=255)
                                      tGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCa                                                                                                             <  4:67757/149‑1 (MQ=255)
                                       gATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGCCAGGCAGCTAAGCCGCCTgacgg                                                                                                            <  1:345301/149‑5 (MQ=255)
                                                aGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTcggcgg                                                                                                   >  3:411615/1‑149 (MQ=255)
                                                    aaCTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTaaa                                                                                               >  1:534273/1‑149 (MQ=255)
                                                        gACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCGGGTTCCTGCCGGAACGTGAGTACCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGccc                                                                                           >  4:230111/1‑149 (MQ=255)
                                                            aaTGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCACAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGt                                                                                       >  4:1509/1‑149 (MQ=255)
                                                             aTGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGCACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGtt                                                                                      <  4:196671/149‑1 (MQ=255)
                                                              tGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTg                                                                                     <  2:436144/149‑1 (MQ=255)
                                                               gCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTgg                                                                                    <  2:544877/149‑1 (MQ=255)
                                                                ccAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGc                                                                                   >  2:478938/1‑149 (MQ=255)
                                                                  aGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCGCGCTCGGCGGTAAAGCCCGCGTTGGCAt                                                                                 >  3:103604/1‑149 (MQ=255)
                                                                   gAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATa                                                                                >  4:69150/1‑149 (MQ=255)
                                                                    aaTCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATaa                                                                               >  2:161721/1‑149 (MQ=255)
                                                                         cTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCa                                                                          <  1:426216/149‑1 (MQ=255)
                                                                               aaaCGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATc                                                                    >  2:268514/1‑149 (MQ=255)
                                                                                       gCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTgg                                                                                    <  3:229544/125‑1 (MQ=255)
                                                                                       gCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTgg                                                                                    >  4:229544/1‑125 (MQ=255)
                                                                                         ttttttGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGcaca                                                          >  4:367342/1‑149 (MQ=255)
                                                                                          tttttGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAg                                                         >  4:502633/1‑149 (MQ=255)
                                                                                                  gcaggcagAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGcc                                                 >  4:570050/1‑149 (MQ=255)
                                                                                                      gcagAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCaa                                             <  1:526140/149‑1 (MQ=255)
                                                                                                       cagAGACTGTCGCGCAAATATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAAc                                            <  1:394608/149‑1 (MQ=255)
                                                                                                        agagACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAAcc                                           >  2:244401/1‑149 (MQ=255)
                                                                                                        agagACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAAcc                                           >  4:432052/1‑149 (MQ=255)
                                                                                                           gaCTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGTGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAg                                        <  4:361437/149‑1 (MQ=255)
                                                                                                                      aaaGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATc                             >  1:165753/1‑149 (MQ=255)
                                                                                                                       aaGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCa                            <  3:387438/149‑1 (MQ=255)
                                                                                                                        aGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCAc                           <  3:69150/149‑1 (MQ=255)
                                                                                                                               gCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGcc                    >  1:540797/1‑149 (MQ=255)
                                                                                                                                  tCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCgag                 <  3:52930/149‑1 (MQ=255)
                                                                                                                                     gtgCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGaaa              <  4:478424/149‑1 (MQ=255)
                                                                                                                                     actcacgTTCCGGCAGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGAt                        <  2:551105/133‑1 (MQ=255)
                                                                                                                                     actcacgTTCCGGCAGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGAt                        >  1:551105/7‑139 (MQ=255)
                                                                                                                                       tcacgTTCCGGCAGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGa            >  1:145761/5‑149 (MQ=255)
                                                                                                                                        ccGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGAt           <  3:355850/149‑1 (MQ=255)
                                                                                                                                          ggTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCa         >  3:298083/1‑149 (MQ=255)
                                                                                                                                            ttCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACt       <  3:31590/149‑1 (MQ=255)
                                                                                                                                                 gCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCgcgc  <  3:502633/149‑1 (MQ=255)
                                                                                                                                                  ccGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGa                               >  2:424346/1‑119 (MQ=255)
                                                                                                                                                  ccGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGa                               <  1:424346/119‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTGAAAAACGGTTTTCACTGCCGTTTTATCACCACGGTTGATCGCCGCAGCTAACTGACCAATGCCAGAATCGCTGCCGAAACGATAGCTTTTTTGCAGCAGGCAGAGACTGTCGCGCAAAGATGCCGCTTCTGTGCCGGTTCCTGCCGGAACGTGAGTCCCCGTCAGGCGGCTTAGCTGCCTGGCACGCTCGGCGGTAAAGCCCGCGTTGGCATAAGCGCAGATATCGCCCAGCACAGCCCCAGCCTCAACCGAGGCCAGTTGATCACGATCGCCGAGAAAGATCACTCGCGC  >  NC_000913/2951301‑2951594

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