Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 209,806 T→C 5.2% L43P (CTG→CCG)  ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913209,8060TC 5.2% 237.1 / 2.2 77L43P (CTG→CCG) ldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base T (42/31);  new base C (2/2);  total (44/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

AGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGC  >  NC_000913/209658‑209954
                                                                                                                                                    |                                                                                                                                                    
aggaTTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCt                                                                                                                                                      >  8:224909/1‑149 (MQ=255)
aggaTTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCt                                                                                                                                                      >  6:182680/1‑149 (MQ=255)
  gaTTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGATGATTTGCTg                                                                                                                                                     >  2:175300/1‑148 (MQ=255)
    ttttCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGaaa                                                                                                                                                  <  6:251708/149‑1 (MQ=255)
    ttttCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCCGaaa                                                                                                                                                  >  4:315601/1‑149 (MQ=255)
         cAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTa                                                                                                                                              <  7:224909/148‑1 (MQ=255)
         cAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTa                                                                                                                                              >  1:311309/1‑148 (MQ=255)
          aggaggAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATc                                                                                                                                            <  1:101856/149‑1 (MQ=255)
          aggaggAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATc                                                                                                                                            <  1:303767/149‑1 (MQ=255)
              ggaACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGc                                                                                                                                        <  5:338387/149‑1 (MQ=255)
              ggaACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGACGAGCCCATCAAAGAACTGGAGGCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGc                                                                                                                                        <  6:163513/149‑1 (MQ=255)
                   acaTGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAAc                                                                                                                                   >  1:303430/1‑149 (MQ=255)
                               ttGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTgc                                                                                                                       >  8:55494/1‑149 (MQ=255)
                                    aTTATGGGACCGCATGGCGTCTTTTATAAAGATGCGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGt                                                                                                                  >  4:27672/1‑149 (MQ=255)
                                    aTTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGt                                                                                                                  >  2:311824/1‑149 (MQ=255)
                                        tGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGAtt                                                                                                              >  7:205649/1‑149 (MQ=255)
                                           gACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGAttttt                                                                                                           <  5:182680/149‑1 (MQ=255)
                                             ccGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCCGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGa                                                                                                         <  3:315601/149‑1 (MQ=255)
                                                cATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTg                                                                                                      <  1:336133/149‑1 (MQ=255)
                                                    gCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGAt                                                                                                  >  6:203088/1‑149 (MQ=255)
                                                     cGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATg                                                                                                 >  6:361805/1‑149 (MQ=255)
                                                      gTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCCCAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGAt                                                                                                  >  2:36988/1‑147 (MQ=255)
                                                      gTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGAt                                                                                                  >  1:240998/1‑147 (MQ=255)
                                                      gTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGa                                                                                                <  1:175300/149‑1 (MQ=255)
                                                       tCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGAt                                                                                                  >  1:161338/1‑146 (MQ=255)
                                                       tCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGa                                                                                                <  2:161338/148‑1 (MQ=255)
                                                         ttttATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTa                                                                                             >  6:26722/1‑149 (MQ=255)
                                                         ttttATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGAGGCGCAAGGCTTTCAGATTATCTCGCCACAAAACAGCGTTGATTTGCCGAAAATTATCGAGCATAACCATCGAATTTGCGGCGTGATTTTTGACTGGCATGAGt                                                                                              >  3:328983/1‑148 (MQ=255)
                                                           ttATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTAc                                                                                            <  1:311824/148‑1 (MQ=255)
                                                           taataAAGATGAGCCAATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGGCACAAAACAGCCTTGATTTGCTGAAATTTATCTAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACAGGGATGAGTACa                                                                                           >  4:109929/3‑149 (MQ=255)
                                                            tataAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAg                                                                                          >  1:148369/1‑149 (MQ=255)
                                                             ataAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTgcg                                                                                                                      >  6:319661/1‑120 (MQ=255)
                                                              taAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGta                                                                                        >  6:175526/1‑148 (MQ=255)
                                                              taAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAAGTTGCGGCGTGATTTTTGACTGGGATGAGTGCAGtc                                                                                        >  8:278496/1‑149 (MQ=255)
                                                              taAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGGTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGTCGTGAttttt                                                                                                           >  5:73773/1‑130 (MQ=255)
                                                              taAAGATGAGCCCATCAAACAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACATCGTTGATTTGCTGAAATTTATCGAGCATAAACATCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACatt                                                                                         >  6:131351/1‑146 (MQ=255)
                                                              caaaCATGCGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAATGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGtc                                                                                        <  2:303430/148‑1 (MQ=255)
                                                                aaGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGtctc                                                                                      <  2:240998/149‑1 (MQ=255)
                                                                   aTGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGAt                                                                                   <  5:26722/149‑1 (MQ=255)
                                                                       gCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGa                                                                                                                                          <  1:102541/90‑1 (MQ=255)
                                                                       gCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGa                                                                                                                                          >  2:102541/1‑90 (MQ=255)
                                                                        cccATCAAAGAACTGGAGTCGGCGCTGGTGGCGAAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATg                                                                              >  2:99604/1‑149 (MQ=255)
                                                                         ccATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCTTGATTTTTGACTGGTAt                                                                                                  >  3:131258/1‑128 (MQ=255)
                                                                         ccATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGAt                                                                                                  <  4:131258/128‑1 (MQ=255)
                                                                         ccATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGt                                                                             <  5:361805/149‑1 (MQ=255)
                                                                          cATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATg                                                                              >  4:129027/1‑147 (MQ=255)
                                                                          cATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTa                                                                            <  2:311309/149‑1 (MQ=255)
                                                                                aGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGata                                                                      <  2:148369/149‑1 (MQ=255)
                                                                                          tCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCt                                                            >  5:144391/1‑149 (MQ=255)
                                                                                            ggCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCtt                                                           >  5:179955/1‑148 (MQ=255)
                                                                                              ctcTGGTTGCGCAAGGCTTTCAGATTATCTGGCCACAATACAGCGTTGATTTGCCGTAATTTATCGAGCATAACAGTCGAATTTGCGGCGTGATTTTTGACTGGGTTGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAAt                                                        <  4:328983/147‑1 (MQ=255)
                                                                                              cgcTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAAt                                                        <  5:175526/149‑1 (MQ=255)
                                                                                                   gtggCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGaa                                                     <  1:36988/147‑1 (MQ=255)
                                                                                                        gcAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTcc                                              >  4:179381/1‑149 (MQ=255)
                                                                                                               tttCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTa                                       >  3:221786/1‑149 (MQ=255)
                                                                                                               tttCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTa                                       >  5:266256/1‑149 (MQ=255)
                                                                                                               tttCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTa                                       >  4:169988/1‑149 (MQ=255)
                                                                                                               tttCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTa                                       >  2:80711/1‑149 (MQ=255)
                                                                                                               tttCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTa                                       >  4:318012/1‑149 (MQ=255)
                                                                                                                      ttATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTATGCCTtc                                <  6:144391/149‑1 (MQ=255)
                                                                                                                         tCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTATGCCTtcat                              >  3:39449/1‑148 (MQ=255)
                                                                                                                             gCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTATGCCTtcatc                             <  4:153277/145‑1 (MQ=255)
                                                                                                                             gCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTATGCCTtcatc                             <  2:155413/145‑1 (MQ=255)
                                                                                                                               cacaAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACAcc                       >  7:228934/1‑149 (MQ=255)
                                                                                                                               cacaAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACAcc                       >  6:326281/1‑149 (MQ=255)
                                                                                                                                acaAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACAcc                       >  7:58145/1‑148 (MQ=255)
                                                                                                                                  aaaaCAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACACCCAc                    <  7:55494/149‑1 (MQ=255)
                                                                                                                                    aaCAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATatcaatc                                                                <  8:324373/103‑1 (MQ=255)
                                                                                                                                    aaCAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATatcaatc                                                                >  7:324373/1‑103 (MQ=255)
                                                                                                                                    aaCAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTAt                                      <  4:322036/129‑1 (MQ=255)
                                                                                                                                    aaCAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTAt                                      >  3:322036/1‑129 (MQ=255)
                                                                                                                                    aaCAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTc                  >  6:289263/1‑149 (MQ=255)
                                                                                                                                         cGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTcc                                              <  1:385669/116‑1 (MQ=255)
                                                                                                                                         cGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTcc                                              >  2:385669/1‑116 (MQ=255)
                                                                                                                                           ttGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACgatg           <  8:228934/149‑1 (MQ=255)
                                                                                                                                           ttGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACgatg           >  2:254707/1‑149 (MQ=255)
                                                                                                                                              aTTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACgatggat        <  6:179955/149‑1 (MQ=255)
                                                                                                                                                 tGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGt      <  4:221786/148‑1 (MQ=255)
                                                                                                                                                  gCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGt      >  1:99481/1‑147 (MQ=255)
                                                                                                                                                    tGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGat                                                                       <  5:289872/80‑1 (MQ=255)
                                                                                                                                                    tGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGat                                                                       >  6:289872/1‑80 (MQ=255)
                                                                                                                                                    tGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTc     >  7:245643/1‑146 (MQ=255)
                                                                                                                                                    tGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTc     <  8:245643/146‑1 (MQ=255)
                                                                                                                                                    tGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGc  >  7:91200/1‑149 (MQ=255)
                                                                                                                                                    tGAAATTTATCGAGCATAACCCTCGAATATGCGGAGTGAGTTTTGACTGGGAAGAGTACATTCTCGAATTATGTAGCGATTTCAATCAGCATAATGAATATCTCCCGCTTAATGAATTCAACAACACACACTCCACCAAGCATGTCAGc  >  8:141832/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGGATTTTCCAGGAGGAACACATGAACATCATTGCCATTATGGGACCGCATGGCGTCTTTTATAAAGATGAGCCCATCAAAGAACTGGAGTCGGCGCTGGTGGCGCAAGGCTTTCAGATTATCTGGCCACAAAACAGCGTTGATTTGCTGAAATTTATCGAGCATAACCCTCGAATTTGCGGCGTGATTTTTGACTGGGATGAGTACAGTCTCGATTTATGTAGCGATATCAATCAGCTTAATGAATATCTCCCGCTTTATGCCTTCATCAACACCCACTCGACGATGGATGTCAGC  >  NC_000913/209658‑209954

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: