Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,106,995 A→T 6.7% intergenic (+25/+25) 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,9950AT 6.7% 568.3 / 20.6 163intergenic (+25/+25)nupG/speCnucleoside transporter/ornithine decarboxylase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (78/73);  new base T (3/8);  total (82/81)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 165 total aligned reads displayed.

CGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAGGAAATA  >  NC_000913/3106847‑3107142
                                                                                                                                                    |                                                                                                                                                    
cGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCa                                                                                                                                                      >  3:1023444/1‑149 (MQ=255)
    tGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAg                                                                                                                                                  >  2:1123615/1‑149 (MQ=255)
      gCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGt                                                                                                                                                >  1:422071/1‑149 (MQ=255)
         gACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑c                                                                                                                                            <  4:184094/149‑1 (MQ=255)
              tATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGtt                                                                                                                                       <  3:417035/149‑1 (MQ=255)
               aTGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCACCTTCTGGCGA‑CCCGttt                                                                                                                                      <  1:175244/149‑1 (MQ=255)
                  gTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGtttttt                                                                                                                                   >  2:702212/1‑149 (MQ=255)
                  gTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGtttttt                                                                                                                                   >  3:697166/1‑149 (MQ=255)
                   ttGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGttttttt                                                                                                                                  >  4:867636/1‑149 (MQ=255)
                    tGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCACCTTCTGGCGA‑CCCGtttttctt                                                                                                                                 <  1:222170/149‑1 (MQ=39)
                     gATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTa                                                                                                                                <  4:864874/149‑1 (MQ=255)
                     gATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTa                                                                                                                                <  4:251254/149‑1 (MQ=255)
                      aTTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCTCTAATGTCGCAAAGATAAACGGGTCACCTTCTGGCGT‑CCCGtttttctttg                                                                                                                               <  2:780343/149‑3 (MQ=255)
                      aTTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTAt                                                                                                                               <  2:747984/149‑1 (MQ=255)
                        tttCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTc                                                                                                                             <  4:697166/149‑1 (MQ=255)
                          tCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCtt                                                                                                                           <  3:966242/149‑1 (MQ=255)
                           cGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTa                                                                                                                          <  3:710489/149‑1 (MQ=255)
                            gCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCTACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTAc                                                                                                                         <  4:1325818/149‑1 (MQ=255)
                             cTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCACCTTCTGGCGA‑CCCGtttttctt                                                                                                                                 <  3:880953/140‑1 (MQ=38)
                             cTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCACCTTCTGGCGA‑CCCGtttttctt                                                                                                                                 >  4:880953/1‑140 (MQ=38)
                              tGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACtt                                                                                                                       >  3:136615/1‑149 (MQ=255)
                              tGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACtt                                                                                                                       >  4:61218/1‑149 (MQ=255)
                                 ttACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCaa                                                                                                                    >  2:724458/1‑149 (MQ=255)
                                  tACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAac                                                                                                                   <  2:56700/149‑1 (MQ=255)
                                  tACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAac                                                                                                                   >  2:398202/1‑149 (MQ=255)
                                     tCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGaa                                                                                                                                                   >  4:468656/1‑115 (MQ=255)
                                     tCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGa                                                                                                                                              >  2:825044/1‑120 (MQ=255)
                                     tCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAAcaca                                                                                                                >  2:1612751/1‑149 (MQ=255)
                                     tCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAAcaca                                                                                                                >  1:1447156/1‑149 (MQ=255)
                                     tCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCACCTTCTGGCGA‑CCCGtttttctttgc                                                                                                                              >  3:1529129/1‑132 (MQ=38)
                                     tCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCACCTTCTGGCGA‑CCCGtttttctttgc                                                                                                                              <  4:1529129/135‑4 (MQ=38)
                                     tCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGACGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAAcaca                                                                                                                >  4:331209/1‑149 (MQ=255)
                                         tGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAAc                                                                                                            <  4:1023444/149‑1 (MQ=255)
                                               tGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACa                                                                                                      >  4:1286144/1‑149 (MQ=255)
                                                 gCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAAc                                                                                                    >  4:1140389/1‑149 (MQ=255)
                                                   cTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGt                                                                                                                                                >  3:1659939/1‑104 (MQ=255)
                                                   cTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCg                                                                                                  >  1:261379/1‑149 (MQ=255)
                                                    ttCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGt                                                                                                 >  2:1527024/1‑149 (MQ=255)
                                                    ttCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGt                                                                                                 >  1:1529215/1‑149 (MQ=255)
                                                     tCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGtt                                                                                                <  3:299108/149‑1 (MQ=255)
                                                         gTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGt                                                                                                                                                <  1:101823/98‑1 (MQ=255)
                                                         gTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGt                                                                                                                                                >  2:101823/1‑98 (MQ=255)
                                                           tCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGc                                                                                          <  2:637144/149‑1 (MQ=255)
                                                           tCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGc                                                                                          >  3:107647/1‑149 (MQ=255)
                                                           tCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGc                                                                                          <  3:1162069/149‑1 (MQ=255)
                                                           tCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGc                                                                                          <  4:1380008/149‑1 (MQ=255)
                                                                      ttCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCg                                                                               <  2:377238/149‑1 (MQ=255)
                                                                         aaaTATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGtt                                                                            >  4:1263983/1‑149 (MQ=255)
                                                                         aaaTATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGtt                                                                            >  2:350891/1‑149 (MQ=255)
                                                                         aaaTATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGtt                                                                            <  1:724458/149‑1 (MQ=255)
                                                                         aaaTATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGtt                                                                            <  1:398202/149‑1 (MQ=255)
                                                                          aaTATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGG‑CCCGTCTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGttt                                                                           <  3:976672/149‑1 (MQ=255)
                                                                           atatAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTc                                                                          <  4:294339/149‑1 (MQ=255)
                                                                             ataAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGc                                                                        <  2:284915/149‑1 (MQ=255)
                                                                             ataAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGc                                                                        <  4:136615/149‑1 (MQ=255)
                                                                              taAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCt                                                                       <  3:1233671/149‑1 (MQ=255)
                                                                              taAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCt                                                                       <  1:36969/149‑1 (MQ=255)
                                                                                aaCACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCtat                                                                     >  3:514816/1‑149 (MQ=255)
                                                                                aaCACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCtat                                                                     <  2:1447156/149‑1 (MQ=255)
                                                                                aaCACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCtat                                                                     >  2:250706/1‑149 (MQ=255)
                                                                                 acacGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTAc                                                                                                       <  3:810059/114‑1 (MQ=255)
                                                                                     gTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAaca                                                                >  3:1004428/1‑149 (MQ=255)
                                                                                         gtgtCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCtt                                                            <  1:1060328/149‑1 (MQ=255)
                                                                                            tCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTc                                                                                                                             <  3:639786/81‑1 (MQ=255)
                                                                                            tCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCa                                                         >  3:274692/1‑149 (MQ=255)
                                                                                              ccGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGc                                                       >  1:993212/1‑149 (MQ=255)
                                                                                               cGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCt                                                      >  4:1115418/1‑149 (MQ=255)
                                                                                                gACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTc                                                     >  2:1720901/1‑149 (MQ=255)
                                                                                                     gCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCg                                                >  3:558439/1‑149 (MQ=255)
                                                                                                       acacaGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAGCCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGaa                                              <  3:346618/149‑1 (MQ=255)
                                                                                                       acacaGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGaa                                              >  2:1682219/1‑149 (MQ=255)
                                                                                                       acacaGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTCTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGaa                                              >  4:326949/1‑149 (MQ=255)
                                                                                                         acaGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGaaaa                                            >  4:1634010/1‑149 (MQ=255)
                                                                                                          caGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAAt                                           >  1:1556566/1‑149 (MQ=255)
                                                                                                            gACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCACCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATcc                                         <  1:1011406/149‑1 (MQ=255)
                                                                                                               ggTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCgg                                      >  2:844384/1‑149 (MQ=255)
                                                                                                                     ccACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGcaacaa                                >  1:1350305/1‑149 (MQ=255)
                                                                                                                     ccACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGcaacaa                                >  2:413103/1‑149 (MQ=255)
                                                                                                                      cACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACaaa                               <  4:717560/149‑1 (MQ=255)
                                                                                                                           ttaaCGCAAAAAAAACCGGTCCGCCA‑AAGGTTACCCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCAc                          <  2:120121/146‑1 (MQ=255)
                                                                                                                           aTTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCAc                          <  2:1294711/149‑1 (MQ=255)
                                                                                                                           aTTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGATTCGCTATAAACACCTTGCAGCTCCGGc                                                 <  3:1447713/126‑1 (MQ=255)
                                                                                                                            ttACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCAcc                         <  1:273519/149‑1 (MQ=255)
                                                                                                                                  aaaGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCcttct                   >  4:1499380/1‑149 (MQ=255)
                                                                                                                                  aaaGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCcttct                   >  4:1432238/1‑149 (MQ=255)
                                                                                                                                  aaaGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCcttct                   >  3:1476865/1‑149 (MQ=255)
                                                                                                                                   aaGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCAc                                                                                            <  4:744849/75‑1 (MQ=255)
                                                                                                                                   aaGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCcttctt                  >  1:1713052/1‑149 (MQ=255)
                                                                                                                                   aaGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCcttctt                  >  1:911062/1‑149 (MQ=255)
                                                                                                                                     aaaaaaCGGGTCACCTTCTGGCGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGc                                                          <  4:1312619/106‑1 (MQ=37)
                                                                                                                                     aaaaaaCGGGTCACCTTCTGGCGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTcc                >  3:997789/2‑149 (MQ=255)
                                                                                                                                      aaaaaCGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCa               <  3:326949/149‑1 (MQ=255)
                                                                                                                                      aaaaaCGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCa               <  1:751870/149‑1 (MQ=255)
                                                                                                                                         aaCGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGaaa                                             >  1:669394/1‑116 (MQ=255)
                                                                                                                                         aaCGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGaaa                                             <  2:669394/116‑1 (MQ=255)
                                                                                                                                            gggTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCaa         <  4:1476865/149‑1 (MQ=255)
                                                                                                                                              gTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAgg       <  4:721108/149‑1 (MQ=255)
                                                                                                                                               tCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACACAACCGAACAACAGTTGCACGCTATCCGCTTAGGTTTCGCTATAAACACCTTTCCGCTCCGGCGACAATCCCGGCGACTAATTCACCCCTTCTTCCAGTCCAAGGa      >  4:225476/1‑149 (MQ=255)
                                                                                                                                                 gCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAGGaaa    >  4:393539/1‑149 (MQ=255)
                                                                                                                                                   cAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAGGAAATa  >  1:126768/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGA‑CCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAGGAAATA  >  NC_000913/3106847‑3107142

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: ‑