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→GGAyfaDtransposase_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

AATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGG  >  NC_000913/2357516‑2357807
                                                                                                                                                      |                                                                                                                                              
aaTGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAggg                                                                                                                                                >  2:767022/1‑150 (MQ=255)
aaTGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAggg                                                                                                                                                >  2:727244/1‑150 (MQ=255)
  tGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCa                                                                                                                                              <  1:622477/150‑1 (MQ=255)
     cagccagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCa                                                                                                                                           >  2:532489/1‑150 (MQ=255)
         cagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAACGGCATAAAAtt                                                                                                                                       >  1:14120/1‑150 (MQ=255)
          agcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtg                                                                                                                                      >  1:641475/1‑150 (MQ=255)
            cgaCAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggc                                                                                                                                    <  1:543459/148‑1 (MQ=255)
            cTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggc                                                                                                                                    >  1:1045888/1‑150 (MQ=255)
             tACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGTATATTGCCGACCGTGCCCCAATCCAAAAGACGGGATAAATGCATATTCCGGACCGCATCCGCAGGAAAGGGCAGCAATTtggtt                                                                                                                                   >  1:1166291/1‑148 (MQ=255)
              agaaaCACTGTTTAATTATCTGGTCCAATGCGGCGATACCTCCCGGTTGACCCGCTTTATTTCGGAGATTGCCTAACGGTACCCACAACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATtggctg                                                                                                                                  <  2:644850/148‑1 (MQ=18)
              aCAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggctg                                                                                                                                  >  1:1230540/1‑150 (MQ=255)
               cAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggctgg                                                                                                                                 <  2:819040/150‑1 (MQ=255)
                aatctcTGTTTAATTATCTGCTGCATTTAGGTGATACCTCGCTTTTCACCCGTTAAATAGCTGAGTCTGACTAACGTTAAACTCAACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTTCGGCAGGAAGGACATCAAATtggctggc                                                                                                                                <  2:1299505/145‑1 (MQ=9)
                aaaCACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggctggc                                                                                                                                <  1:440462/150‑1 (MQ=255)
                aaaCACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggctggc                                                                                                                                <  1:434374/150‑1 (MQ=255)
                 aaCACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCTTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCa                                                                                                                               <  1:284187/150‑1 (MQ=255)
                      tgattAATTAGCTGGGGCTATTCGGCGATACCTCCCGGTCCACCCGTTTTATTGAGGAGTTTGCCTATCGTTCACCACTAAAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAg                                                                                                                          <  1:328462/147‑1 (MQ=18)
                       gTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAgg                                                                                                                         >  1:305015/1‑150 (MQ=255)
                        tttatttaTCTGGTGCATTGCGGCGATACCTTCCGTTGCACCCGCTTTATTGAGGAGAATGCCGATCGTTCACCTCTTAAAAAGGAGAGATTAATGACTATTGCTGAACGGCTTCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGt                                                                                                                        <  2:460521/150‑1 (MQ=18)
                          tAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACTAGCACCACTACAAAAGGAGATATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAt                                                                                                                      <  1:975915/150‑1 (MQ=255)
                           aaTTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATg                                                                                                                     <  2:1324460/150‑1 (MQ=255)
                             ttATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCa                                                                                                                   <  1:910102/150‑1 (MQ=255)
                               aTCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATg                                                                                                                 <  2:147600/150‑1 (MQ=255)
                               aTCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATg                                                                                                                 <  2:348904/150‑1 (MQ=255)
                                tCTGCTGCACTGCGGCGATACCTCCCGTTTCACCCGTGTTATTGAGGAGATTGCCGAACGTTCACCACTACAAACGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGa                                                                                                                <  1:659476/150‑1 (MQ=255)
                                 ctgctgCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGaa                                                                                                               >  2:809577/1‑150 (MQ=255)
                                  tgctgcAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCGGAACGTTCACCACTACAAAAGGAGAGATTAAGGACTATTGCTGAACGGCTACGGCAGGATGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAAc                                                                                                              >  1:212351/1‑150 (MQ=255)
                                      gcAATGCGGCGATACCTCCCCGTTCACCCGATTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGTAGAGATTAATGACTATTGCTCAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGc                                                                                                          <  2:825337/150‑1 (MQ=37)
                                               cgATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAAt                                                                                                 <  1:816659/150‑1 (MQ=255)
                                                  tACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGc                                                                                              <  1:809577/150‑1 (MQ=255)
                                                  tACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATAAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGc                                                                                              <  1:1254541/150‑1 (MQ=255)
                                                   aCCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCt                                                                                             <  2:683650/150‑1 (MQ=255)
                                                     cTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCAAGAACAAGCCCTTCAAATTGCatt                                                                                           >  2:1033900/1‑150 (MQ=255)
                                                      tCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTg                                                                                          <  1:852403/150‑1 (MQ=255)
                                                       cccGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTgc                                                                                         <  1:1224947/150‑1 (MQ=255)
                                                         cGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTgcgc                                                                                       <  1:232699/150‑1 (MQ=255)
                                                          gTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCa                                                                                      >  1:747960/1‑150 (MQ=255)
                                                                cccGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTgg                                                                                <  2:27626/150‑1 (MQ=255)
                                                                cccATTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACATGACATTAAAATTGCTTTGCGCATGCTgt                                                                                >  1:1128839/1‑149 (MQ=255)
                                                                   gTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGAATGAACAAGACATTAAAATTGCTTTGAGCATGCttgtgc                                                                             >  1:606010/1‑145 (MQ=255)
                                                                   gTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGCCATCAAATTGGCAGGAAGGAAGGTATGAATTAACAACCCATTATAATTGATTTGCGCATGCTGGAGc                                                                             >  1:443651/1‑150 (MQ=25)
                                                                       tATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTACTGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAAGTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAggg                                                                         <  1:829772/150‑1 (MQ=255)
                                                                       tATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAggg                                                                         <  1:1320658/150‑1 (MQ=255)
                                                                          tGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGAGCATTGTGGAGCAGGGCtt                                                                      >  1:92036/1‑150 (MQ=255)
                                                                           gaggagATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACATGACATTAAAATTGCATTTAGCATTCATGAGCACGTCttt                                                                     >  2:783/1‑150 (MQ=21)
                                                                               agaTTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCGTTGAAc                                                                 <  2:339602/150‑1 (MQ=255)
                                                                                taTTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTCGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACg                                                                <  2:64795/149‑1 (MQ=255)
                                                                                gaTTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATCCATGAACAAGCCATTTAAAATGCTTTGCGCATGCTGGAGCAGGGCGTTGATCg                                                                >  2:769178/1‑150 (MQ=25)
                                                                                   tGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATtggc                                                                                                                                    <  1:245264/79‑1 (MQ=255)
                                                                                   tGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTga                                                             >  2:635093/1‑150 (MQ=255)
                                                                                   tGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGACATCAAATtggc                                                                                                                                    >  2:245264/1‑79 (MQ=255)
                                                                                    gCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTgag                                                            <  1:698345/150‑1 (MQ=255)
                                                                                    gCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAACGTATCCATGAACAATACAATAAAATTGATTTGAGAATGCTGGAGCAGGGCTTTGATCGTgac                                                            >  2:351906/1‑149 (MQ=25)
                                                                                       gAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGACAtt                                                         >  1:1075638/1‑150 (MQ=255)
                                                                                             tCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAATATTGAACATTACATTAAAATAGCAATGAACATGCTTTAGCAGAGGATTGAACATGACATtgatctt                                                   >  1:50633/1‑145 (MQ=9)
                                                                                              caccacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAgg                                                                                                                         <  1:646450/79‑1 (MQ=255)
                                                                                              caccacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTgg                                                  <  1:1022804/150‑1 (MQ=255)
                                                                                              caccacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGACATCAAATTGGCTGGCAGGAAgg                                                                                                                         >  2:646450/1‑79 (MQ=255)
                                                                                               accacTACCAAAGCCGAGATTAATGACTCTTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGTCAGGAATGTATGCCGGACCAAGCCATTAAAATTGCTTTTCTCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTTGc                                                 <  2:514027/150‑1 (MQ=255)
                                                                                               accacTACAAAATGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGc                                                 <  2:1166504/150‑1 (MQ=255)
                                                                                                ccacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATACAAGAACAATACATTAAACTAGCATTGCTCATTCAGGAGCAGGGCATTGAACCTGACAATGTGCTAGcc                                                >  2:284431/1‑149 (MQ=18)
                                                                                                   ctccAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAGCATGAACAAGCCATTAAAAAAGCATTGCGCATGCTGGAGCAGCTCTTTGAACGTGACATTGTGCTTGCGAcc                                             >  1:992838/4‑149 (MQ=21)
                                                                                                   cTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTTCGAATTCTGGAGCAGTGCTTTGAACGTGACATTGTGCTCGCGAcc                                             >  1:987283/1‑149 (MQ=25)
                                                                                                   cTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGAATGAACAATCCATTAAAATTGCTTTTAGAATTCTGTAGCACTGCATTGAACGTGACAATGTGCTGGCgcac                                             >  1:992835/1‑147 (MQ=18)
                                                                                                      cAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAAcc                                          <  1:4676/150‑1 (MQ=255)
                                                                                                      cAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAGCATGAACAAGCCATTAAAATTGCTATGCGCATGCTGTAGCAGGGCTTTGAACGTGAGATTGAGCTGGCGACAAcc                                          >  1:422100/1‑150 (MQ=255)
                                                                                                      cAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAGCATGAACAAGCCATTAAAAATGCATTGCGCATTGTGGAGCAGGGCTTTGAACGTGGCATTGTGCTGGCGAcccac                                          >  1:648956/1‑146 (MQ=25)
                                                                                                            gagagaTTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTc                                    <  1:387126/150‑1 (MQ=255)
                                                                                                               agaTTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACt                                 <  2:641475/150‑1 (MQ=255)
                                                                                                                gaTTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGTATGAACAATCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGAATGTGCTGGCGCCCACCCACCTCAcag                                >  1:1211377/1‑148 (MQ=25)
                                                                                                                gaTTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTg                                <  2:1166291/150‑1 (MQ=255)
                                                                                                                   tAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATg                             <  2:1230540/150‑1 (MQ=255)
                                                                                                                       taCTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAGGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGAGTGTGCTGGCGACAACCCAACTCACTGATGCTGa                         <  1:198672/149‑1 (MQ=255)
                                                                                                                              gCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTCGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGAAGGTGCAGGGCTTTGAACGTGAGCTTGTGCTGGCGACAAACCACCTAACCGATGCGGATCTTcca                  >  1:291420/1‑149 (MQ=25)
                                                                                                                                  aaCGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGAGCATGCTGGAGCAGGGCTTTGAACGTGACAATGTGCTAGCGACCACCCAACTAACAGAAGCGGATATTccagcca              >  1:1343761/1‑145 (MQ=21)
                                                                                                                                   aCGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTg             <  1:1155511/150‑1 (MQ=255)
                                                                                                                                     ggCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTc           <  2:1073936/150‑1 (MQ=255)
                                                                                                                                     ggCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTc           >  2:642475/1‑150 (MQ=255)
                                                                                                                                         aCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTa       <  2:747960/150‑1 (MQ=255)
                                                                                                                                              aGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAgg  <  1:249476/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                              
AATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGG  >  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: