Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,120,712 G→C 7.5% R67P (CGC→CCC)  USA300HOU_1050 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,120,7120GC 7.5% 188.2 / 2.1 63R67P (CGC→CCC) USA300HOU_1050hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (26/32);  new base C (3/2);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GAAATCTTGCATCTGATTAAAGTTCAACTAGATAATTTAACATTACCTTCATGCCCATTATATGAAGAAGTACTAGATACACAAATGTTTGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAATATATAAGCATATTAAAATGGCGTGAGATATCGACAATTGAGCTGCGCCGATTTTAT  >  CP000730/1120573‑1120852
                                                                                                                                           |                                                                                                                                            
gaAATCTTGCATCTGATTAAAGTTCAACTAGATAATTTAACATTACCTTCATGCCCATTATATGAAGAAGTACTAGATACACAAATGTTTGGACTTCCAAAAGAAGTTGATTTTGCTGGTAAATTAGGTTTAGGTGACcgc                                                                                                                                             >  2:767288/1‑141 (MQ=255)
gaAATCTGGCATCCGATTAAAGTTCAACTAGATAATTTAACATTACCTTCATGCCCATTATATGAAGAAGTACTAGATACACAAATGTTTGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACcgc                                                                                                                                             <  2:574461/141‑1 (MQ=255)
   aTCTTGCATCTGATTAAAGTTCAACTAGATAATTTAACATTACCTTCATGCCCATTATATGAAGAAGTACTAGATACACAAATGTTTGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGaa                                                                                                                                          <  1:208017/141‑1 (MQ=255)
       tGCATCTGATTAAAGTTCAACTAGATAATTTAACATTACCTTCATGCCCATTATATGAAGAAGTACTAGATACACAAATGTTTGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGCTTTAGTTGACCGCGAAGATg                                                                                                                                      <  1:524903/141‑1 (MQ=255)
       tGCATCTGATTAAAGTTCAACTAGATAATTTAACATTACCTTCATGCCCATGATATGAAGAAGTACTAGATACACAAATGTTTGGGCTTCAAAAAGAAGTTGATTTTGCTGTTAAGTTCGGTGTAGTTGACAGCGAAGATg                                                                                                                                      >  1:237166/1‑141 (MQ=255)
           tCTGATTAAAGTTCAACTAGATAATTTAACATTACCTTCATGCCCATTATATGAAGAAGTACTAGATACACAAATGTTTGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGcaa                                                                                                                                  >  2:753720/1‑141 (MQ=255)
                    aGTTCAACTAGATAATTTAACATTACCTTCATGCCCATTATATGAAGAAGTACTAGATACACAAATGTTTGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTAt                                                                                                                         >  1:691753/1‑141 (MQ=255)
                    aGTTCAACTAGATAATTTAACATTACCTTCATGCCCATTATATGAAGAAGTACTAGATACACAAATGTTTGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTAt                                                                                                                         >  2:844042/1‑141 (MQ=255)
                     atacAAATAGAAAATTAAACAACACCTTCATGCCCATTATATAAAGAAGTACGAGATACACAAAAAAAAAAAAAAAAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATg                                                                                                                        <  2:567133/138‑1 (MQ=255)
                      ttCAACTAGATAATTTAACATTACCTTCATGCCCATTATATGAAGAAGTACTAGATACACAAATGTTTGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGGCGCCGACGCTGGCAAGCAAAttggat                                                                                                                       >  1:279439/1‑137 (MQ=255)
                      ttCAACTAGATAATTTAACATTACCTTCATGCCCATTATATGAAGAAGTACTAGATACACAAATGTTTGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCCCGAAGATGGCAAACAAATTATGt                                                                                                                       >  2:112003/1‑141 (MQ=255)
                      ttCAACTAGATAATTTAACATTACCTTCATGCCCATTATATGAAGAAGTAATAGAAACACAAATGTTTGGCATTAATCTTGACGAAAAAGTCTCGGTTAAACGAGCGTTAGTTGACCCCGAAGATGGCAAACAAATTAGGt                                                                                                                       >  2:392933/1‑141 (MQ=255)
                         aaCTAGATAATTTAACATTACCTTCATGCCCATTATATGAAGAAGTACTAGATACACAAATGTTTGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTAc                                                                                                                    >  2:57448/1‑141 (MQ=255)
                           cgagacaagtaaaCAGTACCTGCATGCCCATGATATGACGAAGTACTAGATACACAAATGTTTGGACTTCAAAAAGAAGGTGATTTTGCTGTTAAATTAGGTTTAGTTGACTGCGAAGATGGCAAACAAATTATGTTACGt                                                                                                                  <  2:145377/130‑1 (MQ=255)
                           cTAGATGATTTAACATTACCGTCATGCCCATTATATGAAGAAGTACTAGATACACAAATGTTTGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGt                                                                                                                  <  1:840862/141‑1 (MQ=255)
                                     tAACATTACCTTCATGCCCATTATATGAAGAAGTACTAGATACACAAATGTTTGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGagaaag                                                                                                        <  1:515896/141‑1 (MQ=255)
                                     tAACATTACCTTCATGCCCATTATATGAAGAAGTACTAGATACACAAATGTTTGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGagaaag                                                                                                        >  1:815611/1‑141 (MQ=255)
                                      aaCATTACCTTCATGCCCATTATATGAAGAAGTACTAGATACACAAATGTTTGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGagaaaga                                                                                                       >  1:435874/1‑141 (MQ=255)
                                            aCCTTCATGCCCATTATATGAAGAAGTACTAGATACACAAATGTTTGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTc                                                                                                 <  1:753720/141‑1 (MQ=255)
                                              cTTCATGCCCATTATATGAAGAAGTACTAGATACACAAATGTTTGGCCTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTGTCaa                                                                                               >  1:308523/1‑141 (MQ=255)
                                                    gCCCATTATATGAAGAAGTACTAGATACACAAATGTTTGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGTCCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTAc                                                                                         <  1:767288/141‑1 (MQ=255)
                                                    gCCCATTATATGAAGAAGTACTAGATACACAAATGTTTGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTAc                                                                                         <  2:705655/141‑1 (MQ=255)
                                                     cccATTATATGAAGAAGTACTAGATACACAATTGTTTGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACa                                                                                        <  2:439839/141‑1 (MQ=255)
                                                       cATTATATGAAGAAGTACTAGATACACAAATGTTTGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATg                                                                                      <  2:299308/141‑1 (MQ=255)
                                                        aTTATATGAAGAAGTACTAGATACACAAATGTTTGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGa                                                                                     >  1:244833/1‑141 (MQ=255)
                                                           ataagagaatgTACTAGATACACATATGTTTGTACTTCAAAAAGAAGTTGATTTGGCTGTTAAATTAGGTTTAGGTGACCGCGAAGATGTCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGc                                                                                  <  2:26821/131‑1 (MQ=255)
                                                            tatGAAGAAGTACTAGGTACACAAATGTTTGGAGTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCt                                                                                 <  2:689038/141‑1 (MQ=255)
                                                               aaacaagTTCTGGATACACAACGGTTTGGACTTCAAAAAGAAGTTGCTTTTGCTGTTAAATTACGTTTAGTTGACCGCGACGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCtttt                                                                              <  2:205029/137‑1 (MQ=255)
                                                                   aagTACTAGATACACAAATGTTTGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTacac                                                                          <  2:556845/141‑1 (MQ=255)
                                                                   aagTACTAGATACACAAAGGTTTGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTacac                                                                          <  1:680422/141‑1 (MQ=255)
                                                                       aCTAGATACACAAATGTTTGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGt                                                                      <  2:257957/141‑1 (MQ=255)
                                                                                   aaTGTTTGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAat                                                          >  2:466281/1‑141 (MQ=255)
                                                                                    aTGTTTGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAata                                                         >  2:788299/1‑141 (MQ=255)
                                                                                        ttGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAAtatata                                                     <  1:139667/141‑1 (MQ=255)
                                                                                             cTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAATATATAAGCat                                                >  2:608301/1‑141 (MQ=255)
                                                                                                 aaaacGAAATTGATTTTGCTGTTAAATTAGATTTAGTTGACCGCGGAGATTGCAATCAGATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACATGTTTAAGGTTTAATATATAAGCttatta                                            >  2:417118/1‑141 (MQ=255)
                                                                                                 aaaaaGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAATATATAAGCATATTa                                            >  1:77653/1‑141 (MQ=255)
                                                                                                 aaaaaGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAATATATAAGCATATTa                                            >  2:363343/1‑141 (MQ=255)
                                                                                                 aaaaaGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAATATATAAGAATATTa                                            >  2:416824/1‑141 (MQ=255)
                                                                                                  aaaaaaAGTTGAGTTGGCTGTTAAATTAGGTTTAGTTGACCGGGAAGATGGCAAACAAATTATGTTACGGCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAATATATAAGCATATTaa                                           <  2:237166/141‑1 (MQ=255)
                                                                                                   aaaGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAATATATAAGCATATTaaa                                          <  1:125726/141‑1 (MQ=255)
                                                                                                    aagaagTTGATTTTGCTGGTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAATATATAAGCATATTaaaa                                         <  2:802300/141‑1 (MQ=255)
                                                                                                     agaagTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAATATATAAGCATATTAAAAt                                        >  1:553668/1‑141 (MQ=255)
                                                                                                       aagTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAATATATAAGCATATTAAAATgg                                      >  2:401654/1‑141 (MQ=255)
                                                                                                         gTTGATTTTGCTGTTAAATGAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAATATATAAGCATATTAAAATGGCg                                    <  2:271191/141‑1 (MQ=255)
                                                                                                           tgatcttgccactAAATTAGTCTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTAAGTCCTGAGAAAGAACTTTCAAAATTACAAGAAGCTTTTACACTTGTTTAAGGTTTAATATATAAGCATATTAAAATGGCGTg                                  <  2:298404/129‑1 (MQ=255)
                                                                                                            gATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAATATATAAGCATATTAAAATGGCGTga                                 >  1:574037/1‑141 (MQ=255)
                                                                                                             aTTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAATATATAAGCATATTAAAATGGCGTgag                                <  2:168340/141‑1 (MQ=255)
                                                                                                                ttgttgTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGGTTAATATATACGCATATTAAAATGGCGTGAGata                             <  2:733057/141‑1 (MQ=255)
                                                                                                                ttcgtttttaaTTACTTTGAGTTGACCGTGAAGATTGCAAACAAATTTTGTTACGTCTTGAGAAATAACTTTCAAAATTACATGAAGCCTTTACACGTGTTTAAGGTTTAATATATAAGCATATTAAAATGGCGTGAGata                             <  2:306705/132‑1 (MQ=255)
                                                                                                                ttGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAATATATAAGCATATTAAAATGGCGTGAGata                             <  2:815611/141‑1 (MQ=255)
                                                                                                                  ggtGTTAAATTAGGTGTAGTTGCCCGAGCAGCTGGCAAACAAATCATGGTACGGCGTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAATATATAAGCATATTAAAATGGCGTGAGATATc                           <  2:595464/139‑1 (MQ=255)
                                                                                                                  gCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAATATATAAGCATATTAAAATGGCGTGAGATATc                           >  2:403039/1‑141 (MQ=255)
                                                                                                                   cTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTCAGGTTTAATATATAAGCATATTAAAATGCTGGGAGATATCg                          >  2:151142/1‑141 (MQ=255)
                                                                                                                   cTGTTAAAATCGGTCTATTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAATATATAAGCATATTAAAATGGCGTGAGATATCg                          >  2:840064/1‑141 (MQ=255)
                                                                                                                        aaaTTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAATATATAAGCATATTAAAATGGCGTGAGATATCGACAAt                     >  2:630813/1‑141 (MQ=255)
                                                                                                                             aGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAATATATAAGCATATTAAAATGGCGTGAGATATCGACAATTGAGc                <  1:57448/141‑1 (MQ=255)
                                                                                                                              ggTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGATTTAATATATAAGCATATTAAAATGGCGTGAGATATCGACAATTggggt               >  2:340474/1‑137 (MQ=255)
                                                                                                                                tttAGTTTGCCGCGAAGATGGCAAACAAATTATGTTACGTCTTGATAAAGAACTTTCAATATTACATGAAGCTTTTACACTTGTTTAAGTTTTAATATATAAGCATATTAAAATGGCGTGAGATATCGACAATTGAGCTgc             <  2:847639/141‑1 (MQ=255)
                                                                                                                                       gACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAATATATAAGCATATTAAAATGGCGTGAGATATCGACAATTGAGCTGCGCCGAtt      >  1:503252/1‑141 (MQ=255)
                                                                                                                                       gACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAATATATAAGCATATTAAAATGGCGTGAGATATCGACAATTGAGCTGCGCCGAtt      <  1:605456/141‑1 (MQ=255)
                                                                                                                                       gACCCCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAATATATAAGCATATTAAAATGGCGTGAGATATCGACAATTGAGCTGCGCCGAtt      <  1:112003/141‑1 (MQ=255)
                                                                                                                                       gACCCCGAAGATGGCAAACAAAGTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAATATATAAGCATATTAAAATGGCGTGAGATATCGAAAATTGAGCTGCGCCGAtt      <  1:392933/141‑1 (MQ=255)
                                                                                                                                         ccGCGAAGATGGCAAACAATGTAGGTTACGTCTTGAGAAAGAACTTTCAAAATCACATGAAGCTTTTACACTTGTTGAAGGTTTAATATACACGCAGATTAAAATGGCGTGAGATATCGACAATTGAGCTGCGCCGAtttt    <  1:166505/141‑1 (MQ=255)
                                                                                                                                         ccGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAATATATAAGCATATTAAAATGGCGTGAGCTATCGACAATTGCGCTCGGCCGAtttg    >  2:378139/1‑140 (MQ=255)
                                                                                                                                           gcgAAGATGGCAGACAAATTATGTTACGTCTTGAGAATTGACTTTCGAAAGTACAGGAAGCTTTTACACTTGTCTAAGGTTTAATATATAAGCATATTAAAATGGCGTGAGATAGCGACAATGGATCTGCGCCGAatttct  >  2:568444/1‑139 (MQ=255)
                                                                                                                                           |                                                                                                                                            
GAAATCTTGCATCTGATTAAAGTTCAACTAGATAATTTAACATTACCTTCATGCCCATTATATGAAGAAGTACTAGATACACAAATGTTTGGACTTCAAAAAGAAGTTGATTTTGCTGTTAAATTAGGTTTAGTTGACCGCGAAGATGGCAAACAAATTATGTTACGTCTTGAGAAAGAACTTTCAAAATTACATGAAGCTTTTACACTTGTTTAAGGTTTAATATATAAGCATATTAAAATGGCGTGAGATATCGACAATTGAGCTGCGCCGATTTTAT  >  CP000730/1120573‑1120852

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: