BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,018,260	(A)_8→7	100%	coding (315/444 nt)	fliJ →	flagellar protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,018,253	0	A	.	96.2%	453.9 / 8.9	104	coding (308/444 nt)	fliJ	flagellar protein
Reads supporting (aligned to +/- strand): ref base A (1/2); new base . (43/57); total (44/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 118 total aligned reads displayed.

AGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGG  >  NC_000913/2018154‑2018348

ngcagTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCgaga                                                                                                 >  1:2990259/2‑100 (MQ=255)

agcagTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCgaga                                                                                                 >  1:1997656/1‑100 (MQ=255)

 gcagTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaa                                                                                                >  2:2415254/1‑100 (MQ=255)

 gcagTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaa                                                                                                >  1:2197459/1‑100 (MQ=255)

   ngTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaa                                                                                               >  2:3075508/2‑99 (MQ=255)

   agTTTATCCAGACGCTGGAAGAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaa                                                                                              >  1:2113074/1‑100 (MQ=255)

      ttATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaa                                                                                           >  2:1547859/1‑100 (MQ=255)

       tATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaac                                                                                          <  2:584683/100‑2 (MQ=255)

         tCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaacaa                                                                                        <  1:3299780/100‑4 (MQ=255)

         tCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaacaa                                                                                        >  2:476198/1‑97 (MQ=255)

          ccAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaacaan                                                                                       <  2:2990283/100‑5 (MQ=255)

          ccAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaacaac                                                                                       >  2:764515/1‑96 (MQ=255)

          ccAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaaca                                                                                         <  1:1796407/98‑3 (MQ=255)

              aCGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGttt                                                                                  <  2:2821922/100‑1 (MQ=255)

               cGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCATTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTg                                                                                 <  1:1064657/100‑1 (MQ=255)

                gCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGc                                                                                <  2:574858/100‑1 (MQ=255)

                gCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGc                                                                                >  2:807878/1‑100 (MQ=255)

                 cTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCa                                                                               <  2:500776/100‑1 (MQ=255)

                   ggAAAAAGCCATTAGTCAGCATGGGGAACAACTTAGTCAGAGGACGCAGAAAGTTGACATTGCCCTGAACCGTTGGCGAG‑AAAAAAACAACGTTTGCAg                                                                              <  2:292423/99‑1 (MQ=255)

                   ggAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAgg                                                                             >  2:2992989/1‑100 (MQ=255)

                   ggAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAgg                                                                             >  1:1109676/1‑100 (MQ=255)

                    gAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTCTGc                                                                                >  1:379543/1‑96 (MQ=255)

                     aaaaaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGcc                                                                           <  1:2201504/100‑1 (MQ=255)

                      aaaaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGCCATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCt                                                                          <  2:2916209/100‑1 (MQ=255)

                       aaaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTg                                                                         >  2:127217/1‑100 (MQ=255)

                         aGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGc                                                                       <  2:3214933/100‑1 (MQ=255)

                         aGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGc                                                                       <  2:1135394/100‑1 (MQ=255)

                          gCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCa                                                                      >  1:1918083/1‑100 (MQ=255)

                          gCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCa                                                                      >  1:2470813/1‑100 (MQ=255)

                          gCCATTACACAGCATCGCCAGCAACCTAATCAGCGGACGCAGAAAGTTGACATTGCCCTGAAAAAAAAAAAA‑AAAAAAAACAACGTTTGCAGGCCTGGCa                                                                      <  1:2937513/100‑1 (MQ=255)

                           ccATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAg                                                                     >  1:379233/1‑100 (MQ=255)

                            cATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGa                                                                    >  2:1139924/1‑100 (MQ=255)

                             aTTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGac                                                                   >  1:2233482/1‑100 (MQ=255)

                                aCTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACt                                                                >  2:2044030/1‑100 (MQ=255)

                                aCTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACt                                                                <  1:1553514/100‑1 (MQ=255)

                                aCTCAGCATCGCCAGCAACTGAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACt                                                                >  1:2116153/1‑100 (MQ=255)

                                   cAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCa                                                             >  1:2879594/1‑100 (MQ=255)

                                    aGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAg                                                            <  2:1510607/100‑1 (MQ=255)

                                    aGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAg                                                            <  1:1139913/100‑1 (MQ=255)

                                     gCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAgg                                                           <  1:807871/100‑1 (MQ=255)

                                     gCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAgg                                                           >  2:910721/1‑100 (MQ=255)

                                     gCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAgg                                                           >  1:2881516/1‑100 (MQ=255)

                                      cATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGa                                                          >  2:3160607/1‑100 (MQ=255)

                                       aTCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGaa                                                         <  1:362784/100‑1 (MQ=255)

                                          gCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACgg                                                      >  1:2127577/1‑100 (MQ=255)

                                          gCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACgg                                                      <  2:3056364/100‑1 (MQ=255)

                                          gCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACgg                                                      <  2:3296308/100‑1 (MQ=255)

                                           ccAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGc                                                     <  2:57598/100‑1 (MQ=255)

                                           ccAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGc                                                     <  2:2543212/100‑1 (MQ=255)

                                           ccAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGc                                                     <  2:2351227/100‑1 (MQ=255)

                                           ccAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGc                                                     <  1:2452395/100‑1 (MQ=255)

                                           ccAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGc                                                     <  1:1034910/100‑1 (MQ=255)

                                            cAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGc                                                     <  2:1918497/99‑1 (MQ=255)

                                                 aCTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCa                                               >  1:2980152/1‑100 (MQ=255)

                                                  cTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAc                                              >  1:697950/1‑100 (MQ=255)

                                                   ttAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcg                                             <  2:3296824/100‑1 (MQ=255)

                                                    tAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcgg                                            >  2:3238619/1‑100 (MQ=255)

                                                     aaTCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggc                                           <  1:2870776/100‑1 (MQ=255)

                                                     aaTCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggc                                           <  2:1991268/100‑1 (MQ=255)

                                                     aaTCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggc                                           >  2:1547535/1‑100 (MQ=255)

                                                       tCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggcgg                                         <  1:2559767/100‑1 (MQ=255)

                                                       tCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggcgg                                         >  1:3065543/1‑100 (MQ=255)

                                                       tCAGCGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggcgg                                         <  2:1863561/100‑1 (MQ=255)

                                                        cAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACggcggc                                        <  2:1850800/100‑1 (MQ=255)

                                                         aGTGGACGAAGGAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACggcggc                                        <  1:3251577/99‑1 (MQ=255)

                                                           tGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACt                                     <  2:1524950/100‑1 (MQ=255)

                                                             gACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAAGCCACGGCGGCACTGc                                   >  2:3402400/1‑100 (MQ=255)

                                                              aCGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCt                                  <  2:1109685/100‑1 (MQ=255)

                                                              aCGCAGAAAGAAGACAATGCCCCGAACAGTTGGCGA‑CAAAAAACCAACGTTTGCAGGCCTGGCAGACACTGCAGTAACGGCAATCCACGGCGGCACTGCt                                  <  2:1506428/100‑1 (MQ=255)

                                                               cGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCtt                                 <  1:1261403/100‑1 (MQ=255)

                                                                gCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTg                                <  1:1264075/100‑1 (MQ=255)

                                                                 cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                               <  2:2157666/100‑1 (MQ=255)

                                                                 cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                               <  1:2757492/100‑1 (MQ=255)

                                                                   gAAAGTTGGCAATGCCCTGACCAGCTGGGGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGGCACGGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAg                             <  1:1296934/100‑1 (MQ=255)

                                                                      aGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGaaa                          >  1:3410273/1‑100 (MQ=255)

                                                                      aGTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGaaa                          >  2:1444558/1‑100 (MQ=255)

                                                                       gTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGaaaa                         <  2:1349771/100‑1 (MQ=255)

                                                                       gTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGaaaa                         >  1:407998/1‑100 (MQ=255)

                                                                       gTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGaaaa                         >  2:2120156/1‑100 (MQ=255)

                                                                       gTTGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGaaaa                         <  2:2113091/100‑1 (MQ=255)

                                                                        ttaaCATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCATGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAAc                        <  1:139362/97‑1 (MQ=255)

                                                                        ttGACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAAc                        <  2:2006252/100‑1 (MQ=255)

                                                                          gACATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCg                      >  2:3053846/1‑100 (MQ=255)

                                                                            cATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGcc                    <  2:617995/100‑1 (MQ=255)

                                                                            cATTGCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGcc                    <  1:764509/100‑1 (MQ=255)

                                                                                gCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGa                >  2:863644/1‑100 (MQ=255)

                                                                                gCCCTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGa                >  2:1336856/1‑100 (MQ=255)

                                                                                  ccTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACTGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc              <  1:1897087/100‑1 (MQ=255)

                                                                                  ccTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc              <  1:1791346/100‑1 (MQ=255)

                                                                                  ccTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc              <  1:1169489/100‑1 (MQ=255)

                                                                                  ccTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc              <  1:521930/100‑1 (MQ=255)

                                                                                  ccTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc              >  1:883768/1‑100 (MQ=255)

                                                                                  ccTGAACAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc              <  1:662851/100‑1 (MQ=255)

                                                                                      aaCAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGaa          <  1:2536974/100‑1 (MQ=255)

                                                                                        cAGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGaaaa        <  2:3065567/100‑1 (MQ=255)

                                                                                         aGTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGaaaaa       <  1:1137749/100‑1 (MQ=255)

                                                                                          gTTGGCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAg      >  2:1763802/1‑100 (MQ=255)

                                                                                             ggCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgatg   >  1:296457/1‑100 (MQ=255)

                                                                                              gCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgatgg  >  1:1551723/1‑100 (MQ=255)

                                                                                              gCGAG‑AAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgatg   <  1:3090050/99‑1 (MQ=255)

AGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGG  >  NC_000913/2018154‑2018348

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