Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,108,589 G→T 7.1% N189K (AAC→AAA speC ← ornithine decarboxylase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,108,5890GT 7.1% 289.3 / 10.1 85N189K (AAC→AAAspeCornithine decarboxylase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (39/40);  new base T (3/3);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.55e-01

CGTGAAGCTTCCAGATAGACCGGCGTCGCCCCCGCCTGAATCAGCGCGCCGTGATGATTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAATTTTACGTCAGCGTTACACATATCGGCGCGAAAGACGTTCTCACCAAAGA  >  NC_000913/3108443‑3108734
                                                                                                                                                  |                                                                                                                                                 
cGTGAAGCTTCCAGATAGACCGGCGTCGCCCCCGCCTGAATCAGCGCGCCGTGATGATTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGtt                                                                                                                                                 <  1:894578/149‑1 (MQ=255)
 gTGAAGCTTCCAGATAGACCGGCGTCGCCCCCGCCTGAATCAGCGCGCCGTGATGATTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTc                                                                                                                                                >  1:513290/1‑149 (MQ=255)
  tGAAGCTTCCAGATAGACCGGCGTCGCCCCCGCCTGAATCAGCGCGCCGTGATGATTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCa                                                                                                                                               <  3:185097/149‑1 (MQ=255)
  tGAAGCTTCCAGATAGACCGGCGTCGCCCCCGCCTGAATCAGCGCGCCGTGATGATTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCa                                                                                                                                               <  3:155739/149‑1 (MQ=255)
  tGAAGCTTCCAGATAGACCGGCGTCGCCCCCGCCTGAATCAGCGCGCCGTGATGATTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCa                                                                                                                                               >  4:475030/1‑149 (MQ=255)
          ccAGATAGACCGGCGTCGCCCCCGCCTGAATCAGCGCGCCGTGATGATTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACaaaa                                                                                                                                       >  3:166029/1‑149 (MQ=255)
           cAGATAGACCGGCGTCGCCCCCGCCTGAATCAGCGCGCCGTGATGATTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAAt                                                                                                                                      >  4:345572/1‑149 (MQ=255)
            aGATAGACCGGCGTCGCCCCCGCCTGAATCAGCGCGCCGTGATGATTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATa                                                                                                                                     >  3:628471/1‑149 (MQ=255)
             gATAGACCGGCGTCGCCCCCGCCTGAATCAGCGCGCCGTGATGATTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAg                                                                                                                                    <  1:440520/149‑1 (MQ=255)
             gATAGACCGGCGTCGCCCCCGCCTGAATCAGCGCGCCGTGATGATTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAg                                                                                                                                    >  2:501104/1‑149 (MQ=255)
             gATAGACCGGCGTCGCCCCCGCCTGAATCAGCGCGCCGTGATGATTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAg                                                                                                                                    <  2:672911/149‑1 (MQ=255)
               tAGACCGGCGTCGCCCCCGCCTGAATCAGCGCGCCGTGATGATTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGt                                                                                                                                  <  4:565275/149‑1 (MQ=255)
                aGACCGGCGTCGCCCCCGCCTGAATCAGCGCGCCGTGATGATTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGtt                                                                                                                                 <  2:183295/149‑1 (MQ=255)
                          cGCCCCCGCCTGAATCAGCGCGCCGTGATGATTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCAt                                                                                                                       <  3:268267/149‑1 (MQ=255)
                            cccccGCCTGAATCAGCGCGCCGTGATGATTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGa                                                                                                                     >  2:678773/1‑149 (MQ=255)
                             ccccGCCTGAATCAGCGCGCCGTGATGATTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGaa                                                                                                                    >  2:157624/1‑149 (MQ=255)
                              cccGCCTGAATCAGCGCGCCGTGATGATTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGaaa                                                                                                                   <  1:764766/149‑1 (MQ=255)
                                cGCCTGAATCAGCGCGCCGTGATGATTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGa                                                                                                                 <  4:219043/149‑1 (MQ=255)
                                   cTGAATCAGCGCGCCGTGATGATTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAACGACtt                                                                                                              <  4:470553/149‑1 (MQ=255)
                                   cTGAATCAGCGCGCCGTGATGATTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACtt                                                                                                              >  4:178150/1‑149 (MQ=255)
                                    tGAATCAGCGCGCCGTGATGATTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACttt                                                                                                             <  3:208184/149‑1 (MQ=255)
                                       aTCAGCGCGCCGTGATGATTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGtt                                                                                                                                 >  1:315883/1‑126 (MQ=255)
                                       aTCAGCGCGCCGTGATGATTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGtt                                                                                                                                 <  2:315883/126‑1 (MQ=255)
                                                 cGTGATGATTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTc                                                                                                <  1:504099/149‑1 (MQ=255)
                                                   tgatgaTTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGa                                                                                                                 <  3:141757/130‑1 (MQ=255)
                                                   tgatgaTTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGa                                                                                                                 >  4:141757/1‑130 (MQ=255)
                                                        atTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCTTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTg                                                                                                            >  3:238730/1‑130 (MQ=255)
                                                        atTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCTTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTg                                                                                                            >  3:372594/1‑130 (MQ=255)
                                                        atTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCTTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTg                                                                                                            <  4:238730/130‑1 (MQ=255)
                                                        atTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCTTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTg                                                                                                            <  2:45624/130‑1 (MQ=255)
                                                        atTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCTTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTg                                                                                                            >  1:45624/1‑130 (MQ=255)
                                                        atTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCTTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTg                                                                                                            <  4:372594/130‑1 (MQ=255)
                                                            gACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCttt                                                                                     <  1:157624/149‑1 (MQ=255)
                                                               ttATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGc                                                                                  >  2:87695/1‑149 (MQ=255)
                                                                tATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGcc                                                                                 <  3:597331/149‑1 (MQ=255)
                                                                 aTGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCCAAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCg                                                                                >  4:205831/1‑149 (MQ=255)
                                                                   ggTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGAt                                                                              <  1:607951/149‑1 (MQ=255)
                                                                       gttACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCtt                                                                          <  4:628471/149‑1 (MQ=255)
                                                                          aCGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCAt                                                                       <  4:15615/149‑1 (MQ=255)
                                                                             gTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGaa                                                                    >  4:257103/1‑149 (MQ=255)
                                                                              tCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAAt                                                                   >  1:786782/1‑149 (MQ=255)
                                                                                gAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATaa                                                                 >  2:462800/1‑149 (MQ=255)
                                                                                gAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATaa                                                                 >  3:430811/1‑149 (MQ=255)
                                                                                  agagCACCAGATCGCCACGCGTTAACAGCGCATTCGTTACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGTTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCTGATCCTTCCTGAATAAGc                                                               <  4:256531/149‑1 (MQ=255)
                                                                                   gagCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCa                                                              >  4:451446/1‑149 (MQ=255)
                                                                                   gagCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCa                                                              <  4:712829/149‑1 (MQ=255)
                                                                                     gCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGa                                                            >  1:843359/1‑149 (MQ=255)
                                                                                     gCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGa                                                            <  2:152631/149‑1 (MQ=255)
                                                                                        ccAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCg                                                         >  1:368880/1‑149 (MQ=255)
                                                                                        ccAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCg                                                         >  1:191980/1‑149 (MQ=255)
                                                                                          aGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGcc                                                       <  1:251763/149‑1 (MQ=255)
                                                                                            aTCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCa                                                     >  2:149455/1‑149 (MQ=255)
                                                                                               gCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAAtt                                                  >  2:751838/1‑149 (MQ=255)
                                                                                               gCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAAtt                                                  >  4:619208/1‑149 (MQ=255)
                                                                                                ccACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAAttt                                                 >  1:650548/1‑149 (MQ=255)
                                                                                                 cACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAAtttt                                                >  4:579232/1‑149 (MQ=255)
                                                                                                  aCGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAATTTTa                                               <  3:205831/149‑1 (MQ=255)
                                                                                                   cgcgTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAATTTTAc                                              <  2:251780/149‑1 (MQ=255)
                                                                                                   cgcgTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAATTTTAc                                              <  1:17259/149‑1 (MQ=255)
                                                                                                       ttAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAATTTTACGTCa                                          >  1:174835/1‑149 (MQ=255)
                                                                                                         aaCAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAATTTTACGTCAGc                                        <  4:556225/149‑1 (MQ=255)
                                                                                                             gcgcATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAATTTTACGTCAGCGTTa                                    >  3:331816/1‑149 (MQ=255)
                                                                                                               gcATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAATTTTACGTCAGCGTTaca                                  <  4:713671/149‑1 (MQ=255)
                                                                                                                cATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAATTTTACGTCAGCGTTacac                                 <  3:257103/149‑1 (MQ=255)
                                                                                                                        accacTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAATTTTACGTCAGCGTTACACATATCGgc                         >  1:746559/1‑149 (MQ=255)
                                                                                                                        accacTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAATTTTACGTCAGCGTTACACATATCGgc                         <  4:430811/149‑1 (MQ=255)
                                                                                                                          cacTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCCATTTTACGTCAGCGTTACACATATCGgcgc                       >  2:422542/1‑149 (MQ=255)
                                                                                                                          cacTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAATTTTACGTCAGCGTTACACATATCGgcgc                       >  1:133090/1‑149 (MQ=255)
                                                                                                                          cacTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAATTTTACGTCAGCGTTACACATATCGgcgc                       >  2:58910/1‑149 (MQ=255)
                                                                                                                          cacTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAATTTTACGTCAGCGTTACACATATCGgc                         >  4:6967/1‑147 (MQ=255)
                                                                                                                          cacTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAATTTTACGTCAGCGTTACACATATCGgc                         <  3:6967/147‑1 (MQ=255)
                                                                                                                                 ttCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGCTCCTTCATGAATAAGCAGATCGCCCAATTTTACGTCAGCGTTACACATATCGGCGCGAAAGAc                >  3:542834/1‑149 (MQ=255)
                                                                                                                                 ttCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAAt                                                                   <  4:419733/98‑1 (MQ=255)
                                                                                                                                 ttCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAAt                                                                   >  3:419733/1‑98 (MQ=255)
                                                                                                                                 ttCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAATTTTACGTCAGCGTTACACATATCGGCGCGAAAGAc                >  2:10061/1‑149 (MQ=255)
                                                                                                                                 ttCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAATTTTACGTCAGCGTTACACATATCGGCGCGAAAGAc                <  1:501104/149‑1 (MQ=255)
                                                                                                                                  tCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAATTTTACGTCAGCGTTACACATATCGGCGCGAAAGACg               <  1:678773/149‑1 (MQ=255)
                                                                                                                                    gCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGTCCTTGGCTGCCAATTTCTGCGCATCTTTCACCGCTCCTTCATGAATAAGCAGAACCCCCAATTTTACGTTAGCGTTACACATAACGGCGCGAAAGACGta             <  2:513290/149‑2 (MQ=255)
                                                                                                                                    gCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAATTTTACGTCAGCGTTACACATATCGGCGCGAAAGACGtt             <  3:451446/149‑1 (MQ=255)
                                                                                                                                        ccGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAATTTTACGTCAGCGTTACACATATCGGCGCGAAAGACGTTCTCa         <  3:345572/149‑1 (MQ=255)
                                                                                                                                         cGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGAAGATCGCCCAATTTTACGTCAGCGTTACACATATCGGCGCGAAAGACGTTCTCAc        >  2:4984/1‑149 (MQ=255)
                                                                                                                                              tgCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAATTTTACGTCAGCGTTACACATATCGGCGCGAAAGACGTTCTCACCAAAg   <  4:617260/149‑1 (MQ=255)
                                                                                                                                               gCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAATTTTACGTCAGCGTTACACATATCGGCGCGAAAGACGTTCTCACCAAAGa  >  3:17775/1‑149 (MQ=255)
                                                                                                                                                ccGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGAt                                                           >  2:103416/1‑91 (MQ=255)
                                                                                                                                                ccGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGAt                                                           <  1:103416/91‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CGTGAAGCTTCCAGATAGACCGGCGTCGCCCCCGCCTGAATCAGCGCGCCGTGATGATTCGACTTATGGTTGTTACGGTCGAAGAGCACCAGATCGCCACGCGTTAACAGCGCATTCGTCACCACTTTATTCGCTGCCGATGTGCCGTTCAGCACAAAATAGGTTTTATCGGCATGAAAGACTTTGGCTGCGAATTTCTGCGCATCTTTCGCCGATCCTTCATGAATAAGCAGATCGCCCAATTTTACGTCAGCGTTACACATATCGGCGCGAAAGACGTTCTCACCAAAGA  >  NC_000913/3108443‑3108734

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