Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,022,508 A→G V459A (GTG→GCG)  yccS ← putative transporter, FUSC superfamily inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,022,5080AG97.1% 251.8 / ‑2.1 71V459A (GTG→GCG) yccSputative transporter, FUSC superfamily inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (29/39);  total (31/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.23e-01

TTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGTGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGA  >  NC_000913/1022372‑1022642
                                                                                                                                        |                                                                                                                                      
ttCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGcgca                                                                                                                                      >  1:821746/1‑139 (MQ=255)
     ccAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACca                                                                                                                                 <  1:360026/139‑1 (MQ=255)
        gTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACcacag                                                                                                                              <  2:478222/139‑1 (MQ=255)
                     caGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAgg                                                                                                                        <  1:687009/132‑1 (MQ=255)
                     caGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAgg                                                                                                                        >  2:687009/1‑132 (MQ=255)
                       gTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAAt                                                                                                               >  1:1231974/1‑139 (MQ=255)
                           cACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGcc                                                                                                                  >  1:195682/1‑132 (MQ=255)
                           cACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGcc                                                                                                                  <  2:195682/132‑1 (MQ=255)
                           cACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCg                                                                                                           <  1:973148/139‑1 (MQ=255)
                           cACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCg                                                                                                           <  1:495644/139‑1 (MQ=255)
                           cACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCg                                                                                                           <  1:1101939/139‑1 (MQ=255)
                             ccAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGAt                                                                                                         <  1:300489/139‑1 (MQ=255)
                                aaaGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACcacagca                                                                                                                            <  2:425603/117‑1 (MQ=255)
                                aaaGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACcacagca                                                                                                                            >  1:425603/1‑117 (MQ=255)
                                aaaGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGAtacc                                                                                                      >  1:909862/1‑139 (MQ=255)
                                  aGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGAtaccta                                                                                                    >  1:603613/1‑139 (MQ=255)
                                   gtgtGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGAtacctac                                                                                                   <  1:737309/139‑1 (MQ=255)
                                   gtgtGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGAtacctac                                                                                                   <  2:512264/139‑1 (MQ=255)
                                        atgaACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCg                                                                                              <  2:577660/139‑1 (MQ=255)
                                        atgaACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCg                                                                                              >  1:1032594/1‑139 (MQ=255)
                                         tgaACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGt                                                                                             <  1:435147/139‑1 (MQ=255)
                                              aTCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCaa                                                                                        <  1:1216924/139‑1 (MQ=255)
                                              aTCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCaa                                                                                        <  1:1144368/139‑1 (MQ=255)
                                              aTCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGAACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGACAATGGCGATACCTACCAGCGTACCaa                                                                                        <  2:796994/139‑1 (MQ=255)
                                                cGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAg                                                                                                <  2:1071883/129‑1 (MQ=255)
                                                cGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAg                                                                                                >  1:1071883/1‑129 (MQ=255)
                                                  ttGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGa                                                                                                                       >  1:998546/1‑104 (MQ=255)
                                                  ttGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGa                                                                                                                       <  2:998546/104‑1 (MQ=255)
                                                      aTGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCtt                                                                                >  1:1315755/1‑139 (MQ=255)
                                                          gCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAAc                                                                                                                                   >  1:135100/1‑84 (MQ=255)
                                                          gCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAAc                                                                                                                                   <  2:135100/84‑1 (MQ=255)
                                                           cGTATTGCACGATACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCAAGTGATTGCCCAAACCACAGCACAGGAATGACAATGGCGATACCTACCAGCGTACCAATAATCCTTATCtt                                                                           >  1:965713/1‑139 (MQ=255)
                                                                 gCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCg                                                                     <  1:289452/139‑1 (MQ=255)
                                                                     gTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGAtacctacc                                                                                                  >  1:1254889/1‑106 (MQ=255)
                                                                     gTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGAtacctacc                                                                                                  <  2:1254889/106‑1 (MQ=255)
                                                                       tACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGcc                                                                                                                  <  2:2283/88‑1 (MQ=255)
                                                                       tACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGcc                                                                                                                  >  1:2283/1‑88 (MQ=255)
                                                                       tACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGgcg                                                               >  2:321976/1‑139 (MQ=255)
                                                                         cGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGcgcg                                                             <  1:249095/139‑1 (MQ=255)
                                                                         cGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGcgcg                                                             <  1:146386/139‑1 (MQ=255)
                                                                           gAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGTg                                                           <  2:540826/139‑1 (MQ=255)
                                                                                cAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGAGGCGtt                                                      >  1:411345/1‑139 (MQ=255)
                                                                                 aaaaaaGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGcc                                                                                                                  >  2:565803/1‑78 (MQ=255)
                                                                                 aaaaaaGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGcc                                                                                                                  <  1:565803/78‑1 (MQ=255)
                                                                                 aaaaaaGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATgg                                                                                                             <  1:1160927/83‑1 (MQ=255)
                                                                                 aaaaaaGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATgg                                                                                                             >  2:1160927/1‑83 (MQ=255)
                                                                                 aaaaaaGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGTGGCGTta                                                     >  1:1062319/1‑139 (MQ=255)
                                                                                 aaaaaaGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGTGGCGTta                                                     >  1:933270/1‑139 (MQ=255)
                                                                                        agCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGTGGCGTTATAGTTtg                                              >  1:1105960/1‑139 (MQ=255)
                                                                                          cACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACca                                                                                                                                 >  2:614093/1‑54 (MQ=255)
                                                                                          cACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACca                                                                                                                                 <  1:614093/54‑1 (MQ=255)
                                                                                             gCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGTGGCGTTATAGTTtggctgg                                         <  1:999864/139‑1 (MQ=255)
                                                                                               cGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGTg                                                           >  1:284511/1‑119 (MQ=255)
                                                                                               cGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGTg                                                           <  2:284511/119‑1 (MQ=255)
                                                                                               cGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGTGGCGTTATAGTTTGGCTGGCa                                       >  2:4073/1‑139 (MQ=255)
                                                                                                           gcagcaCCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGTGGCGTTATAGTTTGGCTGGCAGACAAACAAACt                           >  1:50024/1‑139 (MQ=255)
                                                                                                            cagcaCCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGTGGCGTTATAGTTTGGCTGGCAGACAAACAAACtt                          <  2:1105960/139‑1 (MQ=255)
                                                                                                               caCCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCg                                                                                              >  1:622968/1‑68 (MQ=255)
                                                                                                               caCCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCg                                                                                              <  2:622968/68‑1 (MQ=255)
                                                                                                                 ccAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGTGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAg                     >  1:1232448/1‑139 (MQ=255)
                                                                                                                   aGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGTGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCa                   <  2:1231974/139‑1 (MQ=255)
                                                                                                                    gCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCg                                                                                                           <  2:436565/50‑1 (MQ=255)
                                                                                                                    gCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCg                                                                                                           >  1:436565/1‑50 (MQ=255)
                                                                                                                     cTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGTGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAg                     >  1:264864/1‑135 (MQ=255)
                                                                                                                     cTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGTGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAg                     <  2:264864/135‑1 (MQ=255)
                                                                                                                      tGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGTGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGa                <  1:321976/139‑1 (MQ=255)
                                                                                                                       gCCCTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGTGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGAt               >  1:112091/1‑139 (MQ=255)
                                                                                                                         ccTTCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGTGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATcc             <  2:249402/139‑1 (MQ=255)
                                                                                                                           ttCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGTGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCaa           <  2:821746/139‑1 (MQ=255)
                                                                                                                            tCCAGTGATGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGTGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAAt          <  2:603613/139‑1 (MQ=255)
                                                                                                                                    tGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGTGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGa  <  2:555018/139‑1 (MQ=255)
                                                                                                                                    tGGCGCAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGTGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGa  >  1:206679/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                      
TTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGTGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGA  >  NC_000913/1022372‑1022642

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