Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,106,996 G→T 6.2% intergenic (+26/+24) nupG → / ← speC nucleoside transporter/ornithine decarboxylase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,106,9960GT 6.2% 558.4 / 19.5 163intergenic (+26/+24)nupG/speCnucleoside transporter/ornithine decarboxylase, constitutive
Reads supporting (aligned to +/- strand):  ref base G (79/72);  new base T (3/7);  total (82/81)
Fisher's exact test for biased strand distribution p-value = 2.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

Only 100 of 163 total aligned reads displayed.

CTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAGGAAATA  >  NC_000913/3106850‑3107142
                                                                                                                                                  |                                                                                                                                                  
atGGCAGTCCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGCCGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGaa                                                                                                                                                  <  2:837039/148‑1 (MQ=255)
 tGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAg                                                                                                                                                 >  2:1123615/1‑149 (MQ=255)
   gCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGt                                                                                                                                               <  3:1539485/149‑1 (MQ=255)
      gACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGAc                                                                                                                                            <  4:184094/149‑1 (MQ=255)
      gACCGTATGGTTGATCTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGAc                                                                                                                                            >  3:294339/1‑149 (MQ=255)
            aTGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCACCTTCTGGCGACCCGttt                                                                                                                                      <  1:175244/149‑1 (MQ=255)
             tGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGtttt                                                                                                                                     <  3:1497504/149‑1 (MQ=255)
               gTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGtttttt                                                                                                                                   >  2:702212/1‑149 (MQ=255)
               gTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGtttttt                                                                                                                                   >  4:1317633/1‑149 (MQ=255)
               gTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGtttttt                                                                                                                                   >  3:697166/1‑149 (MQ=255)
               gTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGtttttt                                                                                                                                   <  1:25416/149‑1 (MQ=255)
                ttGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGttttttt                                                                                                                                  >  4:867636/1‑149 (MQ=255)
                  gATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTa                                                                                                                                <  4:864874/149‑1 (MQ=255)
                   aTTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCTCTAATGTCGCAAAGATAAACGGGTCACCTTCTGGCGTCCCGtttttctttg                                                                                                                               <  2:780343/149‑3 (MQ=255)
                   aTTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTAt                                                                                                                               <  2:747984/149‑1 (MQ=255)
                     tttCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTc                                                                                                                             <  4:697166/149‑1 (MQ=255)
                      ttCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCt                                                                                                                            >  2:175983/1‑149 (MQ=255)
                        cGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTa                                                                                                                          <  3:1396950/149‑1 (MQ=255)
                         gCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCTACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTAc                                                                                                                         <  4:1325818/149‑1 (MQ=255)
                         gCTGGTTACCCCGTGGTTCGGGAATTCGCGTTCATGGCGTTGTTCAAATATAATCACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCGCTATTTTCGCAAAGAAAAACGGGTCCCCTCAAGGTGACCCGGTTTTTTTATTCtttc                                                                                                                         <  2:1062237/149‑3 (MQ=255)
                          cTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACt                                                                                                                        >  4:346618/1‑149 (MQ=255)
                          cTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCACCTTCTGGCGACCCGtttttctt                                                                                                                                 >  4:880953/1‑140 (MQ=38)
                           tGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACtt                                                                                                                       >  4:61218/1‑149 (MQ=255)
                              ttACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCaa                                                                                                                    >  2:724458/1‑149 (MQ=255)
                               tACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAac                                                                                                                   >  2:398202/1‑149 (MQ=255)
                               tACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAac                                                                                                                   <  1:702212/149‑1 (MQ=255)
                               tACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAac                                                                                                                   <  2:56700/149‑1 (MQ=255)
                                  tCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGaa                                                                                                                                                  <  3:468656/115‑1 (MQ=255)
                                  tCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGaa                                                                                                                                                  >  4:468656/1‑115 (MQ=255)
                                  tCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGa                                                                                                                                             <  1:825044/120‑1 (MQ=255)
                                  tCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAAcaca                                                                                                                >  1:1543614/1‑149 (MQ=255)
                                  tCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAAcaca                                                                                                                >  2:36969/1‑149 (MQ=255)
                                  tCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCACCTTCTGGCGACCCGtttttctttgc                                                                                                                              <  4:1529129/135‑4 (MQ=38)
                                          tCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTa                                                                                                        >  3:1106044/1‑149 (MQ=255)
                                           cTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTAc                                                                                                       >  3:973751/1‑149 (MQ=255)
                                            tGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACa                                                                                                      >  4:1286144/1‑149 (MQ=255)
                                                cTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGt                                                                                                                                               <  4:1659939/104‑1 (MQ=255)
                                                cTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCg                                                                                                  >  1:261379/1‑149 (MQ=255)
                                                 ttCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGt                                                                                                 >  1:1529215/1‑149 (MQ=255)
                                                 ttCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGt                                                                                                 >  3:606210/1‑149 (MQ=255)
                                                    gcgTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTc                                                                                              <  2:56225/149‑1 (MQ=255)
                                                    gcgTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACACAACCGTACAACCGTTTc                                                                                              <  1:359341/149‑1 (MQ=255)
                                                     cgTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCa                                                                                             <  1:175983/149‑1 (MQ=255)
                                                      gTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGt                                                                                                                                               <  1:101823/98‑1 (MQ=255)
                                                      gTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCAc                                                                                            <  4:1394880/149‑1 (MQ=255)
                                                        tCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGc                                                                                          <  2:637144/149‑1 (MQ=255)
                                                            ggCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCTACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATc                                                                                      <  1:1440931/149‑1 (MQ=255)
                                                                   ttCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCg                                                                               <  2:377238/149‑1 (MQ=255)
                                                                    tCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGcc                                                                                         >  3:1359967/1‑138 (MQ=255)
                                                                      aaaTATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGtt                                                                            >  4:1263983/1‑149 (MQ=255)
                                                                      aaaTATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGtt                                                                            <  1:398202/149‑1 (MQ=255)
                                                                      aaaTATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGtt                                                                            <  1:724458/149‑1 (MQ=255)
                                                                       aaTATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGGCCCGTCTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGttt                                                                           <  3:976672/149‑1 (MQ=255)
                                                                       aaTATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGttt                                                                           >  2:1669571/1‑149 (MQ=255)
                                                                        atatAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTc                                                                          <  2:1353181/149‑1 (MQ=255)
                                                                        atatAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTc                                                                          <  4:294339/149‑1 (MQ=255)
                                                                          ataAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGc                                                                        <  4:136615/149‑1 (MQ=255)
                                                                          ataAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGc                                                                        <  2:284915/149‑1 (MQ=255)
                                                                           taAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCt                                                                       <  3:1233671/149‑1 (MQ=255)
                                                                            aaaCACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCta                                                                      >  1:393593/1‑149 (MQ=255)
                                                                             aaCACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCtat                                                                     <  2:1447156/149‑1 (MQ=255)
                                                                             aaCACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCtat                                                                     >  3:514816/1‑149 (MQ=255)
                                                                             aaCACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCtat                                                                     >  2:250706/1‑149 (MQ=255)
                                                                              acacGTTCGTGTCCCGACAGGCATACAGACGGTTAGCCACTAATTATGCAAAGAAACACGGGGCGCCCGAAGGTGCCCCGTTTAATTAATCCTAACAACAACACATAACCGAACAACAGCTTCGCGCCTTCCTCATCGTTTTCGCtata                                                                    >  4:1229926/1‑149 (MQ=255)
                                                                              acacGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTAc                                                                                                       <  3:810059/114‑1 (MQ=255)
                                                                                  gTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAaca                                                                >  3:1004428/1‑149 (MQ=255)
                                                                                      gtgtCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCtt                                                            <  1:1060328/149‑1 (MQ=255)
                                                                                      gtgtCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCtt                                                            <  1:1669571/149‑1 (MQ=255)
                                                                                        gtCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGc                                                          <  4:288028/149‑1 (MQ=255)
                                                                                         tCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCa                                                         <  2:1529215/149‑1 (MQ=255)
                                                                                         tCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCa                                                         >  3:274692/1‑149 (MQ=255)
                                                                                           ccGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAACCGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACAGTTTCACGCCATCCGCATCGGATTCGCTATAAACACCTTGCAGc                                                       >  3:372579/1‑149 (MQ=255)
                                                                                           ccGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGc                                                       >  1:993212/1‑149 (MQ=255)
                                                                                            cGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCt                                                      >  4:1115418/1‑149 (MQ=255)
                                                                                              aCAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTcc                                                    <  4:107647/149‑1 (MQ=255)
                                                                                                    acacaGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGaa                                              >  2:1682219/1‑149 (MQ=255)
                                                                                                    acacaGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTCTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGaa                                              >  4:326949/1‑149 (MQ=255)
                                                                                                      acaGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGaaaa                                            >  3:892409/1‑149 (MQ=255)
                                                                                                         gACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCACCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATcc                                         <  1:1011406/149‑1 (MQ=255)
                                                                                                            ggTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCgg                                      >  2:844384/1‑149 (MQ=255)
                                                                                                                  ccACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTa                                                                                                                          >  3:325177/1‑59 (MQ=255)
                                                                                                                   cACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACaaa                               <  4:1004428/149‑1 (MQ=255)
                                                                                                                   cACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACaaa                               <  4:717560/149‑1 (MQ=255)
                                                                                                                               aaaGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCcttct                   >  3:1476865/1‑149 (MQ=255)
                                                                                                                               aaaGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCcttct                   >  4:1432238/1‑149 (MQ=255)
                                                                                                                                aaGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCAc                                                                                            >  3:744849/1‑75 (MQ=255)
                                                                                                                                aaGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCAc                                                                                            <  4:744849/75‑1 (MQ=255)
                                                                                                                                aaGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCcttctt                  >  2:1229568/1‑149 (MQ=255)
                                                                                                                                aaGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCcttctt                  >  1:911062/1‑149 (MQ=255)
                                                                                                                                  gAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTcc                >  1:639843/1‑149 (MQ=255)
                                                                                                                                   aaaaaCGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCa               >  4:1516922/1‑149 (MQ=255)
                                                                                                                                   aaaaaCGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCa               <  3:326949/149‑1 (MQ=255)
                                                                                                                                     aaaCGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGCCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCatt             >  3:102940/1‑147 (MQ=255)
                                                                                                                                      aaCGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGaaa                                             <  2:669394/116‑1 (MQ=255)
                                                                                                                                         gggTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCaa         <  3:705849/149‑1 (MQ=255)
                                                                                                                                         gggTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCaa         <  4:1476865/149‑1 (MQ=255)
                                                                                                                                           gTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAgg       <  4:721108/149‑1 (MQ=255)
                                                                                                                                            tCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACACAACCGAACAACAGTTGCACGCTATCCGCTTAGGTTTCGCTATAAACACCTTTCCGCTCCGGCGACAATCCCGGCGACTAATTCACCCCTTCTTCCAGTCCAAGGa      >  4:225476/1‑149 (MQ=255)
                                                                                                                                              gCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAGGaaa    >  4:393539/1‑149 (MQ=255)
                                                                                                                                                cAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAGGAAATa  >  1:126768/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAGGAAATA  >  NC_000913/3106850‑3107142

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