Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,357,665 G→A 10.4% G254G (GGG→GGA)  yfaD → transposase_31 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,357,6650GA10.4% 202.4 / 10.2 78G254G (GGG→GGA) yfaDtransposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base G (30/39);  new base A (2/6);  total (33/45)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGG  >  NC_000913/2357516‑2357807
                                                                                                                                                      |                                                                                                                                              
aaTGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAggg                                                                                                                                                >  2:767022/1‑150 (MQ=255)
aaTGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAggg                                                                                                                                                >  2:727244/1‑150 (MQ=255)
  tGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCa                                                                                                                                              <  1:622477/150‑1 (MQ=255)
     cagccagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCa                                                                                                                                           >  2:532489/1‑150 (MQ=255)
         cagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGG‑AACGGCATAAAAtt                                                                                                                                       >  1:14120/1‑150 (MQ=255)
          agcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATtg                                                                                                                                      >  1:641475/1‑150 (MQ=255)
            cgaCAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATtggc                                                                                                                                    <  1:543459/148‑1 (MQ=255)
            cTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATtggc                                                                                                                                    >  1:1045888/1‑150 (MQ=255)
             tACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGTATATTGCCGACCGTGCCCCAATCCAAAAGACGGGATAAATGCATATT‑CCGGACCGCATCCGCAGGAAAGGGCAGCAATTtggtt                                                                                                                                   >  1:1166291/1‑148 (MQ=255)
              agaaaCACTGTTTAATTATCTGGTCCAATGCGGCGATACCTCCCGGTTGACCCGCTTTATTTCGGAGATTGCCTAACGGTACCCACAACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGG‑AAGGACATCAAATtggctg                                                                                                                                  <  2:644850/148‑1 (MQ=18)
              aCAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATtggctg                                                                                                                                  >  1:1230540/1‑150 (MQ=255)
               cAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATtggctgg                                                                                                                                 <  2:819040/150‑1 (MQ=255)
                aatctcTGTTTAATTATCTGCTGCATTTAGGTGATACCTCGCTTTTCACCCGTTAAATAGCTGAGTCTGACTAACGTTAAACTCAACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTTCGGCAGG‑AAGGACATCAAATtggctggc                                                                                                                                <  2:1299505/145‑1 (MQ=9)
                aaaCACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATtggctggc                                                                                                                                <  1:440462/150‑1 (MQ=255)
                aaaCACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATtggctggc                                                                                                                                <  1:434374/150‑1 (MQ=255)
                 aaCACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCTTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCa                                                                                                                               <  1:284187/150‑1 (MQ=255)
                      tgattAATTAGCTGGGGCTATTCGGCGATACCTCCCGGTCCACCCGTTTTATTGAGGAGTTTGCCTATCGTTCACCACTAAAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGG‑AAGGACATCAAATTGGCTGGCAGGAAg                                                                                                                          <  1:328462/147‑1 (MQ=18)
                       gTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAgg                                                                                                                         >  1:305015/1‑150 (MQ=255)
                        tttatttaTCTGGTGCATTGCGGCGATACCTTCCGTTGCACCCGCTTTATTGAGGAGAATGCCGATCGTTCACCTCTTAAAAAGGAGAGATTAATGACTATTGCTGAACGGCTTCGGCAGG‑AAGGACATCAAATTGGCTGGCAGGAAGGt                                                                                                                        <  2:460521/150‑1 (MQ=18)
                          tAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACTAGCACCACTACAAAAGGAGATATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTAt                                                                                                                      <  1:975915/150‑1 (MQ=255)
                           aaTTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATg                                                                                                                     <  2:1324460/150‑1 (MQ=255)
                             ttATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCa                                                                                                                   <  1:910102/150‑1 (MQ=255)
                               aTCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATg                                                                                                                 <  2:147600/150‑1 (MQ=255)
                               aTCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATg                                                                                                                 <  2:348904/150‑1 (MQ=255)
                                tCTGCTGCACTGCGGCGATACCTCCCGTTTCACCCGTGTTATTGAGGAGATTGCCGAACGTTCACCACTACAAACGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGa                                                                                                                <  1:659476/150‑1 (MQ=255)
                                 ctgctgCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGaa                                                                                                               >  2:809577/1‑150 (MQ=255)
                                  tgctgcAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCGGAACGTTCACCACTACAAAAGGAGAGATTAAGGACTATTGCTGAACGGCTACGGCAGG‑ATGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAAc                                                                                                              >  1:212351/1‑150 (MQ=255)
                                      gcAATGCGGCGATACCTCCCCGTTCACCCGATTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGTAGAGATTAATGACTATTGCTCAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGc                                                                                                          <  2:825337/150‑1 (MQ=37)
                                               cgATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAAt                                                                                                 <  1:816659/150‑1 (MQ=255)
                                                  tACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGc                                                                                              <  1:809577/150‑1 (MQ=255)
                                                  tACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATAAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGc                                                                                              <  1:1254541/150‑1 (MQ=255)
                                                   aCCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCt                                                                                             <  2:683650/150‑1 (MQ=255)
                                                     cTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCAAGAACAAGCCCTTCAAATTGCatt                                                                                           >  2:1033900/1‑150 (MQ=255)
                                                      tCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTg                                                                                          <  1:852403/150‑1 (MQ=255)
                                                       cccGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTgc                                                                                         <  1:1224947/150‑1 (MQ=255)
                                                         cGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTgcgc                                                                                       <  1:232699/150‑1 (MQ=255)
                                                          gTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCa                                                                                      >  1:747960/1‑150 (MQ=255)
                                                                cccGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTgg                                                                                <  2:27626/150‑1 (MQ=255)
                                                                cccATTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACATGACATTAAAATTGCTTTGCGCATGCTgt                                                                                >  1:1128839/1‑149 (MQ=255)
                                                                   gTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGAATGAACAAGACATTAAAATTGCTTTGAGCATGCttgtgc                                                                             >  1:606010/1‑145 (MQ=255)
                                                                   gTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGCCATCAAATTGGCAGGAAGGAAGGTATGAATTAACAACCCATTATAATTGATTTGCGCATGCTGGAGc                                                                             >  1:443651/1‑150 (MQ=25)
                                                                       tATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTACTGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAAGTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAggg                                                                         <  1:829772/150‑1 (MQ=255)
                                                                       tATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAggg                                                                         <  1:1320658/150‑1 (MQ=255)
                                                                          tGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGAGCATTGTGGAGCAGGGCtt                                                                      >  1:92036/1‑150 (MQ=255)
                                                                           gaggagATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACATGACATTAAAATTGCATTTAGCATTCATGAGCACGTCttt                                                                     >  2:783/1‑150 (MQ=21)
                                                                               agaTTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCGTTGAAc                                                                 <  2:339602/150‑1 (MQ=255)
                                                                                taTTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTCGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACg                                                                <  2:64795/149‑1 (MQ=255)
                                                                                gaTTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATCCATGAACAAGCCATTTAAAATGCTTTGCGCATGCTGGAGCAGGGCGTTGATCg                                                                >  2:769178/1‑150 (MQ=25)
                                                                                   tGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGG‑AAGGACATCAAATtggc                                                                                                                                    <  1:245264/79‑1 (MQ=255)
                                                                                   tGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTga                                                             >  2:635093/1‑150 (MQ=255)
                                                                                   tGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGACATCAAATtggc                                                                                                                                    >  2:245264/1‑79 (MQ=255)
                                                                                    gCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTgag                                                            <  1:698345/150‑1 (MQ=255)
                                                                                    gCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAACGTATCCATGAACAATACAATAAAATTGATTTGAGAATGCTGGAGCAGGGCTTTGATCGTgac                                                            >  2:351906/1‑149 (MQ=25)
                                                                                       gAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGACAtt                                                         >  1:1075638/1‑150 (MQ=255)
                                                                                             tCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTAATATTGAACATTACATTAAAATAGCAATGAACATGCTTTAGCAGAGGATTGAACATGACATtgatctt                                                   >  1:50633/1‑145 (MQ=9)
                                                                                              caccacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGG‑AAGGACATCAAATTGGCTGGCAGGAAgg                                                                                                                         <  1:646450/79‑1 (MQ=255)
                                                                                              caccacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTgg                                                  <  1:1022804/150‑1 (MQ=255)
                                                                                              caccacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGACATCAAATTGGCTGGCAGGAAgg                                                                                                                         >  2:646450/1‑79 (MQ=255)
                                                                                               accacTACCAAAGCCGAGATTAATGACTCTTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGTCAGGAATGTATGCCGGACCAAGCCATTAAAATTGCTTTTCTCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTTGc                                                 <  2:514027/150‑1 (MQ=255)
                                                                                               accacTACAAAATGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGc                                                 <  2:1166504/150‑1 (MQ=255)
                                                                                                ccacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATACAAGAACAATACATTAAACTAGCATTGCTCATTCAGGAGCAGGGCATTGAACCTGACAATGTGCTAGcc                                                >  2:284431/1‑149 (MQ=18)
                                                                                                   ctccAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTAAGCATGAACAAGCCATTAAAAAAGCATTGCGCATGCTGGAGCAGCTCTTTGAACGTGACATTGTGCTTGCGAcc                                             >  1:992838/4‑149 (MQ=21)
                                                                                                   cTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTTCGAATTCTGGAGCAGTGCTTTGAACGTGACATTGTGCTCGCGAcc                                             >  1:987283/1‑149 (MQ=25)
                                                                                                   cTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGAATGAACAATCCATTAAAATTGCTTTTAGAATTCTGTAGCACTGCATTGAACGTGACAATGTGCTGGCgcac                                             >  1:992835/1‑147 (MQ=18)
                                                                                                      cAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAAcc                                          <  1:4676/150‑1 (MQ=255)
                                                                                                      cAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTAAGCATGAACAAGCCATTAAAATTGCTATGCGCATGCTGTAGCAGGGCTTTGAACGTGAGATTGAGCTGGCGACAAcc                                          >  1:422100/1‑150 (MQ=255)
                                                                                                      cAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTAAGCATGAACAAGCCATTAAAAATGCATTGCGCATTGTGGAGCAGGGCTTTGAACGTGGCATTGTGCTGGCGAcccac                                          >  1:648956/1‑146 (MQ=25)
                                                                                                            gagagaTTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTc                                    <  1:387126/150‑1 (MQ=255)
                                                                                                               agaTTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACt                                 <  2:641475/150‑1 (MQ=255)
                                                                                                                gaTTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGTATGAACAATCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGAATGTGCTGGCGCCCACCCACCTCAcag                                >  1:1211377/1‑148 (MQ=25)
                                                                                                                gaTTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTg                                <  2:1166291/150‑1 (MQ=255)
                                                                                                                   tAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATg                             <  2:1230540/150‑1 (MQ=255)
                                                                                                                       taCTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTAGGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGAGTGTGCTGGCGACAACCCAACTCACTGATGCTGa                         <  1:198672/149‑1 (MQ=255)
                                                                                                                              gCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTCGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGAAGGTGCAGGGCTTTGAACGTGAGCTTGTGCTGGCGACAAACCACCTAACCGATGCGGATCTTcca                  >  1:291420/1‑149 (MQ=25)
                                                                                                                                  aaCGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGAGCATGCTGGAGCAGGGCTTTGAACGTGACAATGTGCTAGCGACCACCCAACTAACAGAAGCGGATATTccagcca              >  1:1343761/1‑145 (MQ=21)
                                                                                                                                   aCGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTg             <  1:1155511/150‑1 (MQ=255)
                                                                                                                                     ggCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTc           <  2:1073936/150‑1 (MQ=255)
                                                                                                                                     ggCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTc           >  2:642475/1‑150 (MQ=255)
                                                                                                                                         aCGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTa       <  2:747960/150‑1 (MQ=255)
                                                                                                                                              aGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAgg  <  1:249476/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                              
AATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGG‑AAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGG  >  NC_000913/2357516‑2357807

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑