Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,987,676 A→G 5.5% intergenic (+258/+38) USA300HOU_1847 → / ← USA300HOU_1848 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,987,6760AG 5.5% 205.3 / 2.2 72intergenic (+258/+38)USA300HOU_1847/USA300HOU_1848hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (37/30);  new base G (2/2);  total (40/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

GCCCCAACACTGAAACTAGCAGAAAGTCAGCTTCTATGAATATAATAAAAAAGCTAGGTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATCATTTTGATCAGAGATTTGTTCCATTTCTTTACCTAATTCTTTTAAATCTTCAAAATCCGTTACC  >  CP000730/1987543‑1987813
                                                                                                                                     |                                                                                                                                         
gCCCCAACACTGAAACTAGCAGAAAGTCAGCTTCTATGAATATAATAAAAAAGCTAGGTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGata                                                                                                                                    >  1:467283/1‑141 (MQ=255)
       cacTGAAACGAGCCGAAAGTCAGCTTCTATGAATATAATAAAAAAGCTAGGTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATGTGATATATGTTCTATACTATACATTATTTACATGATAAATAACt                                                                                                                             <  2:413072/141‑1 (MQ=255)
       aacTGAAACTAGCAGAAAGTCAGCTTCTATGAATATAATAAAAAAGCTAGGTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACt                                                                                                                             <  2:362704/140‑1 (MQ=255)
          tGAAACTAGCAGAAAGTCAGCTTCTATGAATATAATAAAAAAGCTAGGTAACACAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTGCATGATAAATAACTGaa                                                                                                                          >  2:253964/1‑141 (MQ=255)
          tGAAACTAGCAGAAAGTCAGCTTCTATGAATATAATAAAAAAGCTAGGTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTGCATGATAAATAACTGaa                                                                                                                          >  2:261148/1‑141 (MQ=255)
           gAAACTAGCAGAAAGTCAGCTTCTATGAATATAATAAAAAAGCTAGGCAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAat                                                                                                                         <  1:67514/141‑1 (MQ=255)
             aaCTAGCAGAAATGCAGCTTCTTTGAATAGAAAACCAAAGCTAGGTAACAACATGTTACCTAACTTCATATTCAAGTTAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAatat                                                                                                                       <  1:743267/141‑1 (MQ=255)
             aaCTAGCAGAAAGTCAGCTCCTATGAATATAATAAAAAAGCTATGTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAatat                                                                                                                       <  1:645349/141‑1 (MQ=255)
             aaCTAGCAGAAAGGCAGCTTCTATGAATATAATAAAAAAGCTAGGTAACAAAATGATACCTAACTTCATATTCAAGATAATCAAGCCTATTTGATAAATGTTCTATACTATACATTAGTTACATGATAAATAACTGAatat                                                                                                                       >  2:220856/1‑141 (MQ=255)
               cTAGCAGAAAGTCAGCTTCTATGAATATAATAAAAAAGCTAGGTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTa                                                                                                                     <  2:232500/141‑1 (MQ=255)
                 agaaGAAAGCCAACTTCTATAAATATAATAAAAAATCTAGGTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTaca                                                                                                                   <  2:457171/138‑1 (MQ=255)
                 aGCAGAAAGTCAGCTTCTATGAATATAATAAAAAAGCTAGGTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTaca                                                                                                                   <  2:217220/141‑1 (MQ=255)
                 aGCAGAAAGTCAGCTTCTATGAATATAATAAAAAAGCTAGGTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTaca                                                                                                                   >  1:425397/1‑141 (MQ=255)
                     gAAAGTCAGCTTCTATGAATATAATAAAAAAGCTAGGTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAAt                                                                                                               >  1:406485/1‑141 (MQ=255)
                     gAAAGTCAGCTTCTATGAATATAATAAAAAAGCTAGGTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAAt                                                                                                               >  2:220916/1‑141 (MQ=255)
                      aaaGTCAGCTTCTATGAATATAATAAAAAAGCTAGGTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTCCATGATAAATAACTGAATATTACACAAtt                                                                                                              >  2:122612/1‑141 (MQ=255)
                        aGTCAGCTTCTATGAATATAATAAAAAAGCTAGGTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATtat                                                                                                            <  2:335081/141‑1 (MQ=255)
                          tCAGCTTCTATGAATATAATAAAAAAGCTAGGTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTataa                                                                                                          <  1:190916/141‑1 (MQ=255)
                            aGCTTCTATGAATATAATAAAAAAGCTAGGTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTataata                                                                                                        >  1:72413/1‑141 (MQ=255)
                             gCTTCTATGAATATAATAAAAAAGCTAGGTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATAc                                                                                                       >  1:695929/1‑141 (MQ=255)
                                     gAATATAATAAAAAAGCTAGGTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTactga                                                                                               >  2:239512/1‑141 (MQ=255)
                                       atatAATAAAAAAGCTAGGTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTactgact                                                                                             <  1:416689/141‑1 (MQ=255)
                                         ataataAAAAAGCTAGGTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGt                                                                                           <  2:209249/141‑1 (MQ=255)
                                          taataaAAAAGCTAGGTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGtc                                                                                          >  1:279405/1‑141 (MQ=255)
                                          aaataaAAAAGCTAGGTAACAAAAGGATACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGtc                                                                                          <  2:652189/140‑1 (MQ=255)
                                             taaAAAAGCTAGGTAACAAAATGCTACCTAACTTCATATTCAAGAAAATAAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATACCTGAATATTACACAATTATAATACTTTACTGACTGtcttc                                                                                       >  2:180689/1‑141 (MQ=255)
                                              aaaaaaGCTAGCTAAAACATTGCTACCTAATTACATATTCAAGATAATCAAGCCTAGTTGATATATGTTCTAAACTATACATTATTTACATGATAAATAACTGAATATTACACATTTATAAGACTTTACTGACTGtcttct                                                                                      <  1:551237/141‑1 (MQ=255)
                                                gaaaGCTAGGTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACCATTATAATACTTTTCTGACTGTCttcttc                                                                                    >  2:433938/2‑141 (MQ=255)
                                                  aaGCTAGGTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAg                                                                                  <  2:725937/141‑1 (MQ=255)
                                                   aGCTAGGTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATAAATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGAATCTCTCCTTCaag                                                                                 >  2:105643/1‑139 (MQ=255)
                                                    gCTAGGTCAACAAATGCTACCTACATTTTCAATGGAGATATGCAATTTTACTTGATTTAAGTTCTATACGATACATGATTTACTTGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGat                                                                                >  2:153529/1‑140 (MQ=255)
                                                    aatAGGAAACAAAATGCTTTCTAACTTCAGATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTGTACTGACTGTCTTCTTCAgta                                                                                <  2:532657/139‑3 (MQ=255)
                                                       aGGTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTAAGAAttt                                                                             >  2:741886/1‑140 (MQ=255)
                                                         gTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATCATACTTTACTGACTGTCTTCTTCAGAAttctt                                                                           >  1:636384/1‑141 (MQ=255)
                                                           aaCAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAAttctttt                                                                         >  1:707682/1‑141 (MQ=255)
                                                           aaCAAAATGCTACCTAACTGCATATGCAAGATAATAACTCCTATGTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTATACAATTATAATACTTTACTGCTTGTCTTCTTCAGAAttgtttt                                                                         >  1:424210/1‑141 (MQ=255)
                                                            aCAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACGCAATTATAATACTTTACTGACTGTCTTCTTCAGAAttcttttc                                                                        <  2:35160/141‑1 (MQ=255)
                                                             ccaaaTTCTAACTAACTTCATATTCAAAATAATCAATCCTATGTCATATATGTTCTATACTATACATTATTTACCTTATAACTAACGGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAAttcttttct                                                                       <  1:242294/139‑1 (MQ=255)
                                                               aaaTGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATACTTTTCTTg                                                                     >  1:214333/1‑141 (MQ=255)
                                                                  tGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACACTTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTTAta                                                                  >  1:271726/1‑140 (MQ=255)
                                                                   gCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATCa                                                                 <  1:220916/141‑1 (MQ=255)
                                                                        cTAACTTCATATTCAAGATGATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATCATTTTg                                                            <  1:494856/141‑1 (MQ=255)
                                                                         tAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTAACTGTCGTCTTCAGAATTCTTTTCTTGAGCATTTTGa                                                           >  1:270788/1‑141 (MQ=255)
                                                                           aCTTTATATTCAAGATAATGTATTCTAATTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTGCTTCAGAATTCTTTTCTTGATCATTTTGATc                                                         <  1:572873/141‑1 (MQ=255)
                                                                              tCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATTATTTTGATCagg                                                      >  2:90007/1‑140 (MQ=255)
                                                                              tCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATCATTTTGATCaga                                                      <  1:673138/141‑1 (MQ=255)
                                                                                atatTCAAGATAATCCATACATTATAAAATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATCATTTTGATCAgaga                                                    >  1:571315/1‑141 (MQ=255)
                                                                                atatTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATCATTTTGATCAgaga                                                    >  1:491137/1‑141 (MQ=255)
                                                                                     cAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATCATTTTGATCAGAGATTTGt                                               >  1:757491/1‑141 (MQ=255)
                                                                                     cAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATCATTTTGATCAGAGATTTGt                                               >  2:286241/1‑141 (MQ=255)
                                                                                               aatcCTATTTGATATATGTTCTATACGATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATCATTTTGATCAGAGATTTGTTCCATTTCtt                                     >  2:90341/1‑141 (MQ=255)
                                                                                                atcCTACTCGATATATGCTCTATCCTATCCATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACCGACTGTCTTCTTCAGAACTCTTTTCTTGATCATCTTGATCAGAGATTTGTTCCATTTCtta                                    <  2:71225/141‑2 (MQ=255)
                                                                                                   cTATTTGATATATGTTTTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACCGTCTTCGTCGGAATTCTTTTCTTGATCATTTTGATCAGAGATTTGTTCCATTTCTTTAcc                                 <  1:147459/141‑1 (MQ=255)
                                                                                                    tATTTGATATAGGTTGTATAATATACATGATTTACATGATAAATAACTGAAGATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATCATTTTGATCAGAGATTTGTTCCATTTCTTTACCt                                <  2:636384/141‑1 (MQ=255)
                                                                                                     aTTTGATAAATGTTTTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATCATTTTGATCAGAGATTTGTTCCATTTCTTTACCTa                               <  1:24004/141‑1 (MQ=255)
                                                                                                      tttGATATATGTTCTATACTATACATTATTTGCATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTTTTCAGAATTCTTTTCTTGATCATTTTGATCAGAGATTTGTTCCATTTCTTTACCTaa                              <  1:253964/141‑1 (MQ=255)
                                                                                                      tttGATATATGTTCTATACTATACATTATTTGCATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATCATTTTGATCAGAGATTTGTTCCATTTCTTTACCTaa                              <  1:261148/141‑1 (MQ=255)
                                                                                                      tttGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATCATTTTGATCAGAGATTTGTTCCATTTCTTTACCTaa                              <  1:168796/141‑1 (MQ=255)
                                                                                                          atatatGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATCATTTTGATCAGAGATTTGTTCCATTTCTTTACCTAATTCt                          >  1:552878/1‑141 (MQ=255)
                                                                                                             tatGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATCATTTTGATCAGAGATTTGTTCCATTTCTTTACCTAATTCtttt                       >  2:545982/1‑141 (MQ=255)
                                                                                                               tGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATCATTTTGATCAGAGATTTGTTCCATTTCTTTACCTAATTCTTTTaa                     >  1:147502/1‑141 (MQ=255)
                                                                                                               tGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATCATTTTGATCAGAGATTTGTTCCATTTCTTTACCTAATTCTTTTaa                     >  1:411385/1‑141 (MQ=255)
                                                                                                                 ttCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATTATACTTTACTGACTGTCTTCTTCAGAATTCTTATCTTGATCATTTTGATCAGAGATTTGTTCCATTTCTTTACCTAATTCTCTTAAAt                   <  2:488277/141‑1 (MQ=255)
                                                                                                                   ctatactataCATTATTTAGATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATCATTTTGATCAGAGATTTGTTCCATTTCTTTACCTAATTTTTTTAAATCt                 >  1:411330/1‑141 (MQ=255)
                                                                                                                     atactatacATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATCATTTTGATCAGAGATTTGTTCCATTTCTTTACCTAATTCTTTTAAATCTTc               <  2:72413/141‑1 (MQ=255)
                                                                                                                       actatacATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATCATTTTGATCAGAGATTTGTTCCATTTCTTTACCTAATTCTTTTAAATCTTCaa             >  2:326925/1‑141 (MQ=255)
                                                                                                                       actatacATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATCATTTTGATCAGAGATTTGTTCCATTTCTTTACCTAATTCTTTTAAATCTTCaa             >  1:313058/1‑141 (MQ=255)
                                                                                                                        ctatacATTATTTACTTTATAAAGAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATCATTTTGATCAGAGATTTGTTCCATTTCTTTACCTAATTCTTTTAAATCTTCaaa            <  1:56899/141‑1 (MQ=255)
                                                                                                                               ttatcTACATGATAAATAACTGAATATGACACATTTATAATACTTTTCTGACTGTTTTCTTCAGAATTCTTTTCTTGATCATTTTGATCAGAGATTTGTTCCAATTCTTTACCTAATTCTTTTAAATCTTCAAAATCCGtt     >  1:546262/1‑141 (MQ=255)
                                                                                                                                tattTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATCATTTTGATCAGAGATTTGTTCCATTTCTTTACCTAATTGTTTTAAATCTTCAAAATCCGttc    >  2:60350/1‑140 (MQ=255)
                                                                                                                                tattTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATCATTTTGATCAGAGATTTGTTCCATTTCTTTACCTAATTCTTTTAAATCTTCAAAATCCGTTa    >  2:676928/1‑141 (MQ=255)
                                                                                                                                  tttACATGATAAATAACTGAATATTACACAATTATAATGCTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATCATTTTAATCAGAGATTTGTTCCATTTCTTTACCTAATTCTTTTAAATCTACAAAATCCGTTAcc  >  2:410940/1‑141 (MQ=255)
                                                                                                                                     |                                                                                                                                         
GCCCCAACACTGAAACTAGCAGAAAGTCAGCTTCTATGAATATAATAAAAAAGCTAGGTAACAAAATGCTACCTAACTTCATATTCAAGATAATCAATCCTATTTGATATATGTTCTATACTATACATTATTTACATGATAAATAACTGAATATTACACAATTATAATACTTTACTGACTGTCTTCTTCAGAATTCTTTTCTTGATCATTTTGATCAGAGATTTGTTCCATTTCTTTACCTAATTCTTTTAAATCTTCAAAATCCGTTACC  >  CP000730/1987543‑1987813

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: