Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 136,062 T→G K135Q (AAG→CAG)  speE ← spermidine synthase (putrescine aminopropyltransferase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913136,0620TG98.7% 280.6 / ‑4.2 74K135Q (AAG→CAG) speEspermidine synthase (putrescine aminopropyltransferase)
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (39/34);  total (40/34)
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.22e-01

GGATTCAGGCAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTTAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACGCAGCATGGCACCGTCGCCGCC  >  NC_000913/135914‑136209
                                                                                                                                                    |                                                                                                                                                   
ggATTCAGGCAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCTCCAGCTGaa                                                                                                                                                   <  2:587083/151‑1 (MQ=255)
  aTTCAGGCAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGaa                                                                                                                                                   >  2:206337/1‑149 (MQ=255)
   ttCAGGCAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAgc                                                                                                                                                >  1:818527/1‑151 (MQ=255)
     cAGGCAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAgcgc                                                                                                                                              <  2:47751/151‑1 (MQ=255)
          aaCGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGAtc                                                                                                                                         >  2:802796/1‑151 (MQ=255)
            cGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGAtcgt                                                                                                                                       >  1:507728/1‑151 (MQ=255)
             gTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGATGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGAtcgtc                                                                                                                                      <  1:614093/151‑1 (MQ=255)
                tGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATcgtcgt                                                                                                                                    >  2:656809/1‑150 (MQ=255)
                   aGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCt                                                                                                                                >  1:63033/1‑151 (MQ=255)
                       ttCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGACTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGAGAACCAGCTGAAAGCGCGGATCGTCGCAGCTACCg                                                                                                                            >  2:541671/1‑151 (MQ=255)
                         cATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGc                                                                                                                          <  2:797954/151‑1 (MQ=255)
                         cATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGc                                                                                                                          <  1:207856/151‑1 (MQ=255)
                          aTAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCg                                                                                                                         >  1:432652/1‑151 (MQ=255)
                           tAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGt                                                                                                                        >  2:649695/1‑151 (MQ=255)
                            aaaaTGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGtt                                                                                                                       <  1:919662/151‑1 (MQ=255)
                            aaaaTGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGtt                                                                                                                       >  2:395654/1‑151 (MQ=255)
                             aaaTGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTa                                                                                                                      <  1:329197/151‑1 (MQ=255)
                             aaaTGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTa                                                                                                                      <  2:791670/151‑1 (MQ=255)
                                 gCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGt                                                                                                                  >  1:426835/1‑151 (MQ=255)
                                 gCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGt                                                                                                                  >  2:50771/1‑151 (MQ=255)
                                   cGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTg                                                                                                                >  2:583753/1‑151 (MQ=255)
                                   cGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTg                                                                                                                <  1:896431/151‑1 (MQ=255)
                                   cGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTg                                                                                                                >  2:797437/1‑151 (MQ=255)
                                    gAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTgg                                                                                                               >  2:534934/1‑151 (MQ=255)
                                     aaGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTggg                                                                                                              >  1:547803/1‑151 (MQ=255)
                                     aaGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTggg                                                                                                              >  1:17162/1‑151 (MQ=255)
                                        tGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGt                                                                                                             >  1:562461/1‑149 (MQ=255)
                                         gAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGt                                                                                                             >  2:257239/1‑148 (MQ=255)
                                         gAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAg                                                                                                           <  1:336360/150‑1 (MQ=255)
                                          aaaaGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGt                                                                                                             >  1:284660/1‑147 (MQ=255)
                                                  ttCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACg                                                                                                  <  2:604858/150‑1 (MQ=255)
                                                   tCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGc                                                                                                >  1:605256/1‑151 (MQ=255)
                                                    cGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCa                                                                                               <  1:897855/151‑1 (MQ=255)
                                                     gCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAg                                                                                              <  1:769292/151‑1 (MQ=255)
                                                       cGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGa                                                                                             >  2:215497/1‑150 (MQ=255)
                                                           aCCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGAt                                                                                        <  1:901034/151‑1 (MQ=255)
                                                                 aGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGaca                                                                                  <  2:363133/151‑1 (MQ=255)
                                                                  ggATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGacac                                                                                 >  1:37576/1‑151 (MQ=255)
                                                                     tCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCg                                                                              >  1:867754/1‑151 (MQ=255)
                                                                        gTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCAt                                                                           <  1:703110/151‑1 (MQ=255)
                                                                          gCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATc                                                                          <  2:132869/150‑1 (MQ=255)
                                                                          gCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCg                                                                         <  1:656824/151‑1 (MQ=255)
                                                                            aGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGAt                                                                       <  2:384000/151‑1 (MQ=255)
                                                                               cGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTc                                                                    <  1:656436/151‑1 (MQ=255)
                                                                                ggAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTTAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTcc                                                                   >  1:693076/1‑151 (MQ=255)
                                                                                ggAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTcc                                                                   >  2:911508/1‑151 (MQ=255)
                                                                                 gagaTAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTcc                                                                   >  2:211159/1‑150 (MQ=255)
                                                                                 gagaTAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGATGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGAttt                                                                     >  2:854400/1‑148 (MQ=255)
                                                                                             tCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTg                                                      >  1:439003/1‑151 (MQ=255)
                                                                                                  ggTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTc                                                  >  1:96447/1‑150 (MQ=255)
                                                                                                     ctggctggTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACg                                              >  1:127685/1‑151 (MQ=255)
                                                                                                      tggctggTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGt                                             >  1:469196/1‑151 (MQ=255)
                                                                                                       ggctggTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAAc                                               <  1:354814/148‑1 (MQ=255)
                                                                                                           ggTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGttttt                                         <  2:200181/150‑1 (MQ=255)
                                                                                                           ggTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGtttt                                          <  1:337807/149‑1 (MQ=255)
                                                                                                            gTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTAt                                       <  2:730412/151‑1 (MQ=255)
                                                                                                                 aTTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCgg                                  <  1:14218/151‑1 (MQ=255)
                                                                                                                 aTTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCgg                                  <  2:432643/151‑1 (MQ=255)
                                                                                                                    aaCGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGt                                <  1:178338/150‑1 (MQ=255)
                                                                                                                     aCGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGtt                               >  1:360175/1‑150 (MQ=255)
                                                                                                                       gAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGAAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACt                            >  2:156353/1‑151 (MQ=255)
                                                                                                                         aaTTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTc                          >  2:629983/1‑151 (MQ=255)
                                                                                                                            tGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCAc                        <  2:759918/150‑1 (MQ=255)
                                                                                                                             gACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACgc                      >  1:565840/1‑151 (MQ=255)
                                                                                                                              aCGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACgc                      >  2:512449/1‑150 (MQ=255)
                                                                                                                                   aTCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACGCAGCATg                >  1:640813/1‑151 (MQ=255)
                                                                                                                                     cgtcgATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACGCCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACGCAGCATgg               <  1:732460/150‑1 (MQ=255)
                                                                                                                                      gtcgATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACGCAGCATGGc              >  1:851758/1‑150 (MQ=255)
                                                                                                                                      gtcgATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACGCAGCATGGCa             <  1:52098/151‑1 (MQ=255)
                                                                                                                                          aTCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACGCAGCATGGCACCGt         <  1:649710/151‑1 (MQ=255)
                                                                                                                                              ccAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACGCAGCATGGCACCGTcgcc     <  2:597370/151‑1 (MQ=255)
                                                                                                                                               cAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACGCAGCATGGCACCGTcgccg    >  1:86660/1‑151 (MQ=255)
                                                                                                                                                 gCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACGCAGCATGGCACCGTcgccgcc  <  1:502763/151‑1 (MQ=255)
                                                                                                                                                 gCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACGCAGCATGGCACCGTcgccgcc  <  1:534949/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGATTCAGGCAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTTAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACGCAGCATGGCACCGTCGCCGCC  >  NC_000913/135914‑136209

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: