Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,558,855 G→T 5.2% intergenic (‑176/‑3) eutB ← / → intZ ethanolamine ammonia‑lyase, large subunit, heavy chain/CPZ‑55 prophage; putative phage integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,558,8550GT 5.2% 700.3 / 19.4 193intergenic (‑176/‑3)eutB/intZethanolamine ammonia‑lyase, large subunit, heavy chain/CPZ‑55 prophage; putative phage integrase
Reads supporting (aligned to +/- strand):  ref base G (88/95);  new base T (4/6);  total (92/101)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 195 total aligned reads displayed.

ATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGTTTTTGCTGGTGCATCCCAACGGTTCGAAATACTGGCGATTCAGGTATTGGCT  >  NC_000913/2558707‑2559003
                                                                                                                                                    |                                                                                                                                                    
    gCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGtttt                                                                                                                                                  >  4:1428450/1‑149 (MQ=255)
       tttttCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAg                                                                                                                                               >  4:1443347/1‑149 (MQ=255)
          ttCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCaa                                                                                                                                            <  2:464253/149‑1 (MQ=255)
               aGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGa                                                                                                                                          >  2:1134578/1‑146 (MQ=255)
                      tAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTcac                                                                                                                                >  1:1294592/1‑149 (MQ=255)
                       aaCAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTcaca                                                                                                                               >  3:25232/1‑149 (MQ=255)
                            tCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACCTGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCaaa                                                                                                                          >  2:515758/1‑149 (MQ=255)
                            tCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCaaa                                                                                                                          <  3:1403805/149‑1 (MQ=255)
                            tCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCaaa                                                                                                                          <  3:1313002/149‑1 (MQ=255)
                                 tttCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGa                                                                                                                                                <  1:537928/122‑1 (MQ=255)
                                 tttCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGa                                                                                                                                                >  2:537928/1‑122 (MQ=255)
                                 tttCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCa                                                                                                                     <  3:278/149‑1 (MQ=255)
                                  ttCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACGTGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCAt                                                                                                                    <  4:374922/149‑1 (MQ=255)
                                        gTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGAcgc                                                                                                              >  3:1433305/1‑149 (MQ=255)
                                        gTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGAcgc                                                                                                              <  2:599348/149‑1 (MQ=255)
                                           aGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCa                                                                                                                     <  1:1234949/139‑1 (MQ=255)
                                            gCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGca                                                                                                          >  3:1463107/1‑149 (MQ=255)
                                             cccTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGa                                                                                                                 >  2:960051/1‑141 (MQ=255)
                                              ccTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGcaca                                                                                                        >  4:1204013/1‑149 (MQ=255)
                                               cTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTcaca                                                                                                                               >  3:617667/1‑125 (MQ=255)
                                               cTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTcaca                                                                                                                               <  4:617667/125‑1 (MQ=255)
                                                ttGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGc                                                                                                      >  4:370361/1‑149 (MQ=255)
                                                  gTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCg                                                                                                    >  4:209934/1‑149 (MQ=255)
                                                     gtgtAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGTTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAg                                                                                                 <  2:724423/149‑1 (MQ=255)
                                                     gtgtAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAg                                                                                                 <  1:953479/149‑1 (MQ=255)
                                                      tgtAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGc                                                                                                      >  1:1061776/1‑143 (MQ=255)
                                                      tgtAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGc                                                                                                      <  2:1061776/143‑1 (MQ=255)
                                                         aaGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACa                                                                                             >  3:968503/1‑149 (MQ=255)
                                                          aGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACaa                                                                                            >  2:1423962/1‑149 (MQ=255)
                                                            gtggtTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCaa                                                                                                                           <  1:593993/116‑1 (MQ=255)
                                                            gtggtTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCaa                                                                                                                           >  2:593993/1‑116 (MQ=255)
                                                            gtggtTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGc                                                                                          <  3:1193707/149‑1 (MQ=255)
                                                            gtggtTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGc                                                                                          <  3:1210796/149‑1 (MQ=255)
                                                             tggtTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGcc                                                                                         <  4:25232/149‑1 (MQ=255)
                                                                ttGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCta                                                                                      <  3:843445/149‑1 (MQ=255)
                                                                   ccTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGa                                                                                                                                                <  4:1024920/88‑1 (MQ=255)
                                                                   ccTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGa                                                                                                                                                >  3:1024920/1‑88 (MQ=255)
                                                                    cTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGAcgccg                                                                                                            <  3:432985/123‑1 (MQ=255)
                                                                    cTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATaaa                                                                                  >  2:219149/1‑149 (MQ=255)
                                                                     tAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAgg                                                                                                                        <  1:684408/110‑1 (MQ=255)
                                                                     tAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAAc                                                                                 >  4:408329/1‑149 (MQ=255)
                                                                       aaCACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTc                                                                               <  3:370361/149‑1 (MQ=255)
                                                                         cacGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACGCAc                                                                             <  1:515758/149‑1 (MQ=255)
                                                                             cccGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGa                                                                                                                                                <  4:794184/78‑1 (MQ=255)
                                                                             cccGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGAt                                                                         >  4:440503/1‑149 (MQ=255)
                                                                              ccGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGacac                                                                                                                                       <  3:4488/86‑1 (MQ=255)
                                                                               cGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCaa                                                                                                                                            >  3:568571/1‑80 (MQ=255)
                                                                               cGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCaa                                                                                                                                            <  4:568571/80‑1 (MQ=255)
                                                                                gTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCAc                                                                             <  1:401990/142‑1 (MQ=255)
                                                                                gTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCa                                                                      <  1:219149/149‑1 (MQ=255)
                                                                                     tcatcaCGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGAtt                                                                        <  1:488941/142‑1 (MQ=255)
                                                                                     tcatcaCGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGAtt                                                                        >  2:488941/1‑142 (MQ=255)
                                                                                     tcatcaCGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGAtt                                                                        >  2:460677/1‑142 (MQ=255)
                                                                                        tcaCGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGTTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGcc                                                              >  2:385715/1‑149 (MQ=255)
                                                                                        tcaCGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAg                                                                                                                         >  2:654695/1‑90 (MQ=255)
                                                                                        tcaCGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAg                                                                                                                         <  1:654695/90‑1 (MQ=255)
                                                                                        tcaCGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCAc                                                                     <  1:462687/142‑1 (MQ=255)
                                                                                        tcaCGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGcc                                                              >  4:1265113/1‑149 (MQ=255)
                                                                                          aCGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCgg                                                            >  2:302555/1‑149 (MQ=255)
                                                                                            gCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGc                                                                                                                       <  1:863837/88‑1 (MQ=255)
                                                                                            gCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGc                                                                                                                       >  2:863837/1‑88 (MQ=255)
                                                                                             ccGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGc                                                         <  4:968503/149‑1 (MQ=255)
                                                                                              cGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCt                                                        >  4:668912/1‑149 (MQ=255)
                                                                                               gAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTg                                                       <  3:209934/149‑1 (MQ=255)
                                                                                                aaCGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGt                                                      <  1:357825/149‑1 (MQ=255)
                                                                                                   gcgcCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGtttt                                                   <  3:1265113/149‑1 (MQ=255)
                                                                                                     gcCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGc                                                                                                      >  2:1208267/1‑96 (MQ=255)
                                                                                                      ccGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGTTTTTGc                                                <  1:302555/149‑1 (MQ=255)
                                                                                                        ggATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGTTTTTGCTg                                              <  1:1316706/149‑1 (MQ=255)
                                                                                                           tACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATaaa                                                                                  >  2:641766/1‑110 (MQ=255)
                                                                                                             cATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGTTTTTGCTGGTGCa                                         >  2:67835/1‑149 (MQ=255)
                                                                                                              aTGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGAcaccac                                                                                                                                    <  4:1213011/57‑1 (MQ=255)
                                                                                                              aTGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCAc                                                                     >  2:316442/1‑120 (MQ=255)
                                                                                                              aTGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGtt                                                     >  2:264536/1‑136 (MQ=255)
                                                                                                              aTGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGtt                                                     <  1:264536/136‑1 (MQ=255)
                                                                                                                 aTCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCAc                                                                     >  2:1399424/1‑117 (MQ=255)
                                                                                                                 aTCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCgggg                                                          <  3:18101/128‑1 (MQ=255)
                                                                                                                 aTCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCgggg                                                          >  4:18101/1‑128 (MQ=255)
                                                                                                                 aTCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGTTTTTGCTGGTGCATc                                       >  3:404637/1‑147 (MQ=255)
                                                                                                                     gggTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGTTTTTGCTg                                              >  1:1143801/1‑136 (MQ=255)
                                                                                                                        ttATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGTTTTTGCTGGTGCATCCCAACGGtt                              >  3:352041/1‑149 (MQ=255)
                                                                                                                          aTCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAg                                                                                                 >  1:147592/1‑80 (MQ=255)
                                                                                                                            ccAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAg                                                                                                 <  1:963911/78‑1 (MQ=255)
                                                                                                                              aGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGTTTTTGCTGGTGCATCCCAACGGTTCGAAAt                        >  1:1357425/1‑149 (MQ=255)
                                                                                                                               gTCGTTAAATCAAGGTATCCGTTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGTTTTTGCTGGTGCATCCCAAc                                  <  1:112873/138‑1 (MQ=255)
                                                                                                                               gTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGTTTTTGCTGGTGCATCCCAACGGTTCGAAATa                       <  4:1317353/149‑1 (MQ=255)
                                                                                                                                  gTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGTTTTTGCTg                                              >  2:129503/1‑123 (MQ=255)
                                                                                                                                  gTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGTTTTTGCTGGTGCATCCCAACGGTTCGAAATACTg                    <  1:67835/149‑1 (MQ=255)
                                                                                                                                  gTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGTTTTTGCTGGTGCATCCCAACGGTTCGAAATACTg                    >  1:645458/1‑149 (MQ=255)
                                                                                                                                   ttAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAg                                                                                                                              >  4:1464323/1‑42 (MQ=255)
                                                                                                                                   ttAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAg                                                                                                                              <  3:1464323/42‑1 (MQ=255)
                                                                                                                                   ttAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGTTTTTGCTGGTGCATCCCAACGGTTCGAAATACTgg                   <  4:352041/149‑1 (MQ=255)
                                                                                                                                   ttAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGTTTTTGCTGGTGCATCCCAACGGTTCGAAATACTgg                   <  1:141798/149‑1 (MQ=255)
                                                                                                                                   ttAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGTTTTTGCTGGTGCATCCCAACGGTTCGAAATACTgg                   >  3:68774/1‑149 (MQ=255)
                                                                                                                                      aaTCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGTTTTTGCTGGt                                            <  3:949081/121‑1 (MQ=255)
                                                                                                                                      aaTCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGTTTTTGCTGGt                                            >  4:949081/1‑121 (MQ=255)
                                                                                                                                            ggTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCCGGGCTGTTTTTGCTGGTGCATCCCAACGGTTCGAAATACTGGCGATTCAgg          <  2:1249741/149‑1 (MQ=255)
                                                                                                                                                tCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGTTTTTGCTGGTGCATCCCAAc                                  <  4:1182045/121‑1 (MQ=255)
                                                                                                                                                tCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGTTTTTGCTGGTGCATCCCAAc                                  >  3:1182045/1‑121 (MQ=255)
                                                                                                                                                tCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGTTTTTGCTGGTGCATCCCAACGGTTCGAAATACTGGCGATTCAg           <  1:609108/144‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCGTCAGCCCTTGTCGTGTAAGGTGGTTGCCTAAACACGCCCGTTATTCATCACGCCGAACGCGCCGGATACATGATCGGGGTTATCCAGTCGTTAAATCAAGGTATCCGGTTTTGAGCAAGACACCACTCACAGCAAAGGCCATAGACGCCGCACAGCCGCAGGACAAGCCCTATAAACTCACCGATTCACTGACGCCGGGGCTGTTTTTGCTGGTGCATCCCAACGGTTCGAAATACTGGCGATTCAGGTATTGGCT  >  NC_000913/2558707‑2559003

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