Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,413,029 G→C 6.8% S69R (AGC→AGG) ‡ intR ← Rac prophage; integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,413,0290GC 6.8% 248.7 / 8.2 73S69R (AGC→AGG) ‡intRRac prophage; integrase
Reads supporting (aligned to +/- strand):  ref base G (39/29);  new base C (3/2);  total (42/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

ACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAAC  >  NC_000913/1412882‑1413176
                                                                                                                                                   |                                                                                                                                                   
acaAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGc                                                                                                                                                    >  3:597777/1‑149 (MQ=255)
   aaTTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTTCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTgg                                                                                                                                                 >  1:408515/1‑149 (MQ=255)
    attattattCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTgac                                                                                                                                                >  3:579889/1‑147 (MQ=255)
      tattattCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACt                                                                                                                                              <  2:264733/149‑1 (MQ=255)
       attattCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGTACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTg                                                                                                                                             >  3:451690/1‑149 (MQ=255)
        ttattCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTTCCCCAGTCCAGCACTgg                                                                                                                                            <  4:405557/149‑1 (MQ=255)
             cTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCT‑GACTGGGGaaaa                                                                                                                                      <  4:579889/149‑1 (MQ=255)
              tGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTggg                                                                                                                                           >  3:698764/1‑144 (MQ=255)
              tGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTggg                                                                                                                                           <  4:698764/144‑1 (MQ=255)
              tGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGaaaa                                                                                                                                      >  1:349157/1‑149 (MQ=255)
              tGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGaaaa                                                                                                                                      >  4:739006/1‑149 (MQ=255)
                  tgtTTTTTATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTc                                                                                                                                  <  1:377468/149‑1 (MQ=255)
                   gtTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCt                                                                                                                                 <  2:344515/149‑1 (MQ=255)
                    tttttAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTTCCCCAGTCCAGCACTGGGGcaaaaatggt                                                                                                                                <  3:506903/149‑7 (MQ=255)
                    tttttAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCtt                                                                                                                                >  4:733168/1‑149 (MQ=255)
                        tAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAAt                                                                                                                            <  4:653104/149‑1 (MQ=255)
                         aaTTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATa                                                                                                                           <  3:739006/149‑1 (MQ=255)
                                aTTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAAt                                                                                                                    >  3:753589/1‑149 (MQ=255)
                                 tttCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAAtt                                                                                                                   >  4:620936/1‑149 (MQ=255)
                                   tCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCa                                                                                                                 <  4:753589/149‑1 (MQ=255)
                                        ttcttcATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGa                                                                                                            <  3:355050/149‑1 (MQ=255)
                                         tcttcATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAg                                                                                                           <  3:733168/149‑1 (MQ=255)
                                                 gtgtTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTc                                                                                                   >  3:4397/1‑149 (MQ=255)
                                                    tttgttGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATg                                                                                                >  2:271879/1‑149 (MQ=255)
                                                      tgttgtTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGaa                                                                                              <  4:279934/149‑1 (MQ=255)
                                                       gttgttAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAAc                                                                                             >  4:675173/1‑149 (MQ=255)
                                                           ttAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATc                                                                                         >  3:278746/1‑149 (MQ=255)
                                                           ttAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATc                                                                                         >  1:146458/1‑149 (MQ=255)
                                                            tAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCa                                                                                        <  2:605552/149‑1 (MQ=255)
                                                               ctctGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAgcg                                                                                     <  3:620936/149‑1 (MQ=255)
                                                                    tctcTTTAATTTTTGTACAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCtt                                                                                >  2:148055/1‑149 (MQ=255)
                                                                       ctTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTaaa                                                                             <  1:32160/149‑1 (MQ=255)
                                                                        tttAATTTTTGTCTAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAAt                                                                            <  1:252072/149‑1 (MQ=255)
                                                                        tttAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAAt                                                                            >  3:440349/1‑149 (MQ=255)
                                                                          tAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAAttt                                                                          <  3:654582/149‑1 (MQ=255)
                                                                            aTTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTcc                                                                        <  1:678158/149‑1 (MQ=255)
                                                                             tttttGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCg                                                                       >  2:225422/1‑149 (MQ=255)
                                                                                  gTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCtt                                                                  <  4:4397/149‑1 (MQ=255)
                                                                                     cacatatCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCcttctt                                                               <  3:504630/146‑1 (MQ=255)
                                                                                     cAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGaa                                                                                              <  4:235274/118‑1 (MQ=255)
                                                                                     cAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGaa                                                                                              >  3:235274/1‑118 (MQ=255)
                                                                                     cAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCttttt                                                               >  4:289093/1‑149 (MQ=255)
                                                                                     cAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCttttt                                                               >  1:185735/1‑149 (MQ=255)
                                                                                       gATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTaa                                                             >  2:60286/1‑149 (MQ=255)
                                                                                        atatCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAat                                                            <  3:111455/149‑1 (MQ=255)
                                                                                         tatCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAata                                                           >  1:532860/1‑149 (MQ=255)
                                                                                            cacaAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTa                                                        <  4:597777/149‑1 (MQ=255)
                                                                                               aaaGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTg                                                     >  4:489622/1‑149 (MQ=255)
                                                                                               aaaGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGACTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTg                                                     >  1:274961/1‑149 (MQ=255)
                                                                                                  gCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGttt                                                  >  1:365850/1‑149 (MQ=255)
                                                                                                    tctcCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGt                                                <  3:394914/149‑1 (MQ=255)
                                                                                                          aaCGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAAtctc                                          >  3:130341/1‑149 (MQ=255)
                                                                                                             gTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCa                                       <  4:497409/149‑1 (MQ=255)
                                                                                                              tttttATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTa                                               >  4:83733/1‑140 (MQ=255)
                                                                                                              tttttATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTa                                               <  3:83733/140‑1 (MQ=255)
                                                                                                                     gACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTa                               >  1:511641/1‑149 (MQ=255)
                                                                                                                       ctctCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAAt                             <  3:324893/149‑1 (MQ=255)
                                                                                                                       ctctCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAAt                             <  4:634717/149‑1 (MQ=255)
                                                                                                                            gTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTc                        >  3:102733/1‑149 (MQ=255)
                                                                                                                             ttGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATTTTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTcc                       <  4:550740/149‑1 (MQ=255)
                                                                                                                                tCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGa                    >  2:553457/1‑149 (MQ=255)
                                                                                                                                  aagacgttttCCCCAGTCCAGCACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCTCAGCCTTTTAATGTTTCCCGACa                  >  4:255709/7‑149 (MQ=255)
                                                                                                                                  aCCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACa                  >  2:625000/1‑149 (MQ=255)
                                                                                                                                  aCCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACa                  >  2:694723/1‑149 (MQ=255)
                                                                                                                                  aCCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACa                  >  4:464728/1‑149 (MQ=255)
                                                                                                                                  aCCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACa                  >  2:305570/1‑149 (MQ=255)
                                                                                                                                  aCCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACa                  >  3:195131/1‑149 (MQ=255)
                                                                                                                                   ccATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACAt                 >  3:7322/1‑149 (MQ=255)
                                                                                                                                          ttccccAGTCCAGCACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTc          >  2:257682/3‑149 (MQ=255)
                                                                                                                                            ccccAGTCCAGCACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCt        >  4:330277/1‑149 (MQ=255)
                                                                                                                                             cccAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTc       >  1:63080/1‑149 (MQ=255)
                                                                                                                                              ccAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCg      >  3:399444/1‑149 (MQ=255)
                                                                                                                                                  tGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAAc  <  2:511641/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACAAATTATTATTCTGAGTGTTTTTAATTGTGATTTCGTTTTCTTCATTGTGTTTGTTGTTAACTCTGTCTCTTTAATTTTTGTCCAGATATCACAAAGCTCTCCGAACGTTTTTATGACTCTCGTTGTCACCATTTTTGCCCCAGTGCTGGACTGGGGAAAACGTCTTAAATACTCAAATTCACCGGAGTTTATTTCATGAACTATCAGCGCTCTTAAATTTCCGGCCTTTTTAATATTACTGTTTGTAATCTCCCAGCCTTTTAATGTTTCCCGACATCGTTTTCCTCGAAAC  >  NC_000913/1412882‑1413176

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑