Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,764,096 A→T I353N (ATC→AAC)  yfjL ← CP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,764,0960AT98.7% 314.3 / ‑3.6 79I353N (ATC→AAC) yfjLCP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (40/38);  total (41/38)
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.10e-01

ATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATCACCTTCTGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCC  >  NC_000913/2763948‑2764242
                                                                                                                                                    |                                                                                                                                                  
atCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATCACCTTCTGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTc                                                                                                                                                  <  2:1161242/151‑1 (MQ=255)
        ccAACTCGCTGAAGCCAATCCAAAGCTGATCACCTTCTGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTg                                                                                                                                          >  2:436878/1‑151 (MQ=255)
          aaCTCGCTGAAGCCAATCCAAAGCTGATCACCTTCTGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGc                                                                                                                                         >  2:335340/1‑150 (MQ=255)
             tCGCTGAAGCCAATCCAAAGCTGATCACCTTCTGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAAc                                                                                                                                      >  2:759220/1‑150 (MQ=255)
                  gAAGCCAATCCAAAGCTGATCACCTTCTGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCa                                                                                                                                  <  2:699420/149‑1 (MQ=255)
                  gAAGCCAATCCAAAGCTGATCACCTTCTGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCAc                                                                                                                                 <  2:651005/150‑1 (MQ=255)
                         aTCCAAAGCTGATCACCTTCTGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGAccc                                                                                                                          >  2:620237/1‑150 (MQ=255)
                          tCCAAAGCTGATCACCTTCTGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCa                                                                                                                         >  1:675117/1‑150 (MQ=255)
                            cAAAGCTGATCACCTTCTGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCaaa                                                                                                                       >  1:1072038/1‑150 (MQ=255)
                                gCTGATCACCTTCTGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTc                                                                                                                  <  2:670985/151‑1 (MQ=255)
                                     tCACCTTCTGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTg                                                                                                              >  1:1123111/1‑150 (MQ=255)
                                         cTTCTGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTc                                                                                                         >  1:165237/1‑151 (MQ=255)
                                           tCTGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCa                                                                                                        <  1:721995/150‑1 (MQ=255)
                                             tGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATgg                                                                                                     >  1:403004/1‑151 (MQ=255)
                                              ggATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATgg                                                                                                     >  1:480846/1‑150 (MQ=255)
                                              ggATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATg                                                                                                      <  2:516810/149‑1 (MQ=255)
                                                   gcggcgAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGAt                                                                                               <  2:905598/151‑1 (MQ=255)
                                                     ggcgAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGAttt                                                                                             <  2:407523/151‑1 (MQ=255)
                                                      gcgAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTg                                                                                            <  2:838110/151‑1 (MQ=255)
                                                       cgAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTg                                                                                            >  2:502326/1‑150 (MQ=255)
                                                       cgAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGa                                                                                           <  2:358094/151‑1 (MQ=255)
                                                         aaCGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGa                                                                                           >  2:1006007/1‑149 (MQ=255)
                                                             atatGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACgg                                                                                     <  2:329503/151‑1 (MQ=255)
                                                             atatGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACgg                                                                                     >  1:235295/1‑151 (MQ=255)
                                                               atGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGa                                                                                    >  2:834266/1‑150 (MQ=255)
                                                                 gCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAAtt                                                                                 >  1:1154925/1‑151 (MQ=255)
                                                                  cACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTc                                                                                >  2:760214/1‑151 (MQ=255)
                                                                   aCGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCt                                                                               >  1:102057/1‑151 (MQ=255)
                                                                    cGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCtt                                                                              <  2:1200474/151‑1 (MQ=255)
                                                                      ttttCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCttt                                                                             >  1:501838/1‑150 (MQ=255)
                                                                       tttCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCttt                                                                             >  1:1113646/1‑149 (MQ=255)
                                                                        ttCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGt                                                                          >  2:690849/1‑151 (MQ=255)
                                                                          cAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCg                                                                        <  2:171893/151‑1 (MQ=255)
                                                                           aaGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGt                                                                       >  2:23629/1‑151 (MQ=255)
                                                                            aGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGtt                                                                      <  1:121595/151‑1 (MQ=255)
                                                                            aGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGtt                                                                      <  2:712432/151‑1 (MQ=255)
                                                                             gAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTtc                                                                     <  2:373120/151‑1 (MQ=255)
                                                                               aTTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTtctc                                                                   <  2:245055/151‑1 (MQ=255)
                                                                               aTTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTtctc                                                                   <  2:1100706/151‑1 (MQ=255)
                                                                                   tttCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCt                                                                <  2:948333/150‑1 (MQ=255)
                                                                                    ttCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCtt                                                               <  2:893137/150‑1 (MQ=255)
                                                                                      cGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTc                                                            >  1:968352/1‑151 (MQ=255)
                                                                                      cGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTc                                                            >  1:284589/1‑151 (MQ=255)
                                                                                       ggCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCa                                                           <  1:964522/151‑1 (MQ=255)
                                                                                          cccTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCaatg                                                        <  1:573505/151‑1 (MQ=255)
                                                                                            cTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCaatga                                                       >  1:706428/1‑150 (MQ=255)
                                                                                             tGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCaatgaat                                                     >  2:814126/1‑151 (MQ=255)
                                                                                                 ttttCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGGGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTc                                                  >  2:339720/1‑150 (MQ=255)
                                                                                                 ttttCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTc                                                  >  2:912690/1‑150 (MQ=255)
                                                                                                  tttCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCa                                                 <  2:322180/150‑1 (MQ=255)
                                                                                                    tCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGt                                               <  2:843112/150‑1 (MQ=255)
                                                                                                        cGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTGCGTTTATg                                          >  1:850052/1‑151 (MQ=255)
                                                                                                        cGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATg                                          >  2:870691/1‑151 (MQ=255)
                                                                                                        cGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATg                                          >  1:833878/1‑151 (MQ=255)
                                                                                                         gTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGa                                         <  2:1121297/151‑1 (MQ=255)
                                                                                                           ataGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAat                                       >  2:611335/1‑151 (MQ=255)
                                                                                                               cATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCa                                   <  2:399959/151‑1 (MQ=255)
                                                                                                               cATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCa                                   <  2:1184375/151‑1 (MQ=255)
                                                                                                                aTGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATc                                    >  1:540097/1‑149 (MQ=255)
                                                                                                                aTGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATc                                    >  1:229320/1‑149 (MQ=255)
                                                                                                                aTGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAg                                  <  2:12935/151‑1 (MQ=255)
                                                                                                                  gCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGAc                                <  1:538237/151‑1 (MQ=255)
                                                                                                                       gAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTc                           <  1:994546/151‑1 (MQ=255)
                                                                                                                       gAATAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTc                           <  2:1071181/151‑1 (MQ=255)
                                                                                                                        aaTAGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCt                          <  1:486260/151‑1 (MQ=255)
                                                                                                                         aTAGGCTGGCTTTCGTGGTGCGGCCGGTTCGTTGCGTGCAACAGTACTCGACCCAAAAATTCGCTGAAGTCATGGAGCGTTTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCtt                         >  1:375434/1‑151 (MQ=255)
                                                                                                                           aGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCtttt                       <  2:506589/151‑1 (MQ=255)
                                                                                                                           aGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCtttt                       >  1:698976/1‑151 (MQ=255)
                                                                                                                           aGGCTGGCTTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCttt                        <  2:1128359/150‑1 (MQ=255)
                                                                                                                            ggctggctTTCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTg                      >  2:26930/1‑151 (MQ=255)
                                                                                                                                   tttCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCgg                 >  1:188106/1‑149 (MQ=255)
                                                                                                                                    ttCGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCgg                 >  2:536111/1‑148 (MQ=255)
                                                                                                                                      cGTGGTGCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGaaa            <  1:1122887/151‑1 (MQ=255)
                                                                                                                                       gtggtgCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGaaaa           <  1:1264715/151‑1 (MQ=255)
                                                                                                                                          gtgCTGCCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAAcc        >  2:1095592/1‑151 (MQ=255)
                                                                                                                                           tgctgcCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCt       >  2:583245/1‑151 (MQ=255)
                                                                                                                                            gctgcCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCt       <  2:1038142/150‑1 (MQ=255)
                                                                                                                                               gcCGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCa     <  2:252813/149‑1 (MQ=255)
                                                                                                                                                ccGGTTCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATcc  >  1:139677/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ATCAGTCACCAACTCGCTGAAGCCAATCCAAAGCTGATCACCTTCTGGATCGCGGCGAACGATATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCC  >  NC_000913/2763948‑2764242

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