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
*CP00073125,6630A.98.4% 272.6 / ‑3.4 61coding (316/333 nt)USA300HOU_pUSA300HOUMR0031hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (35/25);  minor base C (0/1);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTATGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGAAAACAATAAAT  >  CP000731/25526‑25803
                                                                                                                                         |                                                                                                                                            
gtatGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaa                                                                                                                                           <  1:317681/141‑1 (MQ=255)
  atGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaa                                                                                                                                         <  1:63160/141‑1 (MQ=255)
     ttATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAATGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaatat                                                                                                                                      <  2:220102/141‑4 (MQ=255)
      tATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATaaaaaatata                                                                                                                                     <  1:344656/141‑5 (MQ=255)
       aTTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTAGAAAT‑AAAAAATATAc                                                                                                                                   >  1:404673/1‑141 (MQ=255)
          tACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTCTGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATc                                                                                                                                >  2:122766/1‑141 (MQ=255)
            cacTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTa                                                                                                                              <  1:180887/141‑1 (MQ=255)
            aacTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAAATGGTAT‑CTGTCAAACAAAAAATATACATcga                                                                                                                              <  2:413661/140‑3 (MQ=255)
             acTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATCCATCTaa                                                                                                                             >  1:102057/1‑141 (MQ=255)
             acTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTaa                                                                                                                             >  1:124898/1‑141 (MQ=255)
             acTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTaa                                                                                                                             >  2:87613/1‑141 (MQ=255)
                 gTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAgat                                                                                                                         >  2:70958/1‑139 (MQ=255)
                 gTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGt                                                                                                                         <  1:64407/141‑1 (MQ=255)
                  tGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGtt                                                                                                                        <  1:42789/141‑1 (MQ=255)
                     ttAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTaaa                                                                                                                     >  1:231595/1‑141 (MQ=255)
                      tAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTCTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAAt                                                                                                                    >  2:139999/1‑141 (MQ=255)
                       aaaTTTACATGTTTGTCACAAAGTATTGAAGGATCTTGGTTTTCTCAATATCTAGAAATTTTATTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAAtt                                                                                                                   <  2:55047/141‑1 (MQ=255)
                         aTTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATATATCTAAAGGTTAAATTaa                                                                                                                 >  2:84984/1‑141 (MQ=255)
                          tttACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAAt                                                                                                                >  1:338789/1‑141 (MQ=255)
                              cAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCGAAAGGTTCAATTAATGaaa                                                                                                            >  1:72449/1‑141 (MQ=255)
                              cAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATAAATCTAAAGGTTAAATTAATGaaa                                                                                                            >  1:223173/1‑141 (MQ=255)
                                   ttGTCCCCAAGTATTGAAGGATCTTGGTTTTCTCAATCTCCAGACATTTTCTTTATAATTGTCGCACTATTTCTTATGCTCTTTAATTACTCTACTACAACT‑AAAACATATACATCTACAACTTAAATTAACaaaaaacaa                                                                                                       >  2:103453/1‑138 (MQ=255)
                                    tGTAACAAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACAATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGa                                                                                                      <  2:98097/141‑1 (MQ=255)
                                       cACCAAGTATTGAAGGATCTTGGTTTTCTCAATATGCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAAtt                                                                                                   >  1:212096/1‑141 (MQ=255)
                                         acAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTAATGCTGTTTAATTAGGAAACAACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTag                                                                                                 >  1:399703/2‑141 (MQ=255)
                                          cAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTaga                                                                                                >  1:131472/1‑141 (MQ=255)
                                           aaGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTagag                                                                                               >  1:144349/1‑141 (MQ=255)
                                            aGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATGACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAgaga                                                                                              <  1:61597/141‑1 (MQ=255)
                                               aTTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTTTTAAGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAAGTTAAATTAATGAAAAAAGGAAATGTAGAata                                                                                           >  2:251605/1‑141 (MQ=255)
                                                  gAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAg                                                                                        >  1:332221/1‑141 (MQ=255)
                                                  gAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTACCACTATTCTGCCCATTATATAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAg                                                                                        <  1:97601/141‑1 (MQ=255)
                                                   aaGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGa                                                                                       <  1:153899/141‑1 (MQ=255)
                                                        tCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACa                                                                                  <  1:359536/141‑1 (MQ=255)
                                                        tCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACa                                                                                  <  2:291910/141‑1 (MQ=255)
                                                         cTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaa                                                                                 >  1:262015/1‑141 (MQ=255)
                                                          ttGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaa                                                                                <  1:161278/141‑1 (MQ=255)
                                                           tGGTTTTCTCAATATCCTGCTATTGGCTAAATAATTGTAGAACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACaaaa                                                                               <  2:363757/141‑1 (MQ=255)
                                                             gTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAATAAAAAg                                                                             >  2:290621/1‑141 (MQ=255)
                                                              ttttCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAga                                                                            >  1:287176/1‑140 (MQ=255)
                                                              ttttATTAATATCAAAAAATTTTCATAATAATTGAAGAATAAATTATTATTATGTTAAATTACTATAAAACAAAT‑AAAAAATATACATCTAAAAGTTAAATTAATGAAAAATAGAATTAGAGAATAAAAAAAAAAAAAAgg                                                                            <  2:85977/141‑1 (MQ=255)
                                                                ttCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATCTACTCTTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGaa                                                                          >  1:290107/1‑141 (MQ=255)
                                                                    cAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTtc                                                                      <  1:20421/141‑1 (MQ=255)
                                                                       tataCAGAAAATTTATAAATAATTGTAGAACTATTTCTTATGCTATTAAATTAGTAAACTACAAAT‑AAAAAATATAAATCAAAAGGTTAAATTAATGAAAAATAGAATTAGAGAAAAAAAAAAAAAAAAAGGAATTTCTCa                                                                   <  2:172014/141‑1 (MQ=255)
                                                                        atCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCaa                                                                  <  2:124898/141‑1 (MQ=255)
                                                                         tCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAg                                                                 <  1:245164/141‑1 (MQ=255)
                                                                               aaTTTTATTAATAATTATAACACAATTAATTATGCTGTTTAATTAATATACTAAAAAT‑AAAAAATATACATCTAAAAGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACtt                                                           <  2:163213/141‑1 (MQ=255)
                                                                                    tCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAAGGTTAAga                                                      >  1:50114/1‑140 (MQ=255)
                                                                                     cTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAACACAAAAAAGAATATCTAAAAAAAATATAAGTa                                                     >  2:252769/1‑141 (MQ=255)
                                                                                     cTTTATAATTATAGCACTATTTCTAATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTa                                                     >  1:101692/1‑141 (MQ=255)
                                                                                              tGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGACTTAGAGAATATAGAAAACAAAAAGGATGTTCTCAAGAAACTTTAAGTAAAGAATTaa                                            >  1:265970/1‑141 (MQ=255)
                                                                                              tGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTaa                                            <  1:139999/141‑1 (MQ=255)
                                                                                                  gCACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTCAAAgt                                        >  2:236150/1‑141 (MQ=255)
                                                                                                   cACTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAAAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAgtg                                       >  2:171544/1‑141 (MQ=255)
                                                                                                    aCTATTTCTTATGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACAAAACAAGCTAAATTTAATGAAAAATAGAATTAGACAATATAGAAAACAAAAAGGACATTCTAAAGAAACTTTAAGTAACGAATTACAAatgt                                      >  2:172532/1‑141 (MQ=255)
                                                                                                         ttCTTATGCTGTTTAATTAATATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAgg                                 <  1:282187/141‑1 (MQ=255)
                                                                                                          tCTTATGCTATTTAAATAATATAAAACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGc                                <  2:292384/141‑1 (MQ=255)
                                                                                                              aTGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAAc                            <  1:87613/141‑1 (MQ=255)
                                                                                                               tGCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTCAAGGTTAAATTACTGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTATAGGCAAACt                           >  2:125997/1‑141 (MQ=255)
                                                                                                                gCTGTTTAATTACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTa                          >  1:3347/1‑141 (MQ=255)
                                                                                                                  tagttAATTACTATACTACAAAT‑AAAAAATATAAATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTGCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTAtt                        <  1:248190/138‑1 (MQ=255)
                                                                                                                       aaatactatactaAAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGc                   <  1:134127/138‑1 (MQ=255)
                                                                                                                         ttACTATACTACAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCaa                 >  2:198662/1‑141 (MQ=255)
                                                                                                                          tactatactaCAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAAt                >  1:144321/1‑141 (MQ=255)
                                                                                                                                   acAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGAAAACaa       <  2:63634/141‑1 (MQ=255)
                                                                                                                                    cAAAT‑AAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGAGTCTAGGCAAACTATTAATGCAATAGAcaacaaa      >  2:65548/1‑140 (MQ=255)
                                                                                                                                         aaaaaaaTATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGAAAACaataaat  <  1:122766/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                            
GTATGTTATTTACACTAGTGCTTAAATTTACAGGTTTGTCACCAAGTATTGAAGGATCTTGGTTTTCTCAATATCCAGAAATTTTCTTTATAATTGTAGCACTATTTCTTATGCTGTTTAATTACTATACTACAAATAAAAAAATATACATCTAAAGGTTAAATTAATGAAAAATAGAATTAGAGAATATAGAAAACAAAAAGGAATTTCTCAAGAAACTTTAAGTAAAGAATTAAAAGTGTCTAGGCAAACTATTAATGCAATAGAAAACAATAAAT  >  CP000731/25526‑25803

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