Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,560,330 G→T 16.2% C25F (TGC→TTC)  yffL → CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,560,3300GT16.2% 274.7 / 38.0 108C25F (TGC→TTC) yffLCPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (47/42);  new base T (11/7);  total (58/50)
Fisher's exact test for biased strand distribution p-value = 6.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

Only 100 of 109 total aligned reads displayed.

TTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAA  >  NC_000913/2560234‑2560429
                                                                                                |                                                                                                   
ttCAATATCAAGCACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCg                                                                                                  >  2:3007606/1‑100 (MQ=255)
ttCAATATCAAGAACTTAAATCCATGTTTGCAACGAAATATCCATAGTGCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCg                                                                                                  >  2:605382/1‑100 (MQ=255)
ttCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCg                                                                                                  >  1:338244/1‑100 (MQ=255)
  cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                                <  2:71226/100‑1 (MQ=38)
   aaTATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCgggg                                                                                               <  1:2639775/100‑1 (MQ=255)
    atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggggg                                                                                              <  2:1452196/100‑1 (MQ=38)
    atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCACCATCCCggggg                                                                                              <  2:1580890/100‑1 (MQ=255)
     tatTAAGAACTTAAATCCATGTTTGCAACTAAAGATGCAGTGTTCGAAAATAGGATCAACACAAATAAACGCCCCCGGTATGCTGCAACATGCCGGGGGc                                                                                             <  2:2197576/100‑1 (MQ=255)
     tatCAAGAACTTAAATCCATGTTTGCATCGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgggg                                                                                               <  2:1482719/98‑1 (MQ=38)
     tatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGc                                                                                             <  1:2608539/100‑1 (MQ=255)
     tatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGc                                                                                             <  2:2760649/100‑1 (MQ=255)
       tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                           >  1:16201/1‑100 (MQ=255)
       tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                           <  1:2547103/100‑1 (MQ=255)
       tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                           <  1:1671843/100‑1 (MQ=255)
       tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                           >  2:2215069/1‑100 (MQ=255)
       tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                           >  2:226519/1‑100 (MQ=255)
       tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                           <  2:878116/100‑1 (MQ=255)
        cAAGAACTTAAATCCATGTTTGCAACGAAAGATCCATATTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGt                                                                                          <  1:773013/100‑1 (MQ=255)
        cAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCGGGGGcttt                                                                                          <  1:1186561/100‑4 (MQ=37)
        cAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGt                                                                                          <  2:230208/100‑1 (MQ=255)
         aaGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATATAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTa                                                                                         <  1:2139459/100‑1 (MQ=255)
         aaGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTa                                                                                         <  2:2105753/100‑1 (MQ=255)
          aGATCTTAAATCCATGTTTGCAACGAAATTTCTAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTa                                                                                         <  1:768886/99‑1 (MQ=255)
           gAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATg                                                                                       <  2:207440/100‑1 (MQ=255)
            aaCTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATga                                                                                      <  2:963513/100‑1 (MQ=255)
             aCTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgag                                                                                     <  2:2807323/100‑1 (MQ=255)
              cTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgaga                                                                                    <  2:1897239/100‑1 (MQ=255)
              cTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgaga                                                                                    <  1:975403/100‑1 (MQ=255)
                tAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGaaa                                                                                  <  1:2107541/100‑1 (MQ=255)
                tAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGaa                                                                                   <  1:2289207/99‑1 (MQ=255)
                tAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCAGGAATGCTGCAACATGCCGGGGGCGGTATGAGaaa                                                                                  >  1:1125525/1‑100 (MQ=255)
                    tCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAg                                                                              <  1:419999/100‑1 (MQ=255)
                     ccATGTTTGCAACGAAAGCTCCAGAGTTCGAAAATAGGATCAACTCAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACaa                                                                               >  1:159640/1‑98 (MQ=255)
                      cATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGc                                                                             >  1:1539497/1‑99 (MQ=255)
                      cATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCAAGGGGCGATATGAGAAACAAg                                                                              >  1:2729343/1‑98 (MQ=255)
                       aTGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTa                                                                           <  2:1341916/100‑1 (MQ=255)
                          tttGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAg                                                                        <  1:559575/100‑1 (MQ=255)
                               aaCGAAAGATCGAGAGTTCGAAAATAGGAGCAACACAAATAAAAGGCCCCGGAATGCTGCAAAATGCCGGGGGCGGTATGAGAAACAAGCTAAAGgag                                                                     <  2:829086/98‑1 (MQ=255)
                               aaCGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGaa                                                                   >  1:1752410/1‑100 (MQ=255)
                               aaCGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCACCCGGAATGCTGAAACATGACGGGGGCGGTATGAGAAACAAGCTAAAGGAGaa                                                                   >  1:97838/1‑100 (MQ=255)
                                   aaaGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtc                                                               >  1:1320123/1‑100 (MQ=255)
                                     aGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtctc                                                             <  1:1800455/100‑1 (MQ=255)
                                      gATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCa                                                            >  1:2713404/1‑100 (MQ=255)
                                      gATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCa                                                            >  2:116937/1‑100 (MQ=255)
                                        tCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATg                                                          <  1:2796150/100‑1 (MQ=255)
                                          cAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTc                                                        >  2:1693158/1‑100 (MQ=255)
                                          cAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTc                                                        >  2:1129036/1‑100 (MQ=255)
                                              gTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                                    >  1:1691371/1‑100 (MQ=255)
                                              gTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                                    >  1:1864520/1‑100 (MQ=255)
                                                tCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGatat                                                  <  1:2614471/100‑1 (MQ=255)
                                                 cGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATg                                                 >  2:2029065/1‑100 (MQ=255)
                                                 cGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATg                                                 >  1:147793/1‑100 (MQ=255)
                                                  gAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGAGGGCGGTATGAGAAACAAGATAAAGGAGAATTTCTCATGTCTGATATGc                                                >  2:2068066/1‑100 (MQ=255)
                                                   aaaaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGACGATATGAGAAAAAAAAAAAAGGAGAATTTCTCATGTCTGATATGc                                                >  1:769607/1‑99 (MQ=255)
                                                    aaaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTa                                              <  1:208521/100‑1 (MQ=255)
                                                     aaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAg                                             <  1:1193736/100‑1 (MQ=255)
                                                     aaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAg                                             >  1:973332/1‑100 (MQ=255)
                                                     aaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAg                                             <  1:658077/100‑1 (MQ=255)
                                                      aTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGTAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGc                                            <  1:2573482/100‑1 (MQ=255)
                                                        aGGATCAACACGAATAAAATCCCCCGGAATGCTGCAACATTCCGTGGGCGGTATGAGAAACAAGCTAAAGGATAATTTCTCATGTCTGATATGCTAGc                                            <  2:2876688/98‑1 (MQ=255)
                                                         ggATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGt                                         >  2:1371032/1‑100 (MQ=255)
                                                          gATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTa                                        <  2:810458/100‑1 (MQ=255)
                                                           aTCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGCATTTCTCATGTCTGATATGCTAGCAGTa                                        >  2:2573164/1‑99 (MQ=255)
                                                            tCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGa                                      <  1:149960/100‑1 (MQ=255)
                                                               acacAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACa                                   >  2:2553109/1‑100 (MQ=255)
                                                                  tctcataccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGa                                >  2:1390927/10‑100 (MQ=25)
                                                                  caAATAAAAGCCCCTGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGa                                >  2:1635763/1‑100 (MQ=255)
                                                                   aaaTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaa                               <  1:226521/100‑1 (MQ=255)
                                                                     atctaccgccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaac                             <  1:258244/92‑1 (MQ=37)
                                                                       aaaaGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaa                           >  2:644025/1‑100 (MQ=255)
                                                                         ccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATATCTGATATGATAGCAGTAGAACAAGaaacaaac                         >  1:525121/3‑100 (MQ=37)
                                                                         accgccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacagat                         >  2:258241/5‑97 (MQ=37)
                                                                          cgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAaacaaaca                        >  2:1894762/2‑100 (MQ=38)
                                                                           gCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAacaaacaa                       >  1:727572/1‑100 (MQ=38)
                                                                           gCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTATCATGTCAGATATGCTAGCAGTAGAACAAGAAacaaacaa                       >  1:640718/1‑100 (MQ=37)
                                                                            cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat                      >  2:2844295/1‑100 (MQ=38)
                                                                            cccccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat                      >  2:2437344/1‑100 (MQ=255)
                                                                            cccccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat                      <  1:1289291/100‑1 (MQ=255)
                                                                             ccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatg                     >  2:2491012/1‑100 (MQ=38)
                                                                             ccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatg                     >  1:2937318/1‑100 (MQ=38)
                                                                             ccccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatg                     >  1:228540/1‑100 (MQ=255)
                                                                              cccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatga                    >  1:1292819/1‑100 (MQ=38)
                                                                              cccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatga                    >  2:854282/1‑100 (MQ=255)
                                                                              cccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatga                    <  1:661265/100‑1 (MQ=255)
                                                                              cccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatga                    >  1:807134/1‑100 (MQ=255)
                                                                              cacGGAATGCTGCAACATTCCGGGGGAGGTATGAGAAACAAGCTAAAGCAGAATTTCTCATGTCTGATATGCTATCAGTAGAACAAGAAACAAACAat                      <  2:2743487/96‑1 (MQ=255)
                                                                               ccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatgat                   >  2:1733185/1‑100 (MQ=38)
                                                                               ccGGAATGCTGCAACCTGCCGGGAGAGGTATGAGAAACACGCTAACGGAGAATTCCTCATGTCTGATATGTTAGCAGTAGAACAAGaatcaaactg                       <  2:1763573/96‑3 (MQ=255)
                                                                               ccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatgat                   >  2:1521263/1‑100 (MQ=255)
                                                                                      gctgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTc            >  2:547955/1‑100 (MQ=255)
                                                                                       ctgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCa           >  1:2944616/1‑100 (MQ=255)
                                                                                        tgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCaa          >  2:875598/1‑100 (MQ=255)
                                                                                         gcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAACAGTAGAACAAGAAACAAACAATGATGTTCGTCAAt         >  1:74245/1‑100 (MQ=255)
                                                                                          cAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAAtt        <  2:956803/100‑1 (MQ=255)
                                                                                          cAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAAAAAGAAACAAACAATGATGTTCGTCAAtt        >  1:3040953/1‑100 (MQ=255)
                                                                                           aaCATGCCGGGGGCGGTATGAGAAACGAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAAttt       >  1:1839768/1‑100 (MQ=255)
                                                                                              aTGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTa    <  1:1521281/100‑1 (MQ=255)
                                                                                               tGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCACTTTTTaa   >  1:2161096/1‑100 (MQ=255)
                                                                                               tGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTaa   >  1:1481047/1‑100 (MQ=255)
                                                                                                gCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTaaa  >  2:432687/1‑100 (MQ=255)
                                                                                                |                                                                                                   
TTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAA  >  NC_000913/2560234‑2560429

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: