Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,037,430 A→G 5.1% L226P (CTG→CCG)  recJ ← ssDNA exonuclease, 5' ‑‑> 3'‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,037,4300AG 5.1% 237.9 / 3.2 78L226P (CTG→CCG) recJssDNA exonuclease, 5' ‑‑> 3'‑specific
Reads supporting (aligned to +/- strand):  ref base A (41/32);  new base G (2/2);  total (44/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.54e-01

CATCACGGTTTGCCACTTCAAGCAGCGCTTTAATCCCCGGACGGCACTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGTTAGG  >  NC_000913/3037284‑3037578
                                                                                                                                                  |                                                                                                                                                    
cATCACGGTTTGCCACTTCAAGCAGCGCTTTAATCCCCGGACGGCACTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCcag                                                                                                                                                     <  1:178757/148‑1 (MQ=255)
 aTCACGGTTTGCCTCTTCATTCAGCGATTTAAGCCCCGTACGGCAATTTCCGGCTCGGATGAGACTCATCCCCTGCCAGGTCAGAATGCGATTATTATCGTCCAGCGGGACGAGGTGCGCCTCTGTTTTCAGCGCGACGAGATCctgct                                                                                                                                                   <  4:203520/149‑2 (MQ=255)
 aTCACGGTTTGCCACTTCAAGCAGCGCTTTAATCCCGGGACGGCACTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGTTCAGAATGCGATTATTAGCGTCCAGCGGCACGAGGGCCGCCACTGTCCCCAGCGCGACCAGATCcagca                                                                                                                                                   <  7:63333/149‑1 (MQ=255)
 aTCACGGTTTGCCACTTCAAGCAGCGCTTTAATCCCCGGACGGCACTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCcagca                                                                                                                                                   <  5:151100/149‑1 (MQ=255)
 aTCACGGTTTGCCACTTCAAGCAGCGCTTTAATCCCCGGACGGCACTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCcagca                                                                                                                                                   <  6:11058/149‑1 (MQ=255)
  tCACGGTTTGCCACTTCAAGCAGCGCTTTAATCCCCGGACGGCACTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCcagcag                                                                                                                                                  >  3:322945/1‑149 (MQ=255)
   cACGGTTTGCCACTTCAAGCAGCGCTTTAATCCCCGGACGGCACTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGt                                                                                                                                                 >  4:205112/1‑149 (MQ=255)
      ggTTTGCCACTTCAAGCAGCGCTTTAATCCCCGGACGGCACTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCt                                                                                                                                              >  8:27433/1‑149 (MQ=255)
          tGCCACTTCAAGCAGCGCTTTAATCCCCGGACGGCACTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGc                                                                                                                                            >  1:244024/1‑147 (MQ=255)
           gCCACTTCAAGCAGCGCTTTAATCCCCGGACGGCACTTTCCGGCTCGGATGCGACTCCTCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAg                                                                                                                                         >  5:131974/1‑149 (MQ=255)
                ttCAAGCAGCGCTTTAATCCCCGGACGGCACTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCcag                                                                                                                                                     >  6:228641/1‑132 (MQ=255)
                ttCAAGCAGCGCTTTAATCCCCGGACGGCACTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCcag                                                                                                                                                     <  5:228641/132‑1 (MQ=255)
                  cAAGCAGCGCTTTAATCCCCGGACGGCACTATCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGGCCCCAGCGCGACCAGATCCAGCACTCCTGCCAGGTTAg                                                                                                                                    <  8:114696/147‑1 (MQ=255)
                          gcTTTAATCCCCGGACGGCACTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAAt                                                                                                                                <  2:233147/143‑1 (MQ=255)
                          gcTTTAATCCCCGGACGGCACTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAAt                                                                                                                                >  1:233147/1‑143 (MQ=255)
                                aTCCCCGGACGGCACTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAGTGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTAc                                                                                                                     <  6:365163/148‑1 (MQ=255)
                                 tCCCCGGACGGCACTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGc                                                                                                                   <  7:101458/149‑1 (MQ=255)
                                   cccGGACGGCACTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCtc                                                                                                                 <  1:421182/149‑1 (MQ=255)
                                   cccGGACGGCACTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGAGTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCACCAGCGCGACCAGATCCAGCGGTTCTGCCAGGTTAGGAATTGCGATGTTACGCt                                                                                                                  >  5:95911/1‑148 (MQ=255)
                                     cGGACGGCACTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCtcat                                                                                                               <  8:49084/149‑1 (MQ=255)
                                          ggCACTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAAc                                                                                                          >  6:251482/1‑149 (MQ=255)
                                            cccTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGCCGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAg                                                                                                                                         >  3:315285/3‑116 (MQ=255)
                                            cACTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAg                                                                                                                                         <  4:315285/116‑1 (MQ=255)
                                            cACTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAAc                                                                                                          >  7:246222/1‑147 (MQ=255)
                                               tttCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGcc                                                                                                     >  3:79242/1‑149 (MQ=255)
                                                   cGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGa                                                                                                 <  6:131974/149‑1 (MQ=255)
                                                   cGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGa                                                                                                 <  3:205112/149‑1 (MQ=255)
                                                     gCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAAtt                                                                                                                               >  3:48326/1‑117 (MQ=255)
                                                     gCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAAtt                                                                                                                               <  4:48326/117‑1 (MQ=255)
                                                         ggATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATcgc                                                                                             <  8:246222/147‑1 (MQ=255)
                                                          gATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCCCTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCgcgc                                                                                           >  2:31116/1‑148 (MQ=255)
                                                               gACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAg                                                                                     <  3:302010/149‑1 (MQ=255)
                                                                aCTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAgg                                                                                    >  3:27057/1‑149 (MQ=255)
                                                                     tCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTgcgc                                                                               <  7:27433/149‑1 (MQ=255)
                                                                      ccccTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCt                                                                                                                                              <  6:107294/85‑1 (MQ=255)
                                                                      ccccTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCt                                                                                                                                              >  5:119609/1‑85 (MQ=255)
                                                                      ccccTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCt                                                                                                                                              >  5:107294/1‑85 (MQ=255)
                                                                      ccccTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCt                                                                                                                                              <  6:119609/85‑1 (MQ=255)
                                                                         cTGCCAGGTCAGAATCCGATTATTAGCGTCCAGCGGCAAGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCACCAGTTCTGCCAGGTTAGCAATAGCGATGTTACCCTCATCAAACAATCCCTGATCCCGCAAAAAGGTGCGCacag                                                                           >  1:150607/1‑146 (MQ=255)
                                                                            ccAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCa                                                                        <  4:94482/149‑1 (MQ=255)
                                                                             cAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAg                                                                       <  2:244024/149‑1 (MQ=255)
                                                                                gTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGGTCTGCCAGGTTAGGAATTGCGATGGTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCAt                                                                    >  5:162876/1‑149 (MQ=255)
                                                                                 tCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATc                                                                   <  3:314255/149‑1 (MQ=255)
                                                                                   aGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACGGTCCCCAGCGCGACCAGATCCCGCAGTTCTGCCAGGTTAGGAATTGCGAGGTTACGCTCATCAAACCAGCCCTGATCGCGCAAACAGGTGCGCAGCGCCAGCAt                                                                    >  2:420083/1‑146 (MQ=255)
                                                                                    gAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGa                                                                <  4:79242/149‑1 (MQ=255)
                                                                                     aaTGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGa                                                                >  6:353981/1‑148 (MQ=255)
                                                                                         cGATTATTAGCGTCCAGCTGCACGACGTCCGCCACTGTCGCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGGTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTTATCTCGCAAAAAGGTGCGCAGCGCCAGCATCAGATaaaa                                                           <  2:150607/149‑1 (MQ=255)
                                                                                         cGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATaaaa                                                           <  5:349722/149‑1 (MQ=255)
                                                                                               ttaGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCca                                                      <  4:27057/148‑1 (MQ=255)
                                                                                                  gCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACAcc                                                  >  2:150654/1‑149 (MQ=255)
                                                                                                    gTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCa                                                <  5:251482/149‑1 (MQ=255)
                                                                                                           ggCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGc                                         <  4:322945/149‑1 (MQ=255)
                                                                                                              acgacgTCCGCCACTGTCCCCAGCGCGACCAGATCCGGCAGTTCTGCCAGGTTAGGAATTGCGATGTTAc                                                                                                                     >  4:204461/1‑70 (MQ=255)
                                                                                                              acgacgTCCGCCACTGTCCCCAGCGCGACCAGATCCGGCAGTTCTGCCAGGTTAGGAATTGCGATGTTAc                                                                                                                     <  3:204461/70‑1 (MQ=255)
                                                                                                              acgacgTCCGCCACTGTCCCCAGCGCGACCAGATCCGGCAGTTCTGCCAGGTTAGGAATTGCGATGTTAc                                                                                                                     >  8:18462/1‑70 (MQ=255)
                                                                                                              acgacgTCCGCCACTGTCCCCAGCGCGACCAGATCCGGCAGTTCTGCCAGGTTAGGAATTGCGATGTTAc                                                                                                                     <  7:18462/70‑1 (MQ=255)
                                                                                                                gacgTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGt                                     >  4:345834/1‑148 (MQ=255)
                                                                                                                gacgTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGAACGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGt                                     <  3:345834/148‑1 (MQ=255)
                                                                                                                   gTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGAtt                                 >  6:157569/1‑149 (MQ=255)
                                                                                                                   gTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGAtt                                 >  6:151341/1‑149 (MQ=255)
                                                                                                                    tCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATc                                                                                               >  3:123946/1‑86 (MQ=255)
                                                                                                                    tCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATc                                                                                               <  8:103596/86‑1 (MQ=255)
                                                                                                                    tCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATc                                                                                               >  7:103596/1‑86 (MQ=255)
                                                                                                                    tCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATc                                                                                               <  4:123946/86‑1 (MQ=255)
                                                                                                                     ccGCCACTGTCCCCAGCGCGCCCAGCTCCAGCAGTTCTGCCAGGTTAGGAATTGCGCTGTTACGCTCATCAAACCAGCCCTGCTCGCGCAAAAAGGTGCGCAGCGCCAGCGTCAGATAACACGCCACGCCCCCGCCTGCCGGGGAGTTc                               >  7:376743/1‑149 (MQ=255)
                                                                                                                            tGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACggg                         >  2:30800/1‑148 (MQ=255)
                                                                                                                               ccccAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAAt                     <  2:369941/149‑1 (MQ=255)
                                                                                                                                 ccAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTa                   >  5:284894/1‑149 (MQ=255)
                                                                                                                                  cAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTa                   >  8:225302/1‑148 (MQ=255)
                                                                                                                                        gACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATaaaa                                                           <  8:187792/102‑1 (MQ=255)
                                                                                                                                        gACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATaaaa                                                           >  7:187792/1‑102 (MQ=255)
                                                                                                                                        gACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGa                        <  4:14006/137‑1 (MQ=255)
                                                                                                                                        gACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGa                        >  3:14006/1‑137 (MQ=255)
                                                                                                                                        gACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTcgc            >  2:211048/1‑149 (MQ=255)
                                                                                                                                         aCCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTcgcg           >  8:277417/1‑149 (MQ=255)
                                                                                                                                          ccAGATCCAGCAGTTCTGCCAGGTTAGGAATTGAGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGCCGGGAAATTACAGTcgcg           >  2:73574/1‑148 (MQ=255)
                                                                                                                                            aGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATc                                                                   <  5:222247/90‑1 (MQ=255)
                                                                                                                                            aGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATc                                                                   >  6:222247/1‑90 (MQ=255)
                                                                                                                                             gATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAg       >  2:336867/1‑149 (MQ=255)
                                                                                                                                              aTCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAgg                                                                                    <  1:296095/71‑1 (MQ=255)
                                                                                                                                              aTCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAgg                                                                                    >  2:296095/1‑71 (MQ=255)
                                                                                                                                              aTCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACTGGAAATTACAGTCGCGCAAGt      >  3:48388/1‑149 (MQ=255)
                                                                                                                                               tCCAGCAGTTCTGCCAGGTTAGGAATTGGGATGTTACGCTCATCACACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCCCGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGt      >  2:150778/1‑148 (MQ=255)
                                                                                                                                               tCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCCGCCCTGCTCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGt      >  1:219254/1‑148 (MQ=255)
                                                                                                                                               tCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGt      >  2:173118/1‑148 (MQ=255)
                                                                                                                                               tCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGt      >  3:126623/1‑147 (MQ=255)
                                                                                                                                                 cagcagTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGTTa    <  6:284894/148‑1 (MQ=255)
                                                                                                                                                  agcagTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAAATACAGTCGCGCAAGTTAgg  >  3:291691/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CATCACGGTTTGCCACTTCAAGCAGCGCTTTAATCCCCGGACGGCACTTTCCGGCTCGGATGCGACTCATCCCCTGCCAGGTCAGAATGCGATTATTAGCGTCCAGCGGCACGACGTCCGCCACTGTCCCCAGCGCGACCAGATCCAGCAGTTCTGCCAGGTTAGGAATTGCGATGTTACGCTCATCAAACCAGCCCTGATCGCGCAAAAAGGTGCGCAGCGCCAGCATCAGATAAAACGCCACACCCACGCCTGCCAGTGATTTCGACGGGAAATTACAGTCGCGCAAGTTAGG  >  NC_000913/3037284‑3037578

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