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.9% 278.8 / ‑5.8 87A378V (GCG→GTG) proPproline/glycine betaine transporter
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (44/42);  total (44/43)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTA  >  NC_000913/4331499‑4331769
                                                                                                                                       |                                                                                                                                       
ttGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGAt                                                                                                                                      >  1:213014/1‑139 (MQ=255)
ttGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGAt                                                                                                                                      >  1:850227/1‑139 (MQ=255)
           gtgtTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                           <  2:213014/139‑1 (MQ=255)
           gtgtTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                           <  1:1376054/139‑1 (MQ=255)
              ttGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCata                                                                                                                        >  1:1610619/1‑139 (MQ=255)
                           cGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                           >  2:381508/1‑123 (MQ=255)
                           cGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                           <  1:381508/123‑1 (MQ=255)
                                   cTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcgg                                                                                                   <  2:344868/139‑1 (MQ=255)
                                     gATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                 <  2:850227/139‑1 (MQ=255)
                                     gATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                 >  1:790354/1‑139 (MQ=255)
                                     gATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcgg                                                                                                   <  1:124154/137‑1 (MQ=255)
                                     gATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcgg                                                                                                   >  2:124154/1‑137 (MQ=255)
                                         aaCAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCAt                                                                                             >  2:834207/1‑139 (MQ=255)
                                                acctCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGc                                                                                                                           <  1:1051719/99‑1 (MQ=255)
                                                aCGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACGTTGCCAGTGATGTTCCCGACGc                                                                                                                           >  2:1051719/1‑102 (MQ=255)
                                                aCGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAg                                                                                                              >  1:1079458/1‑115 (MQ=255)
                                                aCGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAg                                                                                                              <  2:1079458/115‑1 (MQ=255)
                                                aCGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCAtt                                                                                            >  1:1475444/1‑133 (MQ=255)
                                                aCGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCAtt                                                                                            <  2:1475444/133‑1 (MQ=255)
                                                aCGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAatat                                                                                      >  2:171628/1‑139 (MQ=255)
                                                aCGTCATCGGCCTGATTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                 <  2:21692/127‑1 (MQ=255)
                                                aCGTCATCGGCCTGATTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcg                                                                                                 >  1:21692/1‑127 (MQ=255)
                                                  gTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATAttt                                                                                    <  1:835845/139‑1 (MQ=255)
                                                     aTCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                           <  2:1593558/97‑1 (MQ=255)
                                                     aTCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                           >  1:1593558/1‑97 (MQ=255)
                                                        ggCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGc                                                                              >  2:449320/1‑139 (MQ=255)
                                                        ggCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGc                                                                              >  1:958982/1‑139 (MQ=255)
                                                                        ggTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                           >  2:48880/1‑78 (MQ=255)
                                                                        ggTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                           <  1:48880/78‑1 (MQ=255)
                                                                            aCTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCa                                                          >  2:1661374/1‑139 (MQ=255)
                                                                                       cGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTggcggc                                               <  1:171628/139‑1 (MQ=255)
                                                                                          tGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCctg                                            >  1:1556479/1‑139 (MQ=255)
                                                                                            aTCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCat                                                                                                                         >  2:467962/1‑60 (MQ=255)
                                                                                            aTCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCat                                                                                                                         <  1:467962/60‑1 (MQ=255)
                                                                                                ttAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcgg                                                                                                   >  2:118881/1‑78 (MQ=255)
                                                                                                ttAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcgg                                                                                                   <  1:118881/78‑1 (MQ=255)
                                                                                                      gCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCa                                                                                              >  1:1624704/1‑77 (MQ=255)
                                                                                                      gCTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCa                                                                                              <  2:1624704/77‑1 (MQ=255)
                                                                                                       cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCatat                                                                                                                       >  1:371137/1‑51 (MQ=255)
                                                                                                       cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCatat                                                                                                                       <  2:371137/51‑1 (MQ=255)
                                                                                                       cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCata                                                                                                                        >  1:589124/1‑50 (MQ=255)
                                                                                                       cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCata                                                                                                                        <  2:589124/50‑1 (MQ=255)
                                                                                                       cTTTACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGc                               >  2:352013/1‑139 (MQ=255)
                                                                                                        tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                           <  2:1667792/46‑1 (MQ=255)
                                                                                                        tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                           >  1:1667792/1‑46 (MQ=255)
                                                                                                        tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCat                                                                                                                         <  2:214833/48‑1 (MQ=255)
                                                                                                        tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCat                                                                                                                         >  1:214833/1‑48 (MQ=255)
                                                                                                        tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTaa                                                                                         >  1:1507542/1‑80 (MQ=255)
                                                                                                        tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTaa                                                                                         <  2:1507542/80‑1 (MQ=255)
                                                                                                        tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCctgg                                           <  2:120475/126‑1 (MQ=255)
                                                                                                        tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCctgg                                           >  1:120475/1‑126 (MQ=255)
                                                                                                        tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCt                              >  2:100312/1‑139 (MQ=255)
                                                                                                        tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCt                              >  2:1426246/1‑139 (MQ=255)
                                                                                                        tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCt                              >  2:783405/1‑139 (MQ=255)
                                                                                                        tttACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCt                              >  2:1571591/1‑139 (MQ=255)
                                                                                                          tACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCg                            <  1:352013/139‑1 (MQ=255)
                                                                                                          tACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCg                            <  1:1571591/139‑1 (MQ=255)
                                                                                                          tACGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCg                            <  1:100312/139‑1 (MQ=255)
                                                                                                            cGGGCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCa                          <  1:834207/139‑1 (MQ=255)
                                                                                                              ggCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGa                        <  2:958982/139‑1 (MQ=255)
                                                                                                               gCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCa                                                                                              <  2:1638294/68‑1 (MQ=255)
                                                                                                               gCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCa                                                                                              >  1:1638294/1‑68 (MQ=255)
                                                                                                               gCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGa                        <  2:422950/138‑1 (MQ=255)
                                                                                                               gCGTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGa                        >  1:422950/1‑138 (MQ=255)
                                                                                                                 gTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                           >  1:1172228/1‑37 (MQ=255)
                                                                                                                 gTTATGGCTTCTACCTTGCCAGTGATGTTCCCGACGc                                                                                                                           <  2:1172228/37‑1 (MQ=255)
                                                                                                                    aTGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTACTATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCt                                                     >  2:1013153/1‑104 (MQ=255)
                                                                                                                    aTGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCt                                                     <  1:1013153/104‑1 (MQ=255)
                                                                                                                     tGGCTTCTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCtgat                 <  1:1661374/139‑1 (MQ=255)
                                                                                                                           cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcggc                                                                                               <  2:118027/55‑1 (MQ=255)
                                                                                                                           cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTggcggcggc                                                                                               >  1:118027/1‑55 (MQ=255)
                                                                                                                           cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCt                                                     <  2:1438887/97‑1 (MQ=255)
                                                                                                                           cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCt                                                     >  1:1438887/1‑97 (MQ=255)
                                                                                                                           cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGAtgcc           <  1:783405/139‑1 (MQ=255)
                                                                                                                           cTACCTTGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGAtgcc           <  2:1636974/139‑1 (MQ=255)
                                                                                                                            tACCTTGCCAGTGATGTTCCCGACGCATATCCGTTa                                                                                                                 >  1:1577533/1‑36 (MQ=255)
                                                                                                                            tACCTTGCCAGTGATGTTCCCGACGCATATCCGTTa                                                                                                                 <  2:1577533/36‑1 (MQ=255)
                                                                                                                                 tGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCaa                                                         >  1:387236/1‑87 (MQ=255)
                                                                                                                                 tGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCaa                                                         <  2:387236/87‑1 (MQ=255)
                                                                                                                                 tGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTg                                                    >  1:285502/1‑92 (MQ=255)
                                                                                                                                 tGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTg                                                    <  2:285502/92‑1 (MQ=255)
                                                                                                                                 tGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTa     >  2:437783/1‑139 (MQ=255)
                                                                                                                                 tGCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTa     >  1:1336004/1‑139 (MQ=255)
                                                                                                                                  gCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGAtgcc           >  1:1715079/1‑132 (MQ=255)
                                                                                                                                  gCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGAtgcc           <  2:1715079/132‑1 (MQ=255)
                                                                                                                                  gCCAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTAt    <  1:449320/139‑1 (MQ=255)
                                                                                                                                    cAGTGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTa  <  1:437783/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
TTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTA  >  NC_000913/4331499‑4331769

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: