BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000731	25,669	(A)_7→6	coding (322/333 nt)	USA300HOU_pUSA300HOUMR0031 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000731	25,663	0	A	.	98.8%	6092.0 / 53.8	1244	coding (316/333 nt)	USA300HOU_pUSA300HOUMR0031	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/15); new base . (678/549); total (678/566)
Fisher's exact test for biased strand distribution p-value = 6.50e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1429 total aligned reads displayed.

ATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAA  >  CP000731/25573‑25753

aTTGTAGTATCTTGGTTTTCTCTATTTCCAGAAATTTTCTTTATAATTTTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATa                                                                                            <  2:3839794/91‑1 (MQ=255)

aTTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATa                                                                                            <  2:4129722/91‑1 (MQ=255)

 ttGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaa                                                                                           >  2:2518025/1‑91 (MQ=255)

 ttGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaa                                                                                           >  1:143458/1‑91 (MQ=255)

 ttGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaa                                                                                           <  1:2771354/91‑1 (MQ=255)

  tGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaa                                                                                          >  2:43359/1‑91 (MQ=255)

  tGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaa                                                                                          >  1:3934975/1‑91 (MQ=255)

   gAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaa                                                                                         >  2:5979951/1‑91 (MQ=255)

   gAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaa                                                                                         >  2:2994200/1‑91 (MQ=255)

    aaGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaa                                                                                        <  1:639984/91‑1 (MQ=255)

    aaGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaa                                                                                        >  1:3107790/1‑91 (MQ=255)

       gATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaata                                                                                     <  2:4504264/91‑3 (MQ=255)

         tCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaatata                                                                                   >  2:2801624/1‑87 (MQ=255)

         tCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaatata                                                                                   <  2:2698806/91‑5 (MQ=255)

           ttGGTTTTCTCAATATCCAGAAATTTTCTCTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACa                                                                                >  2:107882/1‑91 (MQ=255)

             ggTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATc                                                                              >  1:1914411/1‑91 (MQ=255)

             ggTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATc                                                                              >  2:3558749/1‑91 (MQ=255)

             ggTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATc                                                                              >  1:1252781/1‑91 (MQ=255)

             ggTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATc                                                                              >  2:1526637/1‑91 (MQ=255)

               ttttCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTa                                                                            >  1:957324/1‑91 (MQ=255)

               ttttCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTa                                                                            >  2:4804183/1‑91 (MQ=255)

                tttCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTaa                                                                           <  1:4227039/91‑1 (MQ=255)

                tttCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTaa                                                                           >  2:5349960/1‑91 (MQ=255)

                  tctcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAg                                                                         >  1:3554392/1‑91 (MQ=255)

                  tctcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAg                                                                         >  1:2874208/1‑91 (MQ=255)

                  tctcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAg                                                                         <  1:1491475/91‑1 (MQ=255)

                   ctcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAgg                                                                        <  1:5910575/91‑1 (MQ=255)

                   ctcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAgg                                                                        <  2:2120258/91‑1 (MQ=255)

                   ctcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAgg                                                                        <  2:57796/91‑1 (MQ=255)

                   ctcAATATCCAGAAATTTTATTTATAATTGTAGCACTATTTCTTAGGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAgg                                                                        <  1:4012124/91‑1 (MQ=255)

                    tcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGt                                                                       <  2:4910917/91‑1 (MQ=255)

                     cAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGtt                                                                      >  1:4455879/1‑91 (MQ=255)

                        tatCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTAAAAAT‑AAAAAATATACATCTAAAGGTTaaa                                                                   <  2:2590986/91‑1 (MQ=255)

                          tCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAAtt                                                                 <  2:3169321/91‑1 (MQ=255)

                           ccAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTa                                                                >  2:6093391/1‑91 (MQ=255)

                            cAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAGAGGTTAAATTaa                                                               >  2:4209248/1‑91 (MQ=255)

                             aGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAAt                                                              <  1:1361566/91‑1 (MQ=255)

                              gAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTATTATACTACAAAT‑CAAAAATATACATCTAAAGGTTAAATTAATg                                                             >  1:6303215/1‑91 (MQ=255)

                               aaaTTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGa                                                            >  1:1142723/1‑91 (MQ=255)

                               aaaTTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCCAAAGGTTAAATTAATGa                                                            >  1:2648912/1‑91 (MQ=255)

                               aaaTTTTCTTTATAATTGAAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGa                                                            <  1:4970430/91‑1 (MQ=255)

                                 aTTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAATTATACATCTAAAGGTTAAAATAATGaaa                                                          <  2:6298822/91‑1 (MQ=255)

                                    ttCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAAt                                                       <  1:3246913/91‑1 (MQ=255)

                                     tCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATATAAAGGTTAAATTAATGAAAAATa                                                      >  2:5667615/1‑91 (MQ=255)

                                      cTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAg                                                     >  1:1568973/1‑91 (MQ=255)

                                      cTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAg                                                     >  2:6184088/1‑91 (MQ=255)

                                      cTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAg                                                     >  1:2352968/1‑91 (MQ=255)

                                      cTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAg                                                     <  1:226681/91‑1 (MQ=255)

                                      cTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAg                                                     >  1:5092235/1‑91 (MQ=255)

                                      cTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCAAAAGGTTAAATTAATGAAAAATAg                                                     >  1:5192487/1‑91 (MQ=255)

                                       tttATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGa                                                    <  1:6550643/91‑1 (MQ=255)

                                        ttATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGaa                                                   <  2:5022458/91‑1 (MQ=255)

                                        ttATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGaa                                                   >  2:1805183/1‑91 (MQ=255)

                                         tataATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAAt                                                  <  1:3448741/91‑1 (MQ=255)

                                          ataATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAAtt                                                 >  1:703564/1‑91 (MQ=255)

                                           taATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAAt                                                       >  2:1973828/1‑84 (MQ=255)

                                              ttGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTagag                                             <  2:3107790/91‑1 (MQ=255)

                                              ttGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTagag                                             <  2:860295/91‑1 (MQ=255)

                                                gTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGaa                                           <  1:2085658/91‑1 (MQ=255)

                                                gTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGATTAAATTAATGATAAATAGTATTAGAGaa                                           <  1:6056789/91‑1 (MQ=255)

                                                gTAGCAATATTTCTTATGCTGTTTAATTACTATACTAAAAAT‑AAAAAATATACATTTAAAGGAAAAATAAAAAAAAAATAGAATTAGAGaa                                           <  2:5551431/91‑1 (MQ=255)

                                                 tAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGaaa                                          >  2:5138881/1‑90 (MQ=255)

                                                  aGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAata                                         <  2:2059455/91‑1 (MQ=255)

                                                   gCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAatat                                        <  2:3885176/91‑1 (MQ=255)

                                                     aCTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAg                                      >  1:1072276/1‑91 (MQ=255)

                                                     aCTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAg                                      >  1:5769579/1‑91 (MQ=255)

                                                     aCTAGTTCTGATGCTGTGTAATTACTATACTACAAAT‑AAAAAATATACATCTGAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAg                                      >  1:4300719/1‑91 (MQ=255)

                                                       tATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGaa                                    <  1:1001115/91‑1 (MQ=255)

                                                        aTTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGaaa                                   >  1:6201663/1‑91 (MQ=255)

                                                         tttCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGaaaa                                  <  1:2495661/91‑1 (MQ=255)

                                                         tttCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGaaaa                                  <  2:3405098/91‑1 (MQ=255)

                                                         tttCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGaaaa                                  >  1:1716244/1‑91 (MQ=255)

                                                          ttCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAAc                                 >  2:3272060/1‑91 (MQ=255)

                                                          ttCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAAc                                 >  2:1204185/1‑91 (MQ=255)

                                                             ttATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaa                              >  1:861955/1‑91 (MQ=255)

                                                              tATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaaa                             >  1:4833563/1‑91 (MQ=255)

                                                               aTGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaaaa                            >  2:546551/1‑91 (MQ=255)

                                                                 gCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAgg                          >  2:5898262/1‑91 (MQ=255)

                                                                  cTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGa                         >  1:260316/1‑91 (MQ=255)

                                                                  cTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGa                         >  2:1642101/1‑91 (MQ=255)

                                                                   tGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAAAAAAAAAAAAAAAAAGGaa                        <  2:1381634/91‑1 (MQ=255)

                                                                      ttAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAAttt                     <  1:4717957/91‑1 (MQ=255)

                                                                         aTTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTtctc                  <  1:56185/91‑1 (MQ=255)

                                                                          ttACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCa                 >  2:2386560/1‑91 (MQ=255)

                                                                          ttACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAAAATAAAAAACAAAAAGGAATTTCTCa                 <  2:660942/91‑1 (MQ=255)

                                                                           tactatactaCAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCaa                <  1:4141589/91‑1 (MQ=255)

                                                                            actatactaCAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAg               <  2:4708147/91‑1 (MQ=255)

                                                                            actatactaCAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAg               >  1:5304354/1‑91 (MQ=255)

                                                                            actatactaCAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAg               >  2:3702610/1‑91 (MQ=255)

                                                                            actatactaCAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAg               <  2:2891379/91‑1 (MQ=255)

                                                                            actatactaCAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAAAAAAAAAAAAAAAAAGGAATTTCTCAAg               <  2:461424/91‑1 (MQ=255)

                                                                             ctatactaCAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGa              <  1:3735528/91‑1 (MQ=255)

                                                                              tatactaCAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGaa             >  1:5457033/1‑91 (MQ=255)

                                                                                tactacAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAAc           >  2:2262072/1‑91 (MQ=255)

                                                                                tactacAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAAc           <  1:2991776/91‑1 (MQ=255)

                                                                                tactacAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAAc           >  1:3625943/1‑91 (MQ=255)

                                                                                 actacAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACt          >  2:751875/1‑91 (MQ=255)

                                                                                  ctacAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACtt         >  2:4035338/1‑91 (MQ=255)

                                                                                   tacAAAT‑AAAAAATATACATCTAGAGGTTAAATTAGTGAAAGATAGAATTAGAGAATATAGAAATCAAAAAGGAATTTCTCAAGAAACttt        <  1:4547748/91‑1 (MQ=255)

                                                                                     cAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTaa      >  2:6447444/1‑91 (MQ=255)

ATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAA  >  CP000731/25573‑25753

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: `‑`