Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 255,659 A→G C20R (TGT→CGT)  pepD ← aminoacyl‑histidine dipeptidase (peptidase D)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913255,6590AG98.9% 324.6 / ‑4.8 90C20R (TGT→CGT) pepDaminoacyl‑histidine dipeptidase (peptidase D)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (48/41);  minor base T (0/1);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACAGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAACACAAGAC  >  NC_000913/255530‑255794
                                                                                                                                 |                                                                                                                                       
cGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGc                                                                                                                                <  2:1328764/139‑1 (MQ=255)
      cAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAat                                                                                                                          >  1:1401206/1‑139 (MQ=255)
             aCGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCcac                                                                                                                   <  2:175338/139‑1 (MQ=255)
                 aTCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCg                                                                                                               <  2:1401206/139‑1 (MQ=255)
                          ttACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGgag                                                                                                      <  2:640283/139‑1 (MQ=255)
                           tacctaccTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAatat                                                                                                                        >  2:1304330/1‑120 (MQ=255)
                           tacctaccTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAatat                                                                                                                        <  1:1304330/120‑1 (MQ=255)
                                    gATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                            <  2:894605/139‑1 (MQ=255)
                                     aTCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCaaaa                                                                                                                            >  2:931274/1‑106 (MQ=255)
                                     aTCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCaaaa                                                                                                                            <  1:931274/106‑1 (MQ=255)
                                     aTCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAata                                                                                                                         >  2:315824/1‑109 (MQ=255)
                                     aTCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAata                                                                                                                         <  1:315824/109‑1 (MQ=255)
                                     aTCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAg                                                                                                                 <  1:1076112/117‑1 (MQ=255)
                                     aTCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAg                                                                                                                 >  2:1076112/1‑117 (MQ=255)
                                     aTCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCt                                                                                                            >  1:876840/1‑122 (MQ=255)
                                     aTCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCt                                                                                                            <  2:876840/122‑1 (MQ=255)
                                     aTCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGAc                                                                                           >  1:235767/1‑139 (MQ=255)
                                      tCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCaa                                                                                                                              >  1:1411287/1‑103 (MQ=255)
                                      tCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCaa                                                                                                                              <  2:1411287/103‑1 (MQ=255)
                                      tCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACa                                                                                          >  1:1184856/1‑139 (MQ=255)
                                         cgTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGtt                                                                                       <  1:1315118/139‑1 (MQ=255)
                                           ttCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                            <  2:334040/132‑1 (MQ=255)
                                           ttCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                            >  1:334040/1‑132 (MQ=255)
                                           ttCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCa                                                                                     >  2:1688689/1‑139 (MQ=255)
                                                          tttCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACTGATTTTGGCAAAAATAGCCCACAg                                                                                                                 <  1:1033630/96‑1 (MQ=255)
                                                          tttCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAGATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGtctc                                                                      >  2:1112652/1‑139 (MQ=255)
                                                          tttCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAg                                                                                                                 >  2:1033630/1‑96 (MQ=255)
                                                          tttCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGtctc                                                                      >  2:743464/1‑139 (MQ=255)
                                                             ctctTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCtt                                                                   >  1:476253/1‑139 (MQ=255)
                                                                   tGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCaaaa                                                                                                                            >  2:1357308/1‑76 (MQ=255)
                                                                   tGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCaaaa                                                                                                                            <  1:1357308/76‑1 (MQ=255)
                                                                   tGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGt                                                                         >  1:1295410/1‑127 (MQ=255)
                                                                   tGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGt                                                                         <  2:1295410/127‑1 (MQ=255)
                                                                   tGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGAt                                                             <  1:1688689/139‑1 (MQ=255)
                                                                          ccaacAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCtgt                                                                                                          <  1:965415/87‑1 (MQ=255)
                                                                          ccaacAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCtgt                                                                                                          >  2:965415/1‑87 (MQ=255)
                                                                               aaTGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGgag                                                                                                      >  1:41778/1‑86 (MQ=255)
                                                                               aaTGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGgag                                                                                                      <  2:41778/86‑1 (MQ=255)
                                                                                 tGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGCCGATCACCCGCAAAACAg                                               <  1:1112652/139‑1 (MQ=255)
                                                                                     ttCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGAtt                                                                                                                                     >  1:1101135/1‑49 (MQ=255)
                                                                                     ttCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGAtt                                                                                                                                     <  2:1101135/49‑1 (MQ=255)
                                                                                     ttCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGAc                                                                                           >  1:403916/1‑91 (MQ=255)
                                                                                     ttCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGAc                                                                                           <  2:403916/91‑1 (MQ=255)
                                                                                     ttCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTAtt                                           >  2:207876/1‑139 (MQ=255)
                                                                                      tCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGa                                                                                   >  2:409235/1‑98 (MQ=255)
                                                                                      tCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGa                                                                                   <  1:409235/98‑1 (MQ=255)
                                                                                      tCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTg                                          <  1:743464/139‑1 (MQ=255)
                                                                                        aGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAg                                                                                                                 >  2:1582990/1‑66 (MQ=255)
                                                                                        aGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAg                                                                                                                 <  1:1582990/66‑1 (MQ=255)
                                                                                        aGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTGTTGCa                                        >  2:476979/1‑139 (MQ=255)
                                                                                        aGCGAGTTGCTCTTCATGACAGGACGGGTGAGGGATAGAACGGATTTTGGCAAAAATATCCCTCAGCGGATGTGGAGATAATTGAGATAGTTCATACACGTTAAGTCTCCTTGGCGATCACCCGCAAAACAGTATTGCa                                        <  2:667668/139‑1 (MQ=255)
                                                                                         gCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAACAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAg                                       >  1:433690/1‑139 (MQ=255)
                                                                                         gCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAg                                       <  2:476253/139‑1 (MQ=255)
                                                                                            aGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTc                                    <  1:950754/139‑1 (MQ=255)
                                                                                             gTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAg                                                                                                                 >  2:464216/1‑61 (MQ=255)
                                                                                             gTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAg                                                                                                                 <  1:464216/61‑1 (MQ=255)
                                                                                             gTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGa                                                                                   >  1:1618071/1‑91 (MQ=255)
                                                                                             gTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGa                                                                                   <  2:1618071/91‑1 (MQ=255)
                                                                                             gTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTca                                   >  1:1239335/1‑139 (MQ=255)
                                                                                             gTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACAAGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTca                                   >  1:872218/1‑139 (MQ=255)
                                                                                              tttctctTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTCTTGCAGGTcac                                  <  2:433690/136‑1 (MQ=255)
                                                                                                      cATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTcac                                  <  2:893174/131‑1 (MQ=255)
                                                                                                      cATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTcac                                  >  1:893174/1‑131 (MQ=255)
                                                                                                           tAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGaa                     >  2:1284677/1‑139 (MQ=255)
                                                                                                              gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGTTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAt                  >  1:564225/1‑139 (MQ=255)
                                                                                                              gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAt                  >  2:566557/1‑139 (MQ=255)
                                                                                                              gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAt                  >  1:418480/1‑139 (MQ=255)
                                                                                                              gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGACGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAt                  >  2:582461/1‑139 (MQ=255)
                                                                                                                cGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGa                                                                                   <  2:189832/72‑1 (MQ=255)
                                                                                                                cGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGa                                                                                   >  1:189832/1‑72 (MQ=255)
                                                                                                                cGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTaa                                                                           >  2:207412/1‑80 (MQ=255)
                                                                                                                cGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTaa                                                                           <  1:207412/80‑1 (MQ=255)
                                                                                                                cGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgt                >  1:14973/1‑139 (MQ=255)
                                                                                                                cGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgt                >  2:1323861/1‑139 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGAc                                                                                           <  1:881870/61‑1 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGAc                                                                                           >  2:881870/1‑61 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGa                                                                                   <  1:1211329/69‑1 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGa                                                                                   >  2:1211329/1‑69 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTcc                                                                     >  2:128972/1‑83 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTcc                                                                     <  1:128972/83‑1 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGc                                         >  2:614302/1‑111 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGc                                         <  1:614302/111‑1 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgttgt             <  1:1323861/139‑1 (MQ=255)
                                                                                                                    tGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                            >  2:559304/1‑59 (MQ=255)
                                                                                                                    tGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                            <  1:559304/59‑1 (MQ=255)
                                                                                                                      aGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCgg                                                                                                              <  2:129520/39‑1 (MQ=255)
                                                                                                                      aGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCgg                                                                                                              >  1:129520/1‑39 (MQ=255)
                                                                                                                      aGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGACGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCaa          <  1:582461/139‑1 (MQ=255)
                                                                                                                             gAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAACACAAGa   >  2:810973/1‑139 (MQ=255)
                                                                                                                              aaCGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAACACAAGAc  >  2:1249785/1‑139 (MQ=255)
                                                                                                                                 |                                                                                                                                       
CGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACAGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAACACAAGAC  >  NC_000913/255530‑255794

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: