BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	968,412	T→A	5.4%	L16* (TTG→TAG)	lpxK →	lipid A 4'kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	968,412	0	T	A	5.4%	352.7 / 4.5	117	L16* (TTG→TAG)	lpxK	lipid A 4'kinase
Reads supporting (aligned to +/- strand): ref base T (50/60); new base A (4/3); total (54/63)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

CGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCGGTCGTTG  >  NC_000913/968264‑968561

 gTGGTCGTCGAGGATGGTGTCTTTGTGGAACGCGGGAATCATAACGATTTGCTTGAGCTCCGCTGCGTTTTCGAGCATCTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTTTtgg                                                                                                                                                     <  1:1578561/150‑2 (MQ=255)

 gTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATtgc                                                                                                                                                     <  1:125632/150‑1 (MQ=255)

   ggTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATtgctg                                                                                                                                                   <  2:1545933/150‑1 (MQ=255)

    gtcgtcGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATtgctgc                                                                                                                                                  <  2:1341945/150‑1 (MQ=255)

     tcgtcgAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGcc                                                                                                                                                 >  2:242516/1‑150 (MQ=255)

     tcgtcgAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGcc                                                                                                                                                 <  2:1653496/150‑1 (MQ=255)

     gcgtcgAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGGCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGcc                                                                                                                                                 <  2:1034468/149‑1 (MQ=255)

       gtcgAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCATAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCAc                                                                                                                                               <  2:433640/150‑1 (MQ=255)

       ggcgAGGATGGTGTCTTTGTGGATCGCGGTACGCATAACGATTTGGTTGAGCACCGCGGCGTTTACGCGCATCTTCACAAAATGCAGTTTGGACAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTTGCAGCAATAGCCGCCAc                                                                                                                                               <  2:420181/148‑1 (MQ=255)

        tcgAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCAct                                                                                                                                              >  1:1740767/1‑150 (MQ=255)

         cgAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCActc                                                                                                                                             >  2:638337/1‑150 (MQ=255)

          gAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCACCTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCActct                                                                                                                                            >  1:495190/1‑150 (MQ=255)

          gAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCActca                                                                                                                                            >  1:1671380/1‑149 (MQ=255)

          gAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCAGCAATAGCCGCCAcaaa                                                                                                                                            <  1:998343/150‑4 (MQ=255)

            ggATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTcc                                                                                                                                          >  2:71206/1‑150 (MQ=255)

              aTGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTg                                                                                                                                        <  1:1851009/150‑1 (MQ=255)

               tGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCGAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTgg                                                                                                                                       >  1:959339/1‑150 (MQ=255)

                 gtgtCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGtt                                                                                                                                     >  1:361621/1‑150 (MQ=255)

                 gtgtCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTGTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGtt                                                                                                                                     <  2:179201/150‑1 (MQ=255)

                  tgtCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTg                                                                                                                                    <  2:1377926/150‑1 (MQ=255)

                    tCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTa                                                                                                                                  >  1:1477708/1‑150 (MQ=255)

                    tCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTa                                                                                                                                  <  2:1190523/150‑1 (MQ=255)

                      aTTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATg                                                                                                                                >  1:474478/1‑150 (MQ=255)

                         gtgGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGcc                                                                                                                             <  1:871831/150‑1 (MQ=255)

                           ggAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTg                                                                                                                           >  1:906534/1‑150 (MQ=255)

                           ggAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTg                                                                                                                           >  2:1414068/1‑150 (MQ=255)

                            gcaCGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGCTTACGCGCAACTTCACAAAATGCACTTTGGCCAATGATCGAAAAAATCTGGGCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCggg                                                                                                                          <  2:1260621/148‑1 (MQ=255)

                            gAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTgg                                                                                                                          <  1:1619557/150‑1 (MQ=255)

                              aCGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTg                                                                                                                        <  1:839451/150‑1 (MQ=255)

                              aCGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGACTGGtt                                                                                                                        >  1:317130/1‑149 (MQ=255)

                              aCGCGCTACGCATAACGATTTGCTAGAGCACCGCGGCGTTTACGCGAAACTTCCCAAAATGCAGTTTGGCCAATGATCGAAAATAGCTGGTCTGGTGACTCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTg                                                                                                                        <  1:143476/150‑1 (MQ=255)

                                       gCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAGAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGcgcg                                                                                                               >  2:661240/1‑150 (MQ=255)

                                       gCATAACGATTTGCTTGAGCACCGCGGCGGTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGcgcg                                                                                                               <  2:1694419/150‑1 (MQ=255)

                                        cATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAACTGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATAGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGa                                                                                                              >  2:1312834/1‑150 (MQ=255)

                                        cATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAAGCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGa                                                                                                              >  2:1184627/1‑150 (MQ=255)

                                         aTAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTGGTATGGCCTGGTGAGTGGCGCGAt                                                                                                             <  2:1455363/150‑1 (MQ=255)

                                          tAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATc                                                                                                            >  1:725845/1‑150 (MQ=255)

                                             cGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGACTGGTGAGTGGCGCGATCCGt                                                                                                         >  1:845435/1‑150 (MQ=255)

                                               tgttGCTTGAGCTCCGCGGCGTTTACGCGCATCTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTTAGTGGCGCGATCCGTCt                                                                                                       <  1:1829413/148‑1 (MQ=255)

                                               aTTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCt                                                                                                       >  2:1025193/1‑150 (MQ=255)

                                                tttGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCCGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCtt                                                                                                      <  2:1364358/150‑1 (MQ=255)

                                                 ttGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCttt                                                                                                     <  2:1740767/150‑1 (MQ=255)

                                                  tGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTg                                                                                                    >  1:1880709/1‑150 (MQ=255)

                                                   gCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGAGCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCGGGTGAGTGGCGCGATCCGTCTTTGc                                                                                                   <  1:37161/150‑1 (MQ=255)

                                                   gCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGGTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGc                                                                                                   <  1:1139403/150‑1 (MQ=255)

                                                      tGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCtat                                                                                                >  2:874916/1‑150 (MQ=255)

                                                       gAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCtata                                                                                               <  1:751965/150‑1 (MQ=255)

                                                        aGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATaa                                                                                              >  2:797239/1‑150 (MQ=255)

                                                         gCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATaaa                                                                                             <  2:297104/150‑1 (MQ=255)

                                                             cgcgGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATATGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCAGGTTGTATGGCCTGGGGAGTGGCGCGATCCGTCGTTGCTATAAACTaa                                                                                         >  2:1170455/1‑150 (MQ=255)

                                                                 gcgTTTACGGGCAACTTCACAAAAGGCAGTTTGGACAAAGATCGAAAAAATCTGGTCTGGTGACTCCCATTTGTGGCGGCTATTGCTGCCACTCTCCTTGTTGTATGGACTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACt                                                                                     <  2:286942/150‑1 (MQ=255)

                                                                      tACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAgc                                                                                <  1:1368408/150‑1 (MQ=255)

                                                                           gcAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGcc                                                                           <  1:1423846/150‑1 (MQ=255)

                                                                             aaCTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTg                                                                         <  1:409313/150‑1 (MQ=255)

                                                                               cTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATTAACTAAAACTGAAGCGCGCCTGGc                                                                       >  1:277799/1‑150 (MQ=255)

                                                                                ttCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTTGTGAATCCCCTTTGTGGTGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCg                                                                      <  1:1442703/150‑1 (MQ=255)

                                                                                ttCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGTTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCg                                                                      <  2:929444/150‑1 (MQ=255)

                                                                                atCACAAAATCCAGTTACGCCAATGATCGAAAACATCTGGGCTGGTGAATCCCCTTTGTGGGGTCTTTTGCTTCCACTCTCCTGGCTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTACAACTGAAGCGCGCCTGGCg                                                                      <  2:1761588/149‑1 (MQ=255)

                                                                                 tCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGt                                                                     >  2:1323741/1‑150 (MQ=255)

                                                                                  cacaAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTg                                                                    >  1:773821/1‑150 (MQ=255)

                                                                                     aaaaTGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGccc                                                                 <  1:757594/150‑1 (MQ=255)

                                                                                      ataTGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGcccc                                                                <  2:1881875/148‑1 (MQ=255)

                                                                                      aaaTGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGcccc                                                                <  2:1675357/150‑1 (MQ=255)

                                                                                      aaaTGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGcccc                                                                <  2:1686581/150‑1 (MQ=255)

                                                                                       aaTGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGccccc                                                               <  2:1698699/150‑1 (MQ=255)

                                                                                       aaTGCAGTTTGCCCACTGATCGCAAAACTCTGGTCTGGTGAATCCCCTTTGGGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGccccc                                                               <  2:1283341/150‑1 (MQ=255)

                                                                                             gTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGTGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCg                                                         >  1:1388101/1‑150 (MQ=255)

                                                                                              tttGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCTCTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCgg                                                        <  1:701991/150‑1 (MQ=255)

                                                                                              tttGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCgg                                                        <  1:1522798/150‑1 (MQ=255)

                                                                                                tcgCCAATCATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGtt                                                      <  2:151277/148‑1 (MQ=255)

                                                                                                   ccACTGATCGAAAAAAACAGGCCTGGTGAATCCCCTTTCTCGCGGCTATTGCTGCCACTCTCCTGTTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGGTGTc                                                   <  1:115629/150‑1 (MQ=255)

                                                                                                   ccAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTc                                                   >  1:982696/1‑150 (MQ=255)

                                                                                                       tGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtgg                                               <  1:1365758/150‑1 (MQ=255)

                                                                                                        gATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGGACCGGTTGTCGtggt                                              >  1:1764420/1‑150 (MQ=255)

                                                                                                         aTCGAAAAAATCTGGTCTGGTGCATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTCTAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggtt                                             <  2:1747253/150‑1 (MQ=255)

                                                                                                         aTCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggtt                                             >  2:815673/1‑150 (MQ=255)

                                                                                                         aTCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCGATTGCTGCCACTCTCCTGGTTGTAGGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggtt                                             <  2:1137219/150‑1 (MQ=255)

                                                                                                          tCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggttg                                            >  1:1364590/1‑150 (MQ=255)

                                                                                                           cGAAAAAATCTGGTCTGGTGTATCCCCTTTGTGGCGGCTATTGCTGCGACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggttgg                                           <  1:778793/150‑1 (MQ=255)

                                                                                                             aaaaaaTCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCTTACCGGTTGTCTAGGTTGGTa                                         >  2:349386/1‑150 (MQ=255)

                                                                                                                aaaTCTTGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTTGTAAtc                                      <  1:1109180/150‑1 (MQ=255)

                                                                                                                      ggtctggtGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCg                                <  1:1910287/150‑1 (MQ=255)

                                                                                                                       ggctggtGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGc                               <  1:774516/148‑1 (MQ=255)

                                                                                                                           ggtGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAggc                           >  2:1575370/1‑150 (MQ=255)

                                                                                                                              gAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAggcggc                        >  2:153227/1‑150 (MQ=255)

                                                                                                                              gAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAggcggc                        <  2:474478/150‑1 (MQ=255)

                                                                                                                                aTCCCCTTTGTGTCGGCTATTGCTGCCACTCTCCAGGTTGTATGGCATGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGCTAATCTCACCGCACGCGGCaa                      >  1:780164/1‑150 (MQ=255)

                                                                                                                                 tCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTAGCGGTTGTCGTGGTTGGTTATCTCACCGCAGGCGGCAAc                     >  1:1746972/1‑150 (MQ=255)

                                                                                                                                  caggagagtgGCAGCAATAGCCGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACg                    >  2:1855980/8‑150 (MQ=255)

                                                                                                                                  caggagagtgGCAGCAATAGCCGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCATGGCGTGCCCCCGTACCGGTTGTCGTGGGTGGTAATCTCACCGCAGGCGGCAACg                    >  2:239468/8‑150 (MQ=255)

                                                                                                                                   cccTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCTGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACgg                   >  1:987381/1‑150 (MQ=255)

                                                                                                                                   cccTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACgg                   >  1:1534731/1‑150 (MQ=255)

                                                                                                                                   cccTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTAGTCGTGGTTGGTAATCTCACCGCAGGCCGCAACgg                   >  2:1836148/1‑150 (MQ=255)

                                                                                                                                   cccTTTGGGGCGGCTATTGCTGCCACTCTCCTGGGTGTCTGGCCTTGGGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACgg                   <  2:1151001/150‑1 (MQ=255)

                                                                                                                                    ccTTTGTGGCGGCTATTGCTGCCTCTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGa                  >  1:1219175/1‑150 (MQ=255)

                                                                                                                                           ggcggcTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGCAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAAccc           >  2:635936/1‑150 (MQ=255)

                                                                                                                                             cggcTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCGTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCg         >  1:456949/1‑150 (MQ=255)

                                                                                                                                              ggcTATTGCTGCCACTCTCCTGGGTGTATGGCCTGGTGAGGGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCgg        <  2:1767324/150‑1 (MQ=255)

                                                                                                                                                 tATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCGGTCg     <  2:1524203/150‑1 (MQ=255)

                                                                                                                                                    tgctgcCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCGGTCGTTg  >  1:38989/1‑150 (MQ=255)

CGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCGGTCGTTG  >  NC_000913/968264‑968561

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