Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,331,634 C→T A378V (GCG→GTG)  proP → proline/glycine betaine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,331,6340CT100.0% 251.5 / NA 77A378V (GCG→GTG) proPproline/glycine betaine transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (40/37);  total (40/37)

CCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTA  >  NC_000913/4331502‑4331769
                                                                                                                                    |                                                                                                                                       
cccTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGtt                                                                                                                                   >  1:1484702/1‑139 (MQ=255)
  cTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTcc                                                                                                                                 <  1:364105/139‑1 (MQ=255)
       tgtgTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTGCCTGGGAAGGTTCCGACGc                                                                                                                           >  1:80933/1‑139 (MQ=255)
        gtgtTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                           <  2:834010/139‑1 (MQ=255)
               cGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCg                                                                                                                    >  1:450266/1‑139 (MQ=255)
               cGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCg                                                                                                                    >  1:1623443/1‑139 (MQ=255)
                  tCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTa                                                                                                                 <  1:1578268/139‑1 (MQ=255)
                             aTTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTgg                                                                                                      <  1:1480900/139‑1 (MQ=255)
                                cTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcgg                                                                                                   <  2:450266/139‑1 (MQ=255)
                                  gATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                 >  1:814091/1‑139 (MQ=255)
                                      aaCAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                 <  1:85708/135‑1 (MQ=255)
                                      aaCAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                 >  2:85708/1‑135 (MQ=255)
                                          gTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCTTTACAGCGCGCTGGCGGCGGCATTTaa                                                                                         <  2:926666/139‑1 (MQ=255)
                                           tAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCa                                                                                              >  1:1819021/1‑133 (MQ=255)
                                           tAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCa                                                                                              <  2:1819021/133‑1 (MQ=255)
                                           tAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAat                                                                                        >  1:79700/1‑139 (MQ=255)
                                             aCGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                 >  2:1880639/1‑128 (MQ=255)
                                             aCGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                 <  1:1880639/128‑1 (MQ=255)
                                             aCGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAatat                                                                                      >  1:1636842/1‑139 (MQ=255)
                                                  aTCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCat                                                                                                                         >  1:850683/1‑99 (MQ=255)
                                                  aTCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCat                                                                                                                         <  2:850683/99‑1 (MQ=255)
                                                                                         aTCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCctggc                                          >  1:1407005/1‑139 (MQ=255)
                                                                                          tCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCctggct                                         >  1:379350/1‑139 (MQ=255)
                                                                                               aaCTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTcc                                                                                                                                 <  1:1724653/46‑1 (MQ=255)
                                                                                               aaCTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTcc                                                                                                                                 >  2:1724653/1‑46 (MQ=255)
                                                                                               aaCTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCg                                    <  2:1484702/139‑1 (MQ=255)
                                                                                                aCTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGa                                   <  2:1407005/139‑1 (MQ=255)
                                                                                                    cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAgcgcg                                                                                                          <  1:634267/64‑1 (MQ=255)
                                                                                                    cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAgcgcg                                                                                                          >  2:634267/1‑64 (MQ=255)
                                                                                                    cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcgg                                                                                                   <  2:1342723/71‑1 (MQ=255)
                                                                                                    cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcgg                                                                                                   >  1:1342723/1‑71 (MQ=255)
                                                                                                    cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCt                                                     >  1:1553287/1‑117 (MQ=255)
                                                                                                    cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCt                                                     <  2:1553287/117‑1 (MQ=255)
                                                                                                    cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGc                               >  2:705047/1‑139 (MQ=255)
                                                                                                    cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCACTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGc                               >  2:767449/1‑139 (MQ=255)
                                                                                                     tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                           >  2:1447009/1‑46 (MQ=255)
                                                                                                     tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                           <  1:1447009/46‑1 (MQ=255)
                                                                                                     tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                           >  1:115459/1‑46 (MQ=255)
                                                                                                     tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                           <  2:115459/46‑1 (MQ=255)
                                                                                                     tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                 >  1:1531124/1‑72 (MQ=255)
                                                                                                     tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                 <  2:1531124/72‑1 (MQ=255)
                                                                                                     tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcgg                                                                                                   <  1:1769170/70‑1 (MQ=255)
                                                                                                     tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcgg                                                                                                   >  2:1769170/1‑70 (MQ=255)
                                                                                                     tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTg                                                                        <  2:619090/97‑1 (MQ=255)
                                                                                                     tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTg                                                                        >  1:619090/1‑97 (MQ=255)
                                                                                                     tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCctgg                                           >  1:196180/1‑126 (MQ=255)
                                                                                                     tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCctgg                                           <  2:196180/126‑1 (MQ=255)
                                                                                                     tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCt                              >  2:1818356/1‑139 (MQ=255)
                                                                                                     tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCt                              >  2:414617/1‑139 (MQ=255)
                                                                                                         cGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCa                          <  2:379350/139‑1 (MQ=255)
                                                                                                         cGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCa                          <  2:79700/139‑1 (MQ=255)
                                                                                                             cGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAAt                      <  1:705047/139‑1 (MQ=255)
                                                                                                             cGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCACTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAAt                      <  1:767449/139‑1 (MQ=255)
                                                                                                                      ttCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                 <  2:591626/55‑1 (MQ=255)
                                                                                                                      ttCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                 >  1:591626/1‑55 (MQ=255)
                                                                                                                       tCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGAtgc            >  1:164803/1‑139 (MQ=255)
                                                                                                                        cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                 <  2:1518891/53‑1 (MQ=255)
                                                                                                                        cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                 >  1:1518891/1‑53 (MQ=255)
                                                                                                                        cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCa                                                                                              <  2:1826938/56‑1 (MQ=255)
                                                                                                                        cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCa                                                                                              >  1:1826938/1‑56 (MQ=255)
                                                                                                                        cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCAtt                                                                                            >  1:981520/1‑58 (MQ=255)
                                                                                                                        cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCAtt                                                                                            <  2:981520/58‑1 (MQ=255)
                                                                                                                        cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATAtt                                                                                     <  2:538546/65‑1 (MQ=255)
                                                                                                                        cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATAtt                                                                                     >  1:538546/1‑65 (MQ=255)
                                                                                                                        cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGAtgcc           <  2:1623443/139‑1 (MQ=255)
                                                                                                                        cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGAtgcc           <  2:1087652/139‑1 (MQ=255)
                                                                                                                        cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGAtgcc           >  1:1087652/1‑139 (MQ=255)
                                                                                                                              tGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                 <  2:315625/47‑1 (MQ=255)
                                                                                                                              tGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                 >  1:315625/1‑47 (MQ=255)
                                                                                                                              tGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGTCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCCGATGATGCCTGCCTa     >  2:348996/1‑139 (MQ=255)
                                                                                                                              tGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCt                              >  2:104004/1‑114 (MQ=255)
                                                                                                                              tGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCt                              <  1:104004/114‑1 (MQ=255)
                                                                                                                              tGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTa     >  1:1668122/1‑139 (MQ=255)
                                                                                                                               gCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                 <  2:701087/46‑1 (MQ=255)
                                                                                                                               gCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                 >  1:701087/1‑46 (MQ=255)
                                                                                                                               gCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTAt    >  2:556555/1‑139 (MQ=255)
                                                                                                                               gCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTAt    >  2:1865521/1‑139 (MQ=255)
                                                                                                                                 cAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTa  <  1:414617/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                       
CCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTA  >  NC_000913/4331502‑4331769

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