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% 314.9 / ‑5.1 90C20R (TGT→CGT) pepDaminoacyl‑histidine dipeptidase (peptidase D)
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (50/39);  total (50/40)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

CGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACAGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAACACAAGAC  >  NC_000913/255530‑255794
                                                                                                                                 |                                                                                                                                       
cGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAGACCTTTCTCTTTTGCCCAACCAACAATGTCTTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGc                                                                                                                                >  2:852480/1‑139 (MQ=255)
  gTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCaa                                                                                                                              <  1:2641/139‑1 (MQ=255)
                          ttACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGgag                                                                                                      <  2:878043/139‑1 (MQ=255)
                                acctGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTT‑ATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTg                                                                                               >  1:788612/1‑139 (MQ=255)
                                acctGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTT‑ATGATAGGACGGGAGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTg                                                                                               >  1:686747/1‑139 (MQ=255)
                                     aTCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGc                                                                                                                                <  1:803995/102‑1 (MQ=255)
                                     aTCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGc                                                                                                                                >  2:803995/1‑102 (MQ=255)
                                     aTCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGAc                                                                                           >  2:377813/1‑139 (MQ=255)
                                     aTCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGAc                                                                                           >  1:730035/1‑139 (MQ=255)
                                     aTCGCGTTCGACATGGAAACCTTTCTC‑‑TTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCt                                                                                                            <  1:892328/120‑1 (MQ=255)
                                     aTCGCGTTCGACATGGAAACCTTTCTC‑‑TTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCt                                                                                                            >  2:892328/1‑120 (MQ=255)
                                       cgcgTTCGACATGGAAACCTTTCTCATTTGCCCAACCAACAATGTAGTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAg                                                                                         >  2:129582/1‑139 (MQ=255)
                                            tCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAg                                                                                    <  1:131573/139‑1 (MQ=255)
                                                   ggAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGtt                                                                             <  1:377813/139‑1 (MQ=255)
                                                            tctcTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGgag                                                                                                      <  2:929065/105‑1 (MQ=255)
                                                            tctcTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGgag                                                                                                      >  1:929065/1‑105 (MQ=255)
                                                             ctctTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCtt                                                                   >  1:903482/1‑139 (MQ=255)
                                                                    gCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATc                                                            <  2:903482/139‑1 (MQ=255)
                                                                     cccAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                            >  2:125448/1‑106 (MQ=255)
                                                                     cccAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                            <  1:125448/106‑1 (MQ=255)
                                                                          ccaacAATGTATTCAGCGAGTTGCTCTT‑ATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCa                                                     <  2:686747/139‑1 (MQ=255)
                                                                          ccaacAATGTATTCAGCGAGTTGCTCTT‑ATGATAGGACGGGTGAGGAATAGAACAGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCa                                                     <  2:788612/139‑1 (MQ=255)
                                                                            aacacTGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTACTTGTCGATCACCCGCaa                                                    <  1:852480/139‑1 (MQ=255)
                                                                            aacaaTGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCaa                                                    <  2:955342/139‑1 (MQ=255)
                                                                               aaTGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAAc                                                 >  2:28189/1‑139 (MQ=255)
                                                                                    aTTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGc                                                                                                                                >  2:1031796/1‑55 (MQ=255)
                                                                                    aTTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGc                                                                                                                                <  1:1031796/55‑1 (MQ=255)
                                                                                    aTTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAatat                                                                                                                        <  1:242519/63‑1 (MQ=255)
                                                                                    aTTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAatat                                                                                                                        >  2:242519/1‑63 (MQ=255)
                                                                                    aTTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTAt                                            >  1:671064/1‑139 (MQ=255)
                                                                                     ttCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAGTATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTAtt                                           >  2:571615/1‑139 (MQ=255)
                                                                                     ttCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTAtt                                           >  2:760819/1‑139 (MQ=255)
                                                                                     ttCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTAtt                                           >  1:589217/1‑139 (MQ=255)
                                                                                      tCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTg                                          >  2:811300/1‑139 (MQ=255)
                                                                                        aGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCa                                        >  2:438345/1‑139 (MQ=255)
                                                                                           gAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGc                                         <  1:763129/135‑1 (MQ=255)
                                                                                           gAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGAATTGc                                         >  2:763129/1‑135 (MQ=255)
                                                                                             gTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGgag                                                                                                      <  2:335518/72‑1 (MQ=255)
                                                                                             gTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGgag                                                                                                      >  1:335518/1‑72 (MQ=255)
                                                                                             gTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTca                                   >  1:622280/1‑139 (MQ=255)
                                                                                                   ctTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGt                             <  1:28189/139‑1 (MQ=255)
                                                                                                       aTGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGc                         >  1:207798/1‑139 (MQ=255)
                                                                                                              gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTgg                                                                                                        >  1:607499/1‑53 (MQ=255)
                                                                                                              gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTgg                                                                                                        <  2:607499/53‑1 (MQ=255)
                                                                                                              gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTga                                                                                              >  2:250875/1‑63 (MQ=255)
                                                                                                              gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTga                                                                                              <  1:250875/63‑1 (MQ=255)
                                                                                                              gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCt                                                                    <  1:1235827/89‑1 (MQ=255)
                                                                                                              gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCt                                                                    >  2:1235827/1‑89 (MQ=255)
                                                                                                              gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCg                                                               <  1:621169/94‑1 (MQ=255)
                                                                                                              gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCg                                                               >  2:621169/1‑94 (MQ=255)
                                                                                                              gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGc                                         <  1:359763/116‑1 (MQ=255)
                                                                                                              gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGc                                         >  2:72881/1‑116 (MQ=255)
                                                                                                              gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGc                                         >  2:359763/1‑116 (MQ=255)
                                                                                                              gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGc                                         <  1:72881/116‑1 (MQ=255)
                                                                                                              gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAt                  >  1:753868/1‑139 (MQ=255)
                                                                                                                cGGGTGAGGAATAGAACGTATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgt                >  2:717839/1‑139 (MQ=255)
                                                                                                                cGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCtgt                                                                                                          >  1:987245/1‑49 (MQ=255)
                                                                                                                cGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCtgt                                                                                                          <  2:987245/49‑1 (MQ=255)
                                                                                                                cGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCt                                                                    >  1:633107/1‑87 (MQ=255)
                                                                                                                cGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCt                                                                    <  2:633107/87‑1 (MQ=255)
                                                                                                                cGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTcac                                  >  2:231300/1‑121 (MQ=255)
                                                                                                                cGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTcac                                  <  1:231300/121‑1 (MQ=255)
                                                                                                                cGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgt                >  2:867169/1‑139 (MQ=255)
                                                                                                                 gggTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCAc                                                          <  1:54199/96‑1 (MQ=255)
                                                                                                                 gggTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCAc                                                          >  2:54199/1‑96 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCgg                                                                                                              >  2:533229/1‑42 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCgg                                                                                                              <  1:533229/42‑1 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCt                                                                                                            <  2:42045/44‑1 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCt                                                                                                            <  1:370061/44‑1 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCt                                                                                                            >  1:42045/1‑44 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCt                                                                                                            >  2:370061/1‑44 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                            >  1:307895/1‑60 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                            <  2:307895/60‑1 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGTCACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGt                             <  1:42178/123‑1 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGt                             <  2:922458/123‑1 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGt                             >  1:922458/1‑123 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGt                             >  2:42178/1‑123 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgttgt             <  2:207798/139‑1 (MQ=255)
                                                                                                                    tGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCa                                                     >  1:614200/1‑98 (MQ=255)
                                                                                                                    tGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCa                                                     <  2:614200/98‑1 (MQ=255)
                                                                                                                        gAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGc                                         >  2:677385/1‑106 (MQ=255)
                                                                                                                        gAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGc                                         <  1:677385/106‑1 (MQ=255)
                                                                                                                           tAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGc                                         >  1:526259/1‑103 (MQ=255)
                                                                                                                           tAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGc                                         <  2:526259/103‑1 (MQ=255)
                                                                                                                             gAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                            >  2:155323/1‑50 (MQ=255)
                                                                                                                             gAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                            <  1:155323/50‑1 (MQ=255)
                                                                                                                             gAACGGATGTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTAGATCACCCGCAAGACAGTCTTGCAGGTCACAGGGTTAGCCGAAAATGTTGTCAACACAAGa   >  1:1195410/1‑139 (MQ=255)
                                                                                                                              aaCGGATTTTGGCAAAAATATCCCACAGCGGCtgt                                                                                                          <  2:206275/35‑1 (MQ=255)
                                                                                                                              aaCGGATTTTGGCAAAAATATCCCACAGCGGCtgt                                                                                                          >  1:206275/1‑35 (MQ=255)
                                                                                                                              aaCGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAACACAAGAc  >  2:1186288/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: ‑