Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,107,001 T→C 6.0% intergenic (+31/+19) 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,107,0010TC 6.0% 479.3 / 9.2 156intergenic (+31/+19)nupG/speCnucleoside transporter/ornithine decarboxylase, constitutive
Reads supporting (aligned to +/- strand):  ref base T (75/71);  new base C (3/7);  total (78/78)
Fisher's exact test for biased strand distribution p-value = 3.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 162 total aligned reads displayed.

GCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAGGAAATAACGTT  >  NC_000913/3106853‑3107147
                                                                                                                                                    |                                                                                                                                                  
gCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGt                                                                                                                                                    >  1:422071/1‑149 (MQ=255)
   gACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGAc                                                                                                                                                 <  4:184094/149‑1 (MQ=255)
   gACCGTATGGTTGATCTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGAc                                                                                                                                                 >  3:294339/1‑149 (MQ=255)
        tATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGtt                                                                                                                                            <  3:417035/149‑1 (MQ=255)
         aTGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCACCTTCTGGCGACCCGttt                                                                                                                                           <  1:175244/149‑1 (MQ=255)
          tGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGtttt                                                                                                                                          <  3:1497504/149‑1 (MQ=255)
            gTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGtttttt                                                                                                                                        >  2:702212/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)
              tGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCACCTTCTGGCGACCCGtttttctt                                                                                                                                      <  1:222170/149‑1 (MQ=39)
               gATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTa                                                                                                                                     <  4:221260/149‑1 (MQ=255)
               gATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTa                                                                                                                                     <  4:251254/149‑1 (MQ=255)
               gATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTa                                                                                                                                     <  4:864874/149‑1 (MQ=255)
                aTTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCTCTAATGTCGCAAAGATAAACGGGTCACCTTCTGGCGTCCCGtttttctttg                                                                                                                                    <  2:780343/149‑3 (MQ=255)
                aTTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTAt                                                                                                                                    <  2:747984/149‑1 (MQ=255)
                 ttttCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTAtt                                                                                                                                   <  2:977389/149‑1 (MQ=255)
                  tttCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTc                                                                                                                                  <  4:697166/149‑1 (MQ=255)
                   ttCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCt                                                                                                                                 <  1:1038812/149‑1 (MQ=255)
                     cGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTa                                                                                                                               <  3:710489/149‑1 (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                                                                                                                                      <  3:880953/140‑1 (MQ=38)
                       cTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCACCTTCTGGCGACCCGtttttctt                                                                                                                                      >  4:880953/1‑140 (MQ=38)
                            tACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAac                                                                                                                        <  1:702212/149‑1 (MQ=255)
                            tACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAac                                                                                                                        <  2:56700/149‑1 (MQ=255)
                            tACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAac                                                                                                                        >  2:398202/1‑149 (MQ=255)
                               tCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGa                                                                                                                                                  <  1:825044/120‑1 (MQ=255)
                               tCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGa                                                                                                                                                  >  2:825044/1‑120 (MQ=255)
                               tCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAAcaca                                                                                                                     >  2:1612751/1‑149 (MQ=255)
                               tCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAAcaca                                                                                                                     >  1:1447156/1‑149 (MQ=255)
                               tCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCACCTTCTGGCGACCCGtttttctttgc                                                                                                                                   >  3:1529129/1‑132 (MQ=38)
                               tCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGACGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAAcaca                                                                                                                     >  4:331209/1‑149 (MQ=255)
                                        cTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTAc                                                                                                            >  3:973751/1‑149 (MQ=255)
                                           gCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAAc                                                                                                         >  4:1140389/1‑149 (MQ=255)
                                             cTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGt                                                                                                                                                    >  3:1659939/1‑104 (MQ=255)
                                              ttCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGt                                                                                                      >  1:1529215/1‑149 (MQ=255)
                                              ttCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGt                                                                                                      >  3:606210/1‑149 (MQ=255)
                                               tCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGtt                                                                                                     <  3:299108/149‑1 (MQ=255)
                                                   gTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGt                                                                                                                                                    >  2:101823/1‑98 (MQ=255)
                                                   gTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCAc                                                                                                 <  4:1394880/149‑1 (MQ=255)
                                                     tCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGc                                                                                               >  3:107647/1‑149 (MQ=255)
                                                     tCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGc                                                                                               <  3:1162069/149‑1 (MQ=255)
                                                     tCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGc                                                                                               <  4:1380008/149‑1 (MQ=255)
                                                                ttCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCg                                                                                    <  2:377238/149‑1 (MQ=255)
                                                                 tCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGcc                                                                                              <  4:1359967/138‑1 (MQ=255)
                                                                 tCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGcc                                                                                              >  3:1359967/1‑138 (MQ=255)
                                                                   aaaTATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGtt                                                                                 >  4:1263983/1‑149 (MQ=255)
                                                                   aaaTATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGtt                                                                                 >  2:350891/1‑149 (MQ=255)
                                                                       ataAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGc                                                                             <  2:284915/149‑1 (MQ=255)
                                                                        taAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCt                                                                            <  1:36969/149‑1 (MQ=255)
                                                                         aaaCACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCta                                                                           >  1:393593/1‑149 (MQ=255)
                                                                          aaCACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCtat                                                                          >  2:250706/1‑149 (MQ=255)
                                                                          aaCACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCtat                                                                          <  2:1447156/149‑1 (MQ=255)
                                                                               gTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAaca                                                                     >  3:1004428/1‑149 (MQ=255)
                                                                                ttCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAacac                                                                    <  1:980914/149‑1 (MQ=255)
                                                                                   gtgtCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCtt                                                                 <  1:1669571/149‑1 (MQ=255)
                                                                                   gtgtCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCtt                                                                 <  1:1060328/149‑1 (MQ=255)
                                                                                     gtCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGc                                                               <  4:288028/149‑1 (MQ=255)
                                                                                      tCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTc                                                                                                                                  >  4:639786/1‑81 (MQ=255)
                                                                                      tCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCa                                                              <  2:1529215/149‑1 (MQ=255)
                                                                                         cGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCt                                                           >  4:1115418/1‑149 (MQ=255)
                                                                                          gACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTc                                                          >  2:1720901/1‑149 (MQ=255)
                                                                                           aCAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTcc                                                         <  4:3976/149‑1 (MQ=255)
                                                                                           aCAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTcc                                                         <  4:107647/149‑1 (MQ=255)
                                                                                               gCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCg                                                     >  3:558439/1‑149 (MQ=255)
                                                                                                 acacaGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAGCCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGaa                                                   <  3:346618/149‑1 (MQ=255)
                                                                                                 acacaGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGaa                                                   >  2:1682219/1‑149 (MQ=255)
                                                                                                   acaGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGaaaa                                                 >  3:892409/1‑149 (MQ=255)
                                                                                                   acaGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGaaaa                                                 >  4:1634010/1‑149 (MQ=255)
                                                                                                    caGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAAt                                                >  1:1556566/1‑149 (MQ=255)
                                                                                                      gACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCACCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATcc                                              <  1:1011406/149‑1 (MQ=255)
                                                                                                               ccACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTa                                                                                                                               <  4:325177/59‑1 (MQ=255)
                                                                                                               ccACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGcaacaa                                     >  1:1350305/1‑149 (MQ=255)
                                                                                                               ccACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGcaacaa                                     >  2:413103/1‑149 (MQ=255)
                                                                                                                cACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACaaa                                    <  4:717560/149‑1 (MQ=255)
                                                                                                                cACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACaaa                                    <  4:1004428/149‑1 (MQ=255)
                                                                                                                     aTTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGc                                                      >  4:1447713/1‑126 (MQ=255)
                                                                                                                     aTTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCAc                               <  2:1294711/149‑1 (MQ=255)
                                                                                                                     aTTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGATTCGCTATAAACACCTTGCAGCTCCGGc                                                      <  3:1447713/126‑1 (MQ=255)
                                                                                                                      ttACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCAcc                              <  1:273519/149‑1 (MQ=255)
                                                                                                                            aaaGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCcttct                        >  4:1499380/1‑149 (MQ=255)
                                                                                                                             aaGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCAc                                                                                                 <  4:744849/75‑1 (MQ=255)
                                                                                                                             aaGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCAc                                                                                                 >  3:744849/1‑75 (MQ=255)
                                                                                                                             aaGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCcttctt                       >  2:570438/1‑149 (MQ=255)
                                                                                                                             aaGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCcttctt                       >  2:1229568/1‑149 (MQ=255)
                                                                                                                              aGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCttcttc                      <  4:12854/149‑1 (MQ=255)
                                                                                                                               gAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTcc                     >  1:639843/1‑149 (MQ=255)
                                                                                                                               aaaaaaCGGGTCACCTTCTGGCGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGc                                                               <  4:1312619/106‑1 (MQ=37)
                                                                                                                               aaaaaaCGGGTCACCTTCTGGCGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTcc                     >  3:997789/2‑149 (MQ=255)
                                                                                                                                aaaaaCGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCa                    <  1:751870/149‑1 (MQ=255)
                                                                                                                                aaaaaCGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCa                    >  4:1516922/1‑149 (MQ=255)
                                                                                                                                aaaaaCGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCa                    <  3:326949/149‑1 (MQ=255)
                                                                                                                                   aaCGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGaaa                                                  <  2:669394/116‑1 (MQ=255)
                                                                                                                                        gTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAgg            <  4:721108/149‑1 (MQ=255)
                                                                                                                                           gCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAGGaaa         >  2:705850/1‑149 (MQ=255)
                                                                                                                                                  ggTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAGGAAATAACGtt  <  3:1140389/149‑1 (MQ=255)
                                                                                                                                                  ggTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAGGAAATAACGtt  >  1:1017841/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAGGAAATAACGTT  >  NC_000913/3106853‑3107147

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