BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000731	25,663	Δ1 bp	100%	coding (316/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.0%	4691.5 / 56.9	1046	coding (316/333 nt)	USA300HOU_pUSA300HOUMR0031	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/17); new base . (480/545); total (480/566)
Fisher's exact test for biased strand distribution p-value = 2.63e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

Only 100 of 1223 total aligned reads displayed.

ATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAA  >  CP000731/25573‑25753

aTTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATa                                                                                            >  1:2856879/1‑91 (MQ=255)

aTTGAAAGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTAAAAATa                                                                                            <  2:1675200/91‑1 (MQ=255)

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

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

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

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

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

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

   gAAGGATCTTGCTTTTCTCAATCTCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTCATGCTGTTTAATTACTATACTACAAATaaaa                                                                                         <  2:3009390/91‑1 (MQ=255)

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

        aTCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaatat                                                                                    >  1:2788166/1‑88 (MQ=255)

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

          cTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATAc                                                                                 >  1:954573/1‑91 (MQ=255)

           ttGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACa                                                                                >  2:707692/1‑91 (MQ=255)

            tGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTGTTTCTTATTCTGTTTAATTACTATACTACAAAT‑AAAAAATATACAt                                                                               <  1:3152836/91‑1 (MQ=255)

            tGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACAt                                                                               >  2:2832342/1‑91 (MQ=255)

            tGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACAt                                                                               >  1:476783/1‑91 (MQ=255)

            tGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACAt                                                                               <  2:1224827/91‑1 (MQ=255)

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

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

                tttCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑‑AAAAATATACATCTaaa                                                                          >  1:20086/1‑91 (MQ=255)

                 ttCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTaaa                                                                          <  2:1031609/91‑1 (MQ=255)

                 ttCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTaaa                                                                          <  1:897797/91‑1 (MQ=255)

                    tcAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGt                                                                       >  2:3793036/1‑91 (MQ=255)

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

                       atatCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTaa                                                                    <  2:2757501/91‑1 (MQ=255)

                        tatCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTaaa                                                                   >  2:137454/1‑91 (MQ=255)

                        tatCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTaaa                                                                   >  1:1102090/1‑91 (MQ=255)

                        tatCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTaaa                                                                   >  1:3890763/1‑91 (MQ=255)

                            cAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTaa                                                               >  1:732678/1‑91 (MQ=255)

                               aaaTTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGa                                                            <  2:1758944/91‑1 (MQ=255)

                                  ttttCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGaaaa                                                         >  1:335886/1‑91 (MQ=255)

                                   tttCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGaaaaa                                                        <  1:3788251/91‑1 (MQ=255)

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

                                     tCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATa                                                      <  2:3392602/91‑1 (MQ=255)

                                     tCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATa                                                      >  1:2942043/1‑91 (MQ=255)

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

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

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

                                      cTTTATAATAGTAGCACTATTTCTTATGTTGTTTAATTACTATACTACAAAT‑AAAAATTATACATCTAAAGGTTAAATTAATGAAAAATAg                                                     >  1:1862152/1‑91 (MQ=255)

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

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

                                            aaTTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTag                                               >  1:835235/1‑91 (MQ=255)

                                             aTTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTaga                                              <  1:3513768/91‑1 (MQ=255)

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

                                                gTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGaa                                           >  2:452936/1‑91 (MQ=255)

                                                gTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGaa                                           >  2:51730/1‑91 (MQ=255)

                                                gTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGaa                                           >  1:25613/1‑91 (MQ=255)

                                                gTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGaa                                           >  1:3956844/1‑91 (MQ=255)

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

                                                 tAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAat                                          >  1:3606518/1‑91 (MQ=255)

                                                 tAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAat                                          >  1:2610303/1‑91 (MQ=255)

                                                  aGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAAATAATGAAAAATAGAATTAGAGAata                                         >  1:3664227/1‑91 (MQ=255)

                                                  aGCACTATTTCTGATGCTGTTTAATTACTATACGACAAAT‑AAAAAATATACATCTAAAGGTTAAGTAAATGAAAAATAGAATTAGAGAata                                         <  2:3283803/91‑1 (MQ=255)

                                                     aCTCTTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAg                                      >  2:3501785/1‑91 (MQ=255)

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

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

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

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

                                                        aTTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGaaa                                   <  2:1327736/91‑1 (MQ=255)

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

                                                          ttCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGaaaca                                 <  2:2932875/91‑3 (MQ=255)

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

                                                           tCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACa                                >  1:1225909/1‑91 (MQ=255)

                                                            cTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaa                               >  1:1931569/1‑91 (MQ=255)

                                                             ttATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGTAAACaaa                              >  1:1967214/1‑91 (MQ=255)

                                                             ttATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaa                              <  1:3428240/91‑1 (MQ=255)

                                                              tATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaaa                             <  1:168663/91‑1 (MQ=255)

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

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

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

                                                                tGCTGTTAAATTACGATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAg                           <  2:2380946/91‑1 (MQ=255)

                                                                aacTGATTAACTACTACAGCAAAAAT‑AAAAAATGTACCTCTAAAGGTTAAATCAAAGAATTATACTATTAGAGAATATACAAAACAAAAAg                           <  2:3662526/89‑1 (MQ=255)

                                                                 gCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAgg                          <  2:3071881/91‑1 (MQ=255)

                                                                  cTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGa                         <  2:75882/91‑1 (MQ=255)

                                                                  cTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGa                         <  1:1009077/91‑1 (MQ=255)

                                                                  cTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTACAGAATATAGAAAACAAGAAGGa                         >  2:3196169/1‑91 (MQ=255)

                                                                    gTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAAt                       >  2:3426075/1‑91 (MQ=255)

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

                                                                      ttAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAAttt                     >  1:265955/1‑91 (MQ=255)

                                                                        aaTTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAAAAAAAAGGAATTtct                   >  2:1827778/1‑91 (MQ=255)

                                                                         aTTACTATACTAAAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTtctc                  <  2:2524633/91‑1 (MQ=255)

                                                                          ttACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCa                 <  1:3563210/91‑1 (MQ=255)

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

                                                                          ttACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCa                 <  1:2997408/91‑1 (MQ=255)

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

                                                                           tactatactaCAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCaa                >  2:596892/1‑91 (MQ=255)

                                                                            gctatactaCAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATAAAGAAAACAAAAAGGAATTTCTCAAg               <  2:1888114/90‑1 (MQ=255)

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

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

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

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

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

                                                                                 actacAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAAAATAAAAAAAAAAAAGGAATTTCTCAAGAAACt          <  2:632238/91‑1 (MQ=255)

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

                                                                                  ctacAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACtt         <  1:2307891/91‑1 (MQ=255)

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

                                                                                   tacAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACttt        <  1:220182/91‑1 (MQ=255)

                                                                                   aacAAAT‑AAAAAATATACATATAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATAAAAAAAAAAAAAAGGAATTTCTCAAGAAACttt        <  2:3599826/90‑1 (MQ=255)

                                                                                    acAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTa       <  1:2406388/91‑1 (MQ=255)

ATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAA  >  CP000731/25573‑25753

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