Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,106,998 A→T 6.2% intergenic (+28/+22) 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,9980AT 6.2% 550.8 / 18.5 160intergenic (+28/+22)nupG/speCnucleoside transporter/ornithine decarboxylase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (78/72);  new base T (3/7);  total (81/79)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 164 total aligned reads displayed.

CTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAGGAAATAACG  >  NC_000913/3106850‑3107145
                                                                                                                                                    |                                                                                                                                                    
atGGCAGTCCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGCCGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGaa                                                                                                                                                      <  2:837039/148‑1 (MQ=255)
   gCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGt                                                                                                                                                   <  3:1539485/149‑1 (MQ=255)
      gACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGAc                                                                                                                                               <  4:184094/149‑1 (MQ=255)
           tATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGtt                                                                                                                                          <  3:417035/149‑1 (MQ=255)
            aTGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCACCTTCTGGCGACCCGttt                                                                                                                                         <  1:175244/149‑1 (MQ=255)
               gTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGtttttt                                                                                                                                      <  1:25416/149‑1 (MQ=255)
               gTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGtttttt                                                                                                                                      >  4:1317633/1‑149 (MQ=255)
               gTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGtttttt                                                                                                                                      >  3:697166/1‑149 (MQ=255)
               gTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGtttttt                                                                                                                                      >  2:702212/1‑149 (MQ=255)
                ttGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGttttttt                                                                                                                                     >  4:867636/1‑149 (MQ=255)
                  gATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTa                                                                                                                                   <  4:864874/149‑1 (MQ=255)
                   aTTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCTCTAATGTCGCAAAGATAAACGGGTCACCTTCTGGCGTCCCGtttttctttg                                                                                                                                  <  2:780343/149‑3 (MQ=255)
                    ttttCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTAtt                                                                                                                                 <  2:977389/149‑1 (MQ=255)
                     tttCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTc                                                                                                                                <  4:697166/149‑1 (MQ=255)
                      ttCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCt                                                                                                                               >  2:175983/1‑149 (MQ=255)
                       tCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCtt                                                                                                                              <  3:966242/149‑1 (MQ=255)
                        cGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTa                                                                                                                             <  3:710489/149‑1 (MQ=255)
                         gCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCTACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTAc                                                                                                                            <  4:1325818/149‑1 (MQ=255)
                         gCTGGTTACCCCGTGGTTCGGGAATTCGCGTTCATGGCGTTGTTCAAATATAATCACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCGCTATTTTCGCAAAGAAAAACGGGTCCCCTCAAGGTGACCCGGTTTTTTTATTCtttc                                                                                                                            <  2:1062237/149‑3 (MQ=255)
                          cTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCACCTTCTGGCGACCCGtttttctt                                                                                                                                    >  4:880953/1‑140 (MQ=38)
                          cTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCACCTTCTGGCGACCCGtttttctt                                                                                                                                    <  3:880953/140‑1 (MQ=38)
                           tGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACtt                                                                                                                          >  3:136615/1‑149 (MQ=255)
                           tGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACtt                                                                                                                          >  4:61218/1‑149 (MQ=255)
                              ttACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCaa                                                                                                                       >  2:724458/1‑149 (MQ=255)
                               tACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAac                                                                                                                      <  1:702212/149‑1 (MQ=255)
                               tACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAac                                                                                                                      >  2:398202/1‑149 (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)
                                  tCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCACCTTCTGGCGACCCGtttttctttgc                                                                                                                                 >  3:1529129/1‑132 (MQ=38)
                                  tCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGACGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAAcaca                                                                                                                   >  4:331209/1‑149 (MQ=255)
                                      tGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAAc                                                                                                               <  4:1023444/149‑1 (MQ=255)
                                           cTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTAc                                                                                                          >  3:973751/1‑149 (MQ=255)
                                            tGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACa                                                                                                         >  4:1286144/1‑149 (MQ=255)
                                                cTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGt                                                                                                                                                   <  4:1659939/104‑1 (MQ=255)
                                                 ttCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGt                                                                                                    >  2:1527024/1‑149 (MQ=255)
                                                 ttCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGt                                                                                                    >  1:1529215/1‑149 (MQ=255)
                                                  tCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGtt                                                                                                   <  3:299108/149‑1 (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)
                                                        tCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGc                                                                                             >  3:107647/1‑149 (MQ=255)
                                                        tCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGc                                                                                             <  3:1162069/149‑1 (MQ=255)
                                                        tCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGc                                                                                             <  2:637144/149‑1 (MQ=255)
                                                        tCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGc                                                                                             <  4:1380008/149‑1 (MQ=255)
                                                            ggCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCTACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATc                                                                                         <  1:1440931/149‑1 (MQ=255)
                                                                      aaaTATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGtt                                                                               <  1:398202/149‑1 (MQ=255)
                                                                      aaaTATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGtt                                                                               <  1:724458/149‑1 (MQ=255)
                                                                      aaaTATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGtt                                                                               >  4:1263983/1‑149 (MQ=255)
                                                                       aaTATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGttt                                                                              >  2:1669571/1‑149 (MQ=255)
                                                                        atatAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTc                                                                             <  4:294339/149‑1 (MQ=255)
                                                                        atatAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTc                                                                             <  2:1353181/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)
                                                                              acacGTTCGTGTCCCGACAGGCATACAGACGGTTAGCCACTAATTATGCAAAGAAACACGGGGCGCCCGAAGGTGCCCCGTTTAATTAATCCTAACAACAACACATAACCGAACAACAGCTTCGCGCCTTCCTCATCGTTTTCGCtata                                                                       >  4:1229926/1‑149 (MQ=255)
                                                                                  gTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAaca                                                                   >  3:1004428/1‑149 (MQ=255)
                                                                                   ttCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAacac                                                                  <  1:980914/149‑1 (MQ=255)
                                                                                      gtgtCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCtt                                                               <  1:1060328/149‑1 (MQ=255)
                                                                                      gtgtCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCtt                                                               <  1:1669571/149‑1 (MQ=255)
                                                                                        gtCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGc                                                             <  4:288028/149‑1 (MQ=255)
                                                                                         tCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTc                                                                                                                                <  3:639786/81‑1 (MQ=255)
                                                                                         tCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCa                                                            >  3:274692/1‑149 (MQ=255)
                                                                                         tCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCa                                                            <  2:1529215/149‑1 (MQ=255)
                                                                                           ccGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGc                                                          >  1:993212/1‑149 (MQ=255)
                                                                                            cGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCt                                                         >  4:1115418/1‑149 (MQ=255)
                                                                                              aCAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTcc                                                       <  4:107647/149‑1 (MQ=255)
                                                                                                  gCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCg                                                   >  2:273519/1‑149 (MQ=255)
                                                                                                  gCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCg                                                   >  3:558439/1‑149 (MQ=255)
                                                                                                    acacaGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAGCCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGaa                                                 <  3:346618/149‑1 (MQ=255)
                                                                                                    acacaGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGaa                                                 >  2:1682219/1‑149 (MQ=255)
                                                                                                    acacaGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTCTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGaa                                                 >  4:326949/1‑149 (MQ=255)
                                                                                                         gACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCACCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATcc                                            <  1:1011406/149‑1 (MQ=255)
                                                                                                          aCGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATccc                                           <  2:393593/149‑1 (MQ=255)
                                                                                                                  ccACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTa                                                                                                                             >  3:325177/1‑59 (MQ=255)
                                                                                                                  ccACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGcaacaa                                   >  2:682922/1‑149 (MQ=255)
                                                                                                                   cACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACaaa                                  <  4:717560/149‑1 (MQ=255)
                                                                                                                        ttaaCGCAAAAAAAACCGGTCCGCCAAAGGTTACCCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCAc                             <  2:120121/146‑1 (MQ=255)
                                                                                                                        aTTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGATTCGCTATAAACACCTTGCAGCTCCGGc                                                    <  3:1447713/126‑1 (MQ=255)
                                                                                                                         ttACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCAcc                            <  1:273519/149‑1 (MQ=255)
                                                                                                                               aaaGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCcttct                      >  4:1432238/1‑149 (MQ=255)
                                                                                                                               aaaGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCcttct                      >  2:751870/1‑149 (MQ=255)
                                                                                                                               aaaGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCcttct                      >  3:1476865/1‑149 (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:1229568/1‑149 (MQ=255)
                                                                                                                                aaGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCcttctt                     >  2:570438/1‑149 (MQ=255)
                                                                                                                                  gAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTcc                   >  1:639843/1‑149 (MQ=255)
                                                                                                                                  aaaaaaCGGGTCACCTTCTGGCGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTcc                   >  3:997789/2‑149 (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)
                                                                                                                                                    |                                                                                                                                                    
CTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACAGACGGTTAGCCACTAATTACGCAAAGAAAAACGGGTCGCCAGAAGGTGACCCGTTTTTTTTATTCTTACTTCAACACATAACCGTACAACCGTTTCACGCCATCCGCATCGGTTTCGCTATAAACACCTTGCAGCTCCGGCGAAAATCCCGGCAACAAATTCACCCCTTCTTCCAGTGCAAGGAAATAACG  >  NC_000913/3106850‑3107145

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: