Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 699,852 A→G L159L (TTG→CTG)  umpH ← UMP phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913699,8520AG100.0% 213.6 / NA 60L159L (TTG→CTG) umpHUMP phosphatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/26);  total (34/26)

CGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCA  >  NC_000913/699716‑699986
                                                                                                                                        |                                                                                                                                      
cGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAgcgc                                                                                                                                     >  1:1064453/1‑140 (MQ=255)
cGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAgcgc                                                                                                                                     >  1:945551/1‑140 (MQ=255)
cGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAgcgc                                                                                                                                     >  1:798051/1‑140 (MQ=255)
     gTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCcac                                                                                                                                >  1:894243/1‑140 (MQ=255)
        aTCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAg                                                                                                                             <  2:1064453/140‑1 (MQ=255)
        aTCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAg                                                                                                                             >  2:788124/1‑140 (MQ=255)
             cGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGa                                                                                                                        >  1:1415321/1‑140 (MQ=255)
               aCAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATa                                                                                                                      >  2:707260/1‑140 (MQ=255)
                        gTTTCTTCCGAAGGCGCCTGCATTTTGTTTAATGCTGCCCGCATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCTCCCGAAAATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGCGATAAAAACcgcg                                                                                                             <  2:1654122/140‑1 (MQ=255)
                          ttcttcCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGcc                                                                                                           <  2:724912/140‑1 (MQ=255)
                          ttcttcCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGcc                                                                                                           >  1:1654122/1‑140 (MQ=255)
                                cGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTggg                                                                                                     <  2:798051/140‑1 (MQ=255)
                                              ttttGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATtggtg                                                                                       <  2:1092995/140‑1 (MQ=255)
                                                   tttAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGAt                                                                                  <  1:1141607/140‑1 (MQ=255)
                                                      aaTGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATaaa                                                                               >  2:75024/1‑140 (MQ=255)
                                                      aaTGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATaaa                                                                               >  2:454426/1‑140 (MQ=255)
                                                        tgctgcGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATtgg                                                                                         <  2:954741/128‑1 (MQ=255)
                                                        tgctgcGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATtgg                                                                                         >  1:954741/1‑128 (MQ=255)
                                                           tgcGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTg                                                                          >  1:1403822/1‑140 (MQ=255)
                                                             cgcgGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCTGATTGGTGGCGATAAAACGTGCa                                                                        >  1:1459451/1‑140 (MQ=255)
                                                                  atgatCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGtt                                                                   >  2:1318354/1‑140 (MQ=255)
                                                                      tCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGcg                                                               <  1:707260/140‑1 (MQ=255)
                                                                      tCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGcg                                                               <  2:1141585/140‑1 (MQ=255)
                                                                       ccACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGcga                                                              >  1:873110/1‑140 (MQ=255)
                                                                              ctgggctTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATa                                                       >  1:421452/1‑140 (MQ=255)
                                                                                gggctTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAgg                                                     <  1:1138390/140‑1 (MQ=255)
                                                                                gggctTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAgg                                                     >  1:1608459/1‑140 (MQ=255)
                                                                                gggctTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGAGGCGATAAAACGTGCACCGTTAGCGACGAAATAgg                                                     >  2:122045/1‑140 (MQ=255)
                                                                                 ggctTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGgc                                                    >  2:150876/1‑140 (MQ=255)
                                                                                 ggctTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGgc                                                    >  2:454457/1‑140 (MQ=255)
                                                                                           cATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGc                                          >  1:928720/1‑140 (MQ=255)
                                                                                           cATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGc                                          >  1:337577/1‑140 (MQ=255)
                                                                                            aTAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGca                                         >  2:1233908/1‑140 (MQ=255)
                                                                                            aTAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGca                                         <  1:1233908/140‑1 (MQ=255)
                                                                                               gAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatca                                      <  2:1608459/140‑1 (MQ=255)
                                                                                               gAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatca                                      <  1:1318354/140‑1 (MQ=255)
                                                                                               gAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatca                                      <  1:1032945/140‑1 (MQ=255)
                                                                                                aaCGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat                                     >  1:956594/1‑140 (MQ=255)
                                                                                                aaCGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGcatcat                                     >  1:57870/1‑140 (MQ=255)
                                                                                                aaCGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTACCGACGAAATAGGCTGCTTTATGcatcat                                     >  2:115117/1‑140 (MQ=255)
                                                                                                   ggTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTc                                  <  2:1415321/140‑1 (MQ=255)
                                                                                                   ggTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTc                                  <  2:894243/140‑1 (MQ=255)
                                                                                                     tttGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTAt                                            <  2:385594/128‑1 (MQ=255)
                                                                                                     tttGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTccc                                >  1:15439/1‑140 (MQ=255)
                                                                                                     tttGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTccc                                >  1:1183102/1‑140 (MQ=255)
                                                                                                     tttGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTccc                                <  1:75024/140‑1 (MQ=255)
                                                                                                     tttGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTccc                                <  1:454426/140‑1 (MQ=255)
                                                                                                     tttGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTcc                                 <  1:779735/139‑1 (MQ=255)
                                                                                                     tttGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTcc                                 >  2:779735/1‑139 (MQ=255)
                                                                                                     tttGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCCCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTAt                                            >  1:385594/1‑128 (MQ=255)
                                                                                                       tGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTc                                  <  2:158460/136‑1 (MQ=255)
                                                                                                       tGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTc                                  >  1:158460/1‑136 (MQ=255)
                                                                                                       tGCGCCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAg                              >  2:1193048/1‑140 (MQ=255)
                                                                                                        gcgcCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATg                                           <  1:828161/126‑1 (MQ=255)
                                                                                                        gcgcCCGGAGATTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATg                                           >  2:828161/1‑126 (MQ=255)
                                                                                                               gagaTTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTc                                  >  2:1636319/1‑128 (MQ=255)
                                                                                                               gagaTTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTc                                  <  1:1636319/128‑1 (MQ=255)
                                                                                                               gagaTTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGa                      >  2:1579100/1‑140 (MQ=255)
                                                                                                               gagaTTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGa                      >  2:830758/1‑140 (MQ=255)
                                                                                                                agaTTTTCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGaa                     <  2:1183102/140‑1 (MQ=255)
                                                                                                                     ttCTCAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGACGATACAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGt                <  2:1457645/140‑1 (MQ=255)
                                                                                                                      tctcAATCCCTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGtt               <  2:337577/140‑1 (MQ=255)
                                                                                                                               cTGCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCcaaca      <  2:15439/140‑1 (MQ=255)
                                                                                                                                 gCACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAAt    >  2:93578/1‑140 (MQ=255)
                                                                                                                                   acacaGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGAGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCa  <  1:122045/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                      
CGCAGGTTATCGCCGACAATCACCGTTTCTTCCGAATGCGCCTGCATTTTGTTTAATGCTGCGCGGATGATCCACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCA  >  NC_000913/699716‑699986

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