Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 737,692 T→C 5.3% pseudogene (793/1062 nt) ybfL → putative transposase YbfL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913737,6920TC 5.3% 227.1 / 2.1 75pseudogene (793/1062 nt)ybfLputative transposase YbfL
Reads supporting (aligned to +/- strand):  ref base T (37/34);  new base C (2/2);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

TTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACAT  >  NC_000913/737543‑737840
                                                                                                                                                     |                                                                                                                                                    
tttGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGAtt                                                                                                                                                      <  2:173395/150‑1 (MQ=34)
tttGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGAtt                                                                                                                                                      <  2:1287022/150‑1 (MQ=34)
tttGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGAtt                                                                                                                                                      <  1:547497/150‑1 (MQ=34)
   gCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTaa                                                                                                                                                   >  2:416815/1‑150 (MQ=34)
   gCGATGTCCCTGATGAACTTATTGATTTCACGTGTGAATGGAAAGGACTGAAGAAATTATGCGTGCCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGCTATTAAATCAGTTCTGCTGAtttca                                                                                                                                                   >  2:567457/1‑148 (MQ=255)
       tGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGc                                                                                                                                               >  1:1304114/1‑150 (MQ=34)
            cTGAGGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTATCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGaaa                                                                                                                                          >  2:915464/1‑150 (MQ=255)
             tgatgaACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGaaaa                                                                                                                                         >  1:1406367/1‑150 (MQ=39)
              gatgaACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAg                                                                                                                                        <  2:44031/150‑1 (MQ=40)
              gatgaACTTATTGATTTCACATTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGACTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATCTAACCGCAGAAAAg                                                                                                                                        >  2:1332492/1‑150 (MQ=255)
              gaagaACTTATTGATTTCAAATTTTAATGGAAAGTACTGAACAAATTATTCGTCGCATTCTAATTTCAATCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATCTAACCGCAGAAAAg                                                                                                                                        >  2:1327651/4‑150 (MQ=14)
                 gaACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTc                                                                                                                                     >  1:305480/1‑150 (MQ=255)
                   aCTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGc                                                                                                                                   <  1:516406/150‑1 (MQ=255)
                    cTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGcc                                                                                                                                  <  1:911802/150‑1 (MQ=255)
                    cTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGcc                                                                                                                                  <  2:1288393/150‑1 (MQ=255)
                     ttATTGATTTCACGTTTGAATGGAAAGGACTGAATAAATTATGCGTGGCATTCTCCTTTCGGTCAATAATAGCATAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCca                                                                                                                                 <  1:742328/150‑1 (MQ=255)
                      tATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCcac                                                                                                                                <  1:261696/150‑1 (MQ=255)
                          gATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCa                                                                                                                            >  1:1196236/1‑150 (MQ=255)
                          gATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCa                                                                                                                            >  1:1144837/1‑150 (MQ=255)
                             ttCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCTATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATc                                                                                                                         <  2:1304114/150‑1 (MQ=255)
                              tCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATcc                                                                                                                        >  2:1036773/1‑150 (MQ=255)
                               cACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTTCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCg                                                                                                                       <  2:614141/150‑1 (MQ=255)
                                  gTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGaaa                                                                                                                    <  1:416815/150‑1 (MQ=255)
                                       aaTGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACt                                                                                                               <  2:1064137/150‑1 (MQ=255)
                                         tGGAAAGGACTGAAGAAATTATCCGTGGCAGTCTCCTTTCGGACAATAATAGCAGAACAAAAGAAAGAGCCATAAATGACGGTAAGAAATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCTCCACAGAAATCCTAAACCACTgg                                                                                                             >  2:480843/1‑150 (MQ=14)
                                          ggAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGc                                                                                                            >  1:291482/1‑150 (MQ=255)
                                                gACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTgg                                                                                                      >  2:1240132/1‑150 (MQ=255)
                                                gACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTgg                                                                                                      >  1:1028738/1‑150 (MQ=255)
                                                  cTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGgag                                                                                                    >  2:1050501/1‑150 (MQ=255)
                                                    gaagaaATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACg                                                                                                         <  2:946679/143‑1 (MQ=40)
                                                    gaagaaATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACg                                                                                                         >  1:946679/1‑143 (MQ=40)
                                                    gaagaaATTATGCGTGGCAGTCTCCTTTCGGTCAATAAAAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACATCAGAAAAGTTAGCAACAGCAAACCGAAACCACTGGAACGTGGAGaa                                                                                                  >  1:1392864/1‑150 (MQ=14)
                                                      agaaATTATGAGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAAACACTGGCACGTGGAGAATa                                                                                                >  2:515226/1‑150 (MQ=255)
                                                        aaaTTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAACGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAg                                                                                              <  1:945595/150‑1 (MQ=255)
                                                            tATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGc                                                                                          <  1:1230365/150‑1 (MQ=255)
                                                              tGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCAc                                                                                        >  1:238276/1‑150 (MQ=255)
                                                                       gTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTg                                                                               <  2:238276/150‑1 (MQ=255)
                                                                       gTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTg                                                                               >  1:1347071/1‑150 (MQ=255)
                                                                           ccTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACg                                                                           >  1:141692/1‑150 (MQ=255)
                                                                               tCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt                                                                       <  2:1144837/150‑1 (MQ=255)
                                                                               tCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt                                                                       <  1:1050501/150‑1 (MQ=255)
                                                                                  gTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaat                                                                    >  1:26095/1‑150 (MQ=255)
                                                                                     aataatAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaa                                                                 <  2:256546/150‑1 (MQ=40)
                                                                                       taataGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaatg                                                               <  2:1406367/150‑1 (MQ=39)
                                                                                        aataGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAATCTGCACTGGCGTCTGGACGTGGTAatgaatga                                                              <  1:515226/150‑1 (MQ=255)
                                                                                         ataGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAAtgaatgaa                                                             <  2:1266263/150‑1 (MQ=39)
                                                                                          taGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAg                                                            <  2:305480/150‑1 (MQ=39)
                                                                                           aGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAga                                                           >  1:307443/1‑150 (MQ=39)
                                                                                           aGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAga                                                           >  1:43896/1‑150 (MQ=39)
                                                                                           aGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAga                                                           <  2:9819/150‑1 (MQ=39)
                                                                                               gAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAatgaatga                                                              <  1:798852/143‑1 (MQ=39)
                                                                                               gAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAatgaatga                                                              >  2:798852/1‑143 (MQ=39)
                                                                                                 aCAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTg                                                     >  1:623453/1‑150 (MQ=39)
                                                                                                 aCAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTg                                                     <  2:46309/150‑1 (MQ=39)
                                                                                                  cAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGc                                                    >  1:1079534/1‑150 (MQ=39)
                                                                                                   aaaaGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCa                                                   >  1:1133425/1‑150 (MQ=39)
                                                                                                    aaagaaagAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCaa                                                  >  1:1388695/1‑150 (MQ=39)
                                                                                                      agaaagaGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCaaaa                                                >  1:664574/1‑150 (MQ=39)
                                                                                                           gagCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATCTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATaaga                                           <  1:1332492/150‑1 (MQ=255)
                                                                                                           gagCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATCTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATaaga                                           <  1:1327651/150‑1 (MQ=255)
                                                                                                             gCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATaagaag                                         <  2:1196236/150‑1 (MQ=38)
                                                                                                             gCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATaagaag                                         >  1:893729/1‑150 (MQ=38)
                                                                                                              ccAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAagaaga                                        >  1:477056/1‑150 (MQ=37)
                                                                                                                       aCGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAAcgcc                               <  2:1028738/150‑1 (MQ=38)
                                                                                                                          gTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCa                            <  2:1014478/150‑1 (MQ=39)
                                                                                                                             aGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGaa                         <  2:1347071/150‑1 (MQ=39)
                                                                                                                              tatatTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATTAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAat                        <  2:367662/149‑1 (MQ=255)
                                                                                                                                 attatATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAattat                     <  1:34257/150‑1 (MQ=39)
                                                                                                                                 attatATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAattat                     >  2:107641/1‑150 (MQ=39)
                                                                                                                                   tataTCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTAttt                   <  1:1095447/150‑1 (MQ=39)
                                                                                                                                    atatCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTAtttt                  >  2:1034448/1‑150 (MQ=39)
                                                                                                                                      atCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCa                <  1:211911/150‑1 (MQ=39)
                                                                                                                                      atCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCa                >  1:1064716/1‑150 (MQ=39)
                                                                                                                                       tCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAg               >  2:1423217/1‑150 (MQ=39)
                                                                                                                                        cAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAgg              >  2:28500/1‑150 (MQ=39)
                                                                                                                                        cAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAgg              <  2:141692/150‑1 (MQ=39)
                                                                                                                                         aGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAggg             >  2:575055/1‑150 (MQ=39)
                                                                                                                                                 tgATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGca     <  1:555844/150‑1 (MQ=39)
                                                                                                                                                   aTTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGcaca   >  1:313445/1‑150 (MQ=39)
                                                                                                                                                    tttAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACAt  >  1:602254/1‑150 (MQ=39)
                                                                                                                                                     |                                                                                                                                                    
TTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGACTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACAT  >  NC_000913/737543‑737840

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