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,6340CT98.7% 242.9 / NA 76A378V (GCG→GTG) proPproline/glycine betaine transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (38/37);  minor base G (0/1);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTA  >  NC_000913/4331500‑4331766
                                                                                                                                      |                                                                                                                                    
tGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATg                                                                                                                                  >  1:378726/1‑139 (MQ=255)
  cccTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGtt                                                                                                                                >  2:1073646/1‑139 (MQ=255)
  cccTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGtt                                                                                                                                >  1:306221/1‑139 (MQ=255)
  cccTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGtt                                                                                                                                >  2:242628/1‑139 (MQ=255)
    cTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTcc                                                                                                                              <  2:378726/139‑1 (MQ=255)
       tttGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGa                                                                                                                           >  1:472451/1‑139 (MQ=255)
          gtgatgtctattccGGGGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGGCGGGTTACTGATGCTGGCGGTGGTCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGGGATGTTCCCGACGc                                                                                                                        <  1:1238741/127‑1 (MQ=255)
            gtTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTccc                                                                                                                             >  2:141391/1‑132 (MQ=255)
            gtTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTccc                                                                                                                             <  1:141391/132‑1 (MQ=255)
            gtTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCat                                                                                                                      <  1:1073646/139‑1 (MQ=255)
            gtTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCat                                                                                                                      >  1:491183/1‑139 (MQ=255)
            gtTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCat                                                                                                                      <  1:368501/139‑1 (MQ=255)
            gtTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCat                                                                                                                      <  2:191891/139‑1 (MQ=255)
              tGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCatat                                                                                                                    >  2:248190/1‑139 (MQ=255)
                 cGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCg                                                                                                                 >  1:1305125/1‑139 (MQ=255)
                          cGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAgcgc                                                                                                        >  2:1245556/1‑139 (MQ=255)
                                  cTGTTTAACAGTAACGTCATCGGACTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcgg                                                                                                <  2:1449133/139‑1 (MQ=255)
                                  cTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcgg                                                                                                <  1:1245556/139‑1 (MQ=255)
                                    gATTAACAGTAACGCCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCACTggcggcg                                                                                              <  2:1079221/139‑1 (MQ=255)
                                         aCAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCAtt                                                                                         <  1:248190/139‑1 (MQ=255)
                                          cAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGAATATCCGTTACAGCGCGCTggcgg                                                                                                >  1:669086/1‑131 (MQ=255)
                                          cAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGAATATCCGTTACAGCGCGCTggcgg                                                                                                <  2:669086/131‑1 (MQ=255)
                                            gTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTaa                                                                                      <  1:474611/139‑1 (MQ=255)
                                               aCGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAatat                                                                                   >  2:599144/1‑139 (MQ=255)
                                               aCGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAatat                                                                                   >  1:543484/1‑139 (MQ=255)
                                                    aTCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTgg                                                                                                   <  2:892450/118‑1 (MQ=255)
                                                    aTCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTgg                                                                                                   >  1:892450/1‑118 (MQ=255)
                                                    aTCGGCCTGATTTTTGCCAGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTCCgg                                                                              >  2:1059602/1‑139 (MQ=255)
                                                         ccTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTg                                                                         >  1:497949/1‑139 (MQ=255)
                                                                  tGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGt                                                                <  2:1305125/139‑1 (MQ=255)
                                                                   gCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTc                                                               >  1:568216/1‑139 (MQ=255)
                                                                    ccGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCt                                                              <  2:568216/139‑1 (MQ=255)
                                                                                 gCTGGCGGTGATCCTTAACTGCTTTACGGGCGATATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTg                                                 >  1:958590/1‑139 (MQ=255)
                                                                                  cTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcgg                                                                                                >  2:159467/1‑91 (MQ=255)
                                                                                  cTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcgg                                                                                                <  1:159467/91‑1 (MQ=255)
                                                                                  cTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTgg                                                >  2:1482560/1‑139 (MQ=255)
                                                                                     gCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCat                                                                                                                      <  1:986034/66‑1 (MQ=255)
                                                                                     gCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCat                                                                                                                      >  2:986034/1‑66 (MQ=255)
                                                                                      cGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGACATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTggcggc                                            >  2:854460/1‑139 (MQ=255)
                                                                                       ggTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTg                                                 <  2:384182/133‑1 (MQ=255)
                                                                                       ggTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTg                                                 >  1:384182/1‑133 (MQ=255)
                                                                                         tGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCctg                                         >  2:625042/1‑139 (MQ=255)
                                                                                           aTCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                              <  1:1273213/84‑1 (MQ=255)
                                                                                           aTCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                              >  2:1273213/1‑84 (MQ=255)
                                                                                            tCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCctggct                                      <  2:472451/139‑1 (MQ=255)
                                                                                                     gCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAg                             <  1:1491590/139‑1 (MQ=255)
                                                                                                     gCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAg                             >  1:1590840/1‑139 (MQ=255)
                                                                                                      cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTccc                                                                                                                             >  1:810326/1‑42 (MQ=255)
                                                                                                      cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTccc                                                                                                                             <  2:810326/42‑1 (MQ=255)
                                                                                                      cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                        >  2:349904/1‑47 (MQ=255)
                                                                                                      cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                        <  1:349904/47‑1 (MQ=255)
                                                                                                       tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTccc                                                                                                                             >  1:402116/1‑41 (MQ=255)
                                                                                                       tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTccc                                                                                                                             <  2:402116/41‑1 (MQ=255)
                                                                                                       tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                              >  1:1198492/1‑72 (MQ=255)
                                                                                                       tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                              <  2:1198492/72‑1 (MQ=255)
                                                                                                        ttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTc                          >  2:1343513/1‑139 (MQ=255)
                                                                                                         tACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCg                         <  1:625042/139‑1 (MQ=255)
                                                                                                         tACGGGCGCTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTTATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCg                         <  2:958590/139‑1 (MQ=255)
                                                                                                          aCGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGc                        <  1:1482560/139‑1 (MQ=255)
                                                                                                           cGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACGGCGCGCTGGCGGCGGCATTTAATATTTCGGGTCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCa                       <  1:1059602/139‑1 (MQ=255)
                                                                                                           cGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCa                       <  2:497949/139‑1 (MQ=255)
                                                                                                           cGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCa                       <  1:599144/139‑1 (MQ=255)
                                                                                                           cGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCa                       <  1:242628/139‑1 (MQ=255)
                                                                                                                 ttATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGc                            >  2:770885/1‑128 (MQ=255)
                                                                                                                 ttATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGc                            <  1:770885/128‑1 (MQ=255)
                                                                                                                   aTGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGa                                >  2:368765/1‑122 (MQ=255)
                                                                                                                   aTGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGa                                <  1:368765/122‑1 (MQ=255)
                                                                                                                         tCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGAtgc         <  2:543484/139‑1 (MQ=255)
                                                                                                                          cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcgg                                                                                                >  2:1217270/1‑51 (MQ=255)
                                                                                                                          cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcgg                                                                                                <  1:1217270/51‑1 (MQ=255)
                                                                                                                          cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGAtgcc        >  2:200025/1‑139 (MQ=255)
                                                                                                                            aCCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCa                       <  1:195274/122‑1 (MQ=255)
                                                                                                                            aCCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCa                       >  2:195274/1‑122 (MQ=255)
                                                                                                                                tGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTa  >  1:652034/1‑139 (MQ=255)
                                                                                                                                 gCCAGTGATGTTCCCGACGCATATCCGTTACAgcgcg                                                                                                       <  1:857122/37‑1 (MQ=255)
                                                                                                                                 gCCAGTGATGTTCCCGACGCATATCCGTTACAgcgcg                                                                                                       >  2:857122/1‑37 (MQ=255)
                                                                                                                                      |                                                                                                                                    
TGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTA  >  NC_000913/4331500‑4331766

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