Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,777,573 A→G 6.1% T41A (ACG→GCG) ‡ ypjF → CP4‑57 prophage; toxin of the YpjF‑YfjZ toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,777,5730AG 6.1% 247.2 / 5.3 82T41A (ACG→GCG) ‡ypjFCP4‑57 prophage; toxin of the YpjF‑YfjZ toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base A (39/38);  new base G (2/3);  total (41/41)
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.26e-01

CAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTT  >  NC_000913/2777424‑2777720
                                                                                                                                                     |                                                                                                                                                   
cAACGCTAACTAATTTCACGAGAGCACGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGaca                                                                                                                                                     >  2:1451888/1‑150 (MQ=255)
cAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGaca                                                                                                                                                     >  2:1743055/1‑150 (MQ=255)
cAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGaca                                                                                                                                                     >  1:667815/1‑150 (MQ=255)
  aCGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACg                                                                                                                                                   <  1:1356468/150‑1 (MQ=255)
        actaaTTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTc                                                                                                                                             >  2:722759/1‑150 (MQ=255)
           aaTTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGACTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGGCGTGTCGTTCAGt                                                                                                                                          <  1:1302355/150‑1 (MQ=255)
            aTTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGACAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGGCGTGTCGTTCAGtg                                                                                                                                         <  1:1583327/150‑1 (MQ=255)
              ttCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGtgat                                                                                                                                       >  1:798892/1‑150 (MQ=255)
              ttCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGACAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGtgat                                                                                                                                       <  1:390164/150‑1 (MQ=255)
               tCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGtgatg                                                                                                                                      >  1:115041/1‑150 (MQ=255)
                cACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGtgatga                                                                                                                                     >  2:734391/1‑150 (MQ=255)
                      agcaagcaTGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGt                                                                                                                               <  1:427125/150‑1 (MQ=255)
                       gcaagcaTGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGtt                                                                                                                              <  1:787513/150‑1 (MQ=255)
                       gcaagcaTGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGtt                                                                                                                              <  1:660032/150‑1 (MQ=255)
                        caagcaTGAACACTCTACCTGCTACAATTTTGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGGCGTGTCGTTCAGtgtca                                                                                                                                      <  1:38674/141‑4 (MQ=255)
                        caagcaTGAACACTCTACCTGCTACAATTTTGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGGCGTGTCGTTCAGtgtca                                                                                                                                      >  2:38674/1‑138 (MQ=255)
                        caagcaTGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGtta                                                                                                                             >  2:1479940/1‑150 (MQ=255)
                        caagcaTGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGtta                                                                                                                             >  1:531529/1‑150 (MQ=255)
                         aagcaTGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGttat                                                                                                                            <  1:1753061/150‑1 (MQ=255)
                                            gCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCt                                                                                                         >  1:1917469/1‑150 (MQ=255)
                                              tACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTgg                                                                                                       <  2:502159/150‑1 (MQ=255)
                                               aCAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCATGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGt                                                                                                      >  2:1100131/1‑150 (MQ=255)
                                                      cGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCActc                                                                                               <  2:531529/150‑1 (MQ=255)
                                                      cGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCActc                                                                                               <  1:672980/150‑1 (MQ=255)
                                                        cAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTg                                                                                             >  1:1221724/1‑150 (MQ=255)
                                                        cAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTg                                                                                             >  2:199413/1‑150 (MQ=255)
                                                         cggcggcGATGCCCTGCCTGTCGCCTGGGGCTGTCTGGCAAATGCTACTGACCCGCCGGCTGGAACAACACTATGGCGTGACACTTACCGACACGCCGTGCAGTGATGAAACTGTTATTAAGTAACATATCGATGCTGGTATCACTGTgg                                                                                            <  2:1718605/149‑1 (MQ=255)
                                                          ggcggcGAAGCCCTGCCTGTCGCCAGTGGCTGTTTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCt                                                                                                         >  2:52983/1‑136 (MQ=255)
                                                          ggcggcGAAGCCCTGCCTGTCGCCAGTGGCTGTTTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCt                                                                                                         <  1:52983/136‑1 (MQ=255)
                                                             ggcgAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAATGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGa                                                                                        >  1:1625647/1‑150 (MQ=255)
                                                               cgAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTTAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATg                                                                                      <  1:1490638/150‑1 (MQ=255)
                                                               cgAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTTAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATg                                                                                      <  1:1483882/150‑1 (MQ=255)
                                                                gAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGc                                                                                     <  2:126953/150‑1 (MQ=255)
                                                                 aaGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCa                                                                                    <  2:667815/150‑1 (MQ=255)
                                                                      ctgcctgTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCaa                                                                               <  1:834924/150‑1 (MQ=255)
                                                                       tgcctgTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAAt                                                                              >  2:1549480/1‑150 (MQ=255)
                                                                        gcctgTCGCCAGTGGCTGTCTGGCAAATGGTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAAtt                                                                             >  2:478478/1‑150 (MQ=255)
                                                                            gTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCt                                                                         >  1:361353/1‑150 (MQ=255)
                                                                            gTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCt                                                                         <  1:734391/150‑1 (MQ=255)
                                                                               gCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCtggt                                                                      >  1:230764/1‑150 (MQ=255)
                                                                               gCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCtggt                                                                      >  2:1052870/1‑150 (MQ=255)
                                                                               gCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCtggt                                                                      >  1:1083027/1‑150 (MQ=255)
                                                                                      gCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAg                                                               <  1:580265/150‑1 (MQ=255)
                                                                                      gCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAg                                                               <  1:1743055/150‑1 (MQ=255)
                                                                                        tGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGTTGCAGTCAATTTTCTGGTGGAAAAGTa                                                             >  1:92378/1‑150 (MQ=255)
                                                                                        tGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTa                                                             <  1:1881389/150‑1 (MQ=255)
                                                                                        tGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTa                                                             >  2:1280091/1‑150 (MQ=255)
                                                                                         gTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTAc                                                            <  2:1300484/150‑1 (MQ=255)
                                                                                               cAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTg                                                      >  1:922855/1‑150 (MQ=255)
                                                                                                 aaTGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGt                                                    >  1:1730274/1‑150 (MQ=255)
                                                                                                 aaTGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGt                                                    <  2:1103389/150‑1 (MQ=255)
                                                                                                        cTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATc                                             <  1:1394160/150‑1 (MQ=255)
                                                                                                        cTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATc                                             >  1:1410262/1‑150 (MQ=255)
                                                                                                            cacGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATc                                         >  2:848644/1‑150 (MQ=255)
                                                                                                            cacGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATc                                         >  1:1440991/1‑150 (MQ=255)
                                                                                                                ccTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACag                                     >  1:1310654/1‑150 (MQ=255)
                                                                                                                   gctgGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAgg                                  >  2:463010/1‑150 (MQ=255)
                                                                                                                      ggAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGAtt                               <  1:249200/150‑1 (MQ=255)
                                                                                                                        aacaacACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGAtttt                             <  2:1772945/150‑1 (MQ=255)
                                                                                                                         acaacaCTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTc                            >  2:475041/1‑150 (MQ=255)
                                                                                                                          caacaCTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCg                           >  1:184593/1‑150 (MQ=255)
                                                                                                                           aacaCTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGt                          <  1:209014/150‑1 (MQ=255)
                                                                                                                           aacaCTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGt                          >  2:677207/1‑150 (MQ=255)
                                                                                                                            acacTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTg                         <  2:1083027/150‑1 (MQ=255)
                                                                                                                             cacTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTgg                        <  2:1917469/150‑1 (MQ=255)
                                                                                                                             cacTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTgg                        <  1:722759/150‑1 (MQ=255)
                                                                                                                               cTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGca                      >  2:1662280/1‑150 (MQ=255)
                                                                                                                               cTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGca                      <  1:463010/150‑1 (MQ=255)
                                                                                                                                tATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGcaa                     >  1:316835/1‑150 (MQ=255)
                                                                                                                                 aTGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGcaac                    >  1:1233271/1‑150 (MQ=255)
                                                                                                                                   ggCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCTATCACAGAGGATTTTCGTGGcaacaa                  >  1:1821610/1‑150 (MQ=255)
                                                                                                                                   ggCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGcaacaa                  >  2:1576385/1‑150 (MQ=255)
                                                                                                                                     ccTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAacaaca                <  2:798892/150‑1 (MQ=255)
                                                                                                                                      cTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAg               >  1:1811045/1‑150 (MQ=255)
                                                                                                                                        gACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGtc             <  1:1830823/150‑1 (MQ=255)
                                                                                                                                          cacTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGtctc           <  1:1576385/150‑1 (MQ=255)
                                                                                                                                          cacTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGtctc           <  2:922855/150‑1 (MQ=255)
                                                                                                                                           acTGAACGGCGTGTCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGAAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTcc          >  2:430738/1‑150 (MQ=255)
                                                                                                                                            cTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCa         >  1:985185/1‑150 (MQ=255)
                                                                                                                                             tGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGTTGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCat        <  2:92378/150‑1 (MQ=255)
                                                                                                                                             tGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCat        >  1:1604243/1‑150 (MQ=255)
                                                                                                                                               aaCGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGGTGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCatat      <  2:1238022/150‑1 (MQ=255)
                                                                                                                                                 cGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCatatat    >  1:1858744/1‑150 (MQ=255)
                                                                                                                                                 cGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCatatat    <  1:1549480/150‑1 (MQ=255)
                                                                                                                                                  gACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATAtt   <  1:1451888/150‑1 (MQ=255)
                                                                                                                                                   acacGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATAttt  >  2:1022133/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTT  >  NC_000913/2777424‑2777720

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