Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,643,899 G→A G428G (GGC→GGT prlC ← oligopeptidase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,643,8990GA100.0% 324.5 / NA 100G428G (GGC→GGTprlColigopeptidase A
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (52/48);  total (52/48)

CGTGACCGAACTCGTGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAA  >  NC_000913/3643763‑3644032
                                                                                                                                        |                                                                                                                                     
cGTGACCGAACTCGTGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGAcc                                                                                                                                     >  2:1144493/1‑139 (MQ=255)
 gTGACCGAACTCGTGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCt                                                                                                                                    <  1:1127052/139‑1 (MQ=255)
  tGACCGAACTCGTGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTa                                                                                                                                   >  1:871135/1‑139 (MQ=255)
   gACCGAACTCGTGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTAc                                                                                                                                  >  2:1018648/1‑139 (MQ=255)
   gACCGAACTCGTGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTAc                                                                                                                                  >  1:1379040/1‑139 (MQ=255)
              tGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTcatcca                                                                                                                       >  1:430274/1‑139 (MQ=255)
                  agagGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCAtccagcca                                                                                                                   >  1:1068710/1‑139 (MQ=255)
                     gggTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGc                                                                                                                >  2:254639/1‑139 (MQ=255)
                     gggTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGCTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGc                                                                                                                >  1:1406552/1‑139 (MQ=255)
                      ggTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGcc                                                                                                               >  1:293043/1‑139 (MQ=255)
                      ggTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGcc                                                                                                               >  2:602351/1‑139 (MQ=255)
                      ggTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGcc                                                                                                               >  2:714518/1‑139 (MQ=255)
                       gTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGccc                                                                                                              >  2:878583/1‑139 (MQ=255)
                          aTCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGc                                                                                                                                <  2:1480576/118‑1 (MQ=255)
                          aTCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGc                                                                                                                                >  1:1480576/1‑118 (MQ=255)
                          aTCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTca                                                                                                                           >  1:1170169/1‑123 (MQ=255)
                          aTCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTGTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTca                                                                                                                           <  2:1170169/123‑1 (MQ=255)
                           tCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCccgcc                                                                                                          >  1:719687/1‑139 (MQ=255)
                            cACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCccgccg                                                                                                         <  2:293043/139‑1 (MQ=255)
                                      tgagtgaACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGtttt                                                                                               >  2:1351560/1‑139 (MQ=255)
                                      tgagtgaACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGtttt                                                                                               <  1:254639/139‑1 (MQ=255)
                                       gagtgaACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTc                                                                                              >  1:84810/1‑139 (MQ=255)
                                       gagtgaACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTc                                                                                              >  2:1235474/1‑139 (MQ=255)
                                        agtgaACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCa                                                                                             <  2:1379040/139‑1 (MQ=255)
                                        agtgaACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCa                                                                                             <  1:1104187/139‑1 (MQ=255)
                                         gtgaACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCAc                                                                                            >  2:216495/1‑139 (MQ=255)
                                          tgaACAGCGCCGGTTTACCATGTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACg                                                                                           <  2:971393/139‑1 (MQ=255)
                                             aCAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGc                                                                                        >  2:180288/1‑139 (MQ=255)
                                             aCAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGc                                                                                        >  1:1356425/1‑139 (MQ=255)
                                                gcgcCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATa                                                                                     >  1:1417842/1‑139 (MQ=255)
                                                 cgcCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGtttt                                                                                               <  2:118354/128‑1 (MQ=255)
                                                 cgcCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGtttt                                                                                               >  1:118354/1‑128 (MQ=255)
                                                 cgcCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCAt                                                                                      >  2:779713/1‑137 (MQ=255)
                                                 cgcCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCAt                                                                                      <  1:779713/137‑1 (MQ=255)
                                                 cgcCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATAc                                                                                    >  1:204454/1‑139 (MQ=255)
                                                   ccGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAg                                                                                  >  1:1417740/1‑139 (MQ=255)
                                                      gTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAGAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACg                                                                                           <  2:1020337/127‑1 (MQ=255)
                                                      gTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACg                                                                                           >  1:1020337/1‑127 (MQ=255)
                                                      gTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGAt                                                                                <  2:1653625/138‑1 (MQ=255)
                                                      gTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGAt                                                                                >  1:1653625/1‑138 (MQ=255)
                                                      gTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATc                                                                               <  2:719687/139‑1 (MQ=255)
                                                      gTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATc                                                                               <  1:262694/139‑1 (MQ=255)
                                                           ccATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATa                                                                                     >  1:1002761/1‑128 (MQ=255)
                                                           ccATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTGGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATa                                                                                     <  2:1002761/128‑1 (MQ=255)
                                                              tttACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGa                                                                       <  2:1417842/139‑1 (MQ=255)
                                                                tACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCAtccatcca                                                                                                                   <  2:1562067/93‑1 (MQ=255)
                                                                tACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCAtccatcca                                                                                                                   >  1:1562067/1‑93 (MQ=255)
                                                                tACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAg                                                                     >  2:262695/1‑139 (MQ=255)
                                                                 aCCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGtttt                                                                                               >  1:1270807/1‑112 (MQ=255)
                                                                 aCCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGtttt                                                                                               <  2:1270807/112‑1 (MQ=255)
                                                                 aCCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGc                                                                    >  2:56989/1‑139 (MQ=255)
                                                                     ggCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATc                                                                               >  1:1117282/1‑124 (MQ=255)
                                                                     ggCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATc                                                                               <  2:1117282/124‑1 (MQ=255)
                                                                         gTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCgcgc                                                            >  1:75713/1‑139 (MQ=255)
                                                                          ttGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCa                                                           <  1:262695/139‑1 (MQ=255)
                                                                               gTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCg                                                      >  2:578534/1‑139 (MQ=255)
                                                                                  aCAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTa                                                                  >  1:81937/1‑124 (MQ=255)
                                                                                  aCAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTa                                                                  <  2:81937/124‑1 (MQ=255)
                                                                                   cAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTAt                                                  <  1:1144493/139‑1 (MQ=255)
                                                                                   cAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTAt                                                  <  2:819761/139‑1 (MQ=255)
                                                                                    aaGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACg                                                                                           >  1:308873/1‑97 (MQ=255)
                                                                                    aaGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACg                                                                                           <  2:308873/97‑1 (MQ=255)
                                                                                     aGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTAttt                                                <  1:56989/139‑1 (MQ=255)
                                                                                     aGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTAttt                                                <  1:1351467/139‑1 (MQ=255)
                                                                                       tCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTc                                              >  2:1324106/1‑139 (MQ=255)
                                                                                        cAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACg                                                                                           >  1:400491/1‑93 (MQ=255)
                                                                                        cAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACg                                                                                           <  2:400491/93‑1 (MQ=255)
                                                                                        cAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGtt                                                        >  1:21228/1‑128 (MQ=255)
                                                                                        cAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGtt                                                        <  2:21228/128‑1 (MQ=255)
                                                                                          aaTACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGAt                                                                                >  2:54165/1‑102 (MQ=255)
                                                                                          aaTACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGAt                                                                                <  1:54165/102‑1 (MQ=255)
                                                                                             aCGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATa                                        <  2:871135/139‑1 (MQ=255)
                                                                                             aCGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATa                                        <  1:518064/139‑1 (MQ=255)
                                                                                             aCGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATa                                        <  2:1474504/139‑1 (MQ=255)
                                                                                                 gACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATa                                                                                     >  2:495624/1‑90 (MQ=255)
                                                                                                 gACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATa                                                                                     <  1:495624/90‑1 (MQ=255)
                                                                                                 gACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGt                                    >  1:1631189/1‑139 (MQ=255)
                                                                                                   ccGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTc                                  >  2:472642/1‑139 (MQ=255)
                                                                                                    cGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCg                                 <  1:134374/139‑1 (MQ=255)
                                                                                                       tttttGCAAAGAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCgaag                              <  2:1540844/139‑1 (MQ=255)
                                                                                                                 gAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACAt                    <  2:84810/139‑1 (MQ=255)
                                                                                                                 gAACCATCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACAt                    >  1:1163056/1‑139 (MQ=255)
                                                                                                                       tCAGCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGGAGAAACGTACATCCGGAt              >  1:1481074/1‑139 (MQ=255)
                                                                                                                          gCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGcc           <  1:878583/139‑1 (MQ=255)
                                                                                                                          gCTTTACGCATCTGACCTACGCAGTCATCCATCCACGCCCCCCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGcc           <  2:1241184/139‑1 (MQ=255)
                                                                                                                            tttACGCATCTGACCTACGCAGTCATCCATCCACGCCCTGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAg         <  1:504456/139‑1 (MQ=255)
                                                                                                                            tttACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGAt                                                                                >  2:314411/1‑68 (MQ=255)
                                                                                                                            tttACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGAt                                                                                <  1:314411/68‑1 (MQ=255)
                                                                                                                            tttACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAg         <  1:1235474/139‑1 (MQ=255)
                                                                                                                              tACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCa                                                                                             <  1:380213/53‑1 (MQ=255)
                                                                                                                              tACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCa                                                                                             >  2:380213/1‑53 (MQ=255)
                                                                                                                              tACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACg                                                                                           >  1:796881/1‑55 (MQ=255)
                                                                                                                              tACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACg                                                                                           <  2:796881/55‑1 (MQ=255)
                                                                                                                              tACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCa          >  1:1569324/1‑136 (MQ=255)
                                                                                                                              tACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCa          <  2:1569324/136‑1 (MQ=255)
                                                                                                                              tACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGAc       <  1:524684/139‑1 (MQ=255)
                                                                                                                              tACGCATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGAc       >  2:229946/1‑139 (MQ=255)
                                                                                                                                  cATCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCa   <  2:740644/139‑1 (MQ=255)
                                                                                                                                   aTCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCaa  <  1:216495/139‑1 (MQ=255)
                                                                                                                                   aTCTGACCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCaa  <  2:1406552/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                     
CGTGACCGAACTCGTGGAAGAGGGTGATCACTTCGTCGTGAGTGAACAGCGCCGGTTTACCATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAA  >  NC_000913/3643763‑3644032

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: