Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,857,867 C→T 17.3% L169L (TTG→TTA)  USA300HOU_1728 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,857,8670CT17.3% 141.3 / 12.8 65L169L (TTG→TTA) USA300HOU_1728hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/25);  new base T (3/9);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCTTGTGCTTCTCTTTTAGCTTCTTGAATCTCTT  >  CP000730/1857731‑1858003
                                                                                                                                        |                                                                                                                                        
cttTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTTAAat                                                                                                                                      >  1:310014/1‑141 (MQ=255)
    tttctctgctTGTTGTTTTTTCACTTCCACTTTTTTTGGTTCTTCATTTTTTACTGCTGTGTTTTTGTCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAattatt                                                                                                                                  <  2:495253/136‑1 (MQ=17)
    tgtgggtgattgtggTTCTTTCAATTCTACTGTTGCTGGTCTTTCCTTTTATGCTGCTGCACTTTTCTCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCCTGTGCCTCTTGTGATGTATCACTCAAattatt                                                                                                                                  <  2:495271/132‑1 (MQ=12)
    tATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAattatt                                                                                                                                  <  1:330171/141‑1 (MQ=255)
           gctTGTTGTTCTTTCAATGCTACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACtt                                                                                                                           <  1:900805/141‑1 (MQ=255)
            ctTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTTTGGTGTTTCACTCGAGTTATTTGctgttt                                                                                                                          >  2:617816/1‑136 (MQ=17)
                ttgtGCTTTCACTTCTACTGCTTCTGCTCCTTCATTTGTTGCTGCTGTGCTTCTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGcttc                                                                                                                      <  2:536606/141‑1 (MQ=17)
                   tgCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGCGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGcttcttc                                                                                                                   <  2:1018406/141‑1 (MQ=255)
                    gCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGTATCTCTTGTGCCTCTTGTGATGTATCAATCAAATTATTTGCACTTGcttcttct                                                                                                                  >  2:947227/1‑141 (MQ=255)
                    gCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGcttcttct                                                                                                                  >  1:791222/1‑141 (MQ=255)
                      tttCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCttt                                                                                                                >  1:169395/1‑141 (MQ=255)
                         cAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTAt                                                                                                             >  1:974338/1‑141 (MQ=255)
                              ccACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCa                                                                                                                                         >  1:261454/1‑108 (MQ=255)
                              ccACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCa                                                                                                                                         <  2:261454/108‑1 (MQ=255)
                                   gCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCTCTGCttgt                                                                                                   >  2:201558/1‑141 (MQ=255)
                                        tGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGc                                                                                              <  1:977325/141‑1 (MQ=255)
                                        tGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGc                                                                                              >  2:1090050/1‑141 (MQ=255)
                                                          aGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCttt                                                                            >  2:400855/1‑141 (MQ=255)
                                                           gCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCtttt                                                                           >  1:999794/1‑141 (MQ=255)
                                                             tgtgtttttGTCTTTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGa                                                                         <  2:127358/136‑1 (MQ=255)
                                                              gttgttgtgtCTTTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGAt                                                                        <  2:231384/132‑1 (MQ=11)
                                                              gtgtttttGTCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGAt                                                                        <  2:483788/137‑1 (MQ=255)
                                                                gctattgtCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATATATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTc                                                                      <  1:995367/134‑1 (MQ=255)
                                                                    ttgtctttctCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTa                                                                  <  2:834268/132‑1 (MQ=2)
                                                                        cTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTg                                                              <  1:374279/141‑1 (MQ=255)
                                                                             tCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGTTTTTTTTTTTTTTACa                                                         >  2:768056/1‑141 (MQ=14)
                                                                             tCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACa                                                         >  1:957102/1‑141 (MQ=255)
                                                                                gCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATGTATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGcc                                                      <  1:1128883/141‑1 (MQ=21)
                                                                                 cTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGACGTGTCACTTAAATTATTTGCACTTGCTTCTTCTTTTATTGCTGCTTGTTGTACTTTTAATGCCGATGCTTTTGACTCTTCATTTGATACAGCca                                                     >  1:732027/1‑141 (MQ=11)
                                                                                  ttGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCCTTCAATGCCTCTGGTTTTGATT‑TTTAtttttataaagaatc                                                   >  1:613624/1‑129 (MQ=14)
                                                                                     tgCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATATATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACt                                                 <  2:430897/141‑1 (MQ=255)
                                                                                       cttcttTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACttt                                               <  1:1193272/141‑1 (MQ=255)
                                                                                         tcttTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTa                                             <  2:58579/141‑1 (MQ=255)
                                                                                         tcttTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTa                                             <  1:871969/141‑1 (MQ=255)
                                                                                          cttTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTtgttgt                                                                                                <  2:317105/89‑1 (MQ=255)
                                                                                          cttTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTtgttgt                                                                                                >  1:317105/1‑89 (MQ=255)
                                                                                          cttTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGcc                                                      <  2:785706/131‑1 (MQ=255)
                                                                                          cttTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGcc                                                      >  1:785706/1‑131 (MQ=255)
                                                                                           ttttttAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATc                                           <  1:365248/141‑1 (MQ=255)
                                                                                                 aGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCCCCGCTATTGATTCTTTATTTGCTACAGGCACACTTTTAGCCGTTTc                                     >  2:82398/1‑141 (MQ=17)
                                                                                                 aGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTc                                     <  1:101014/141‑1 (MQ=255)
                                                                                                 aGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTc                                     <  1:223356/141‑1 (MQ=255)
                                                                                                   cttcttGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTTTTACAGCCACACTTTTATCCGCTTTTg                                   >  2:1045889/1‑141 (MQ=255)
                                                                                                   cttcttGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTg                                   <  1:470359/141‑1 (MQ=255)
                                                                                                   cttcttGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTg                                   >  2:52444/1‑141 (MQ=255)
                                                                                                     tcttGAATCTCTTGTGCCTCTTGTGATGTATCACTCCAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCt                                 >  1:599906/1‑141 (MQ=255)
                                                                                                     tcttGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGc                                                                                              >  1:407092/1‑80 (MQ=39)
                                                                                                        tGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCTtgt                              >  1:820827/1‑141 (MQ=255)
                                                                                                           aTCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCTTGTGCt                           >  2:57560/1‑141 (MQ=255)
                                                                                                             ctctTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTTTTTATTTTATACAGCCACTCTTTTATCCTCTTCTGCTTGTGCTtc                         >  2:1010304/1‑141 (MQ=255)
                                                                                                              tctTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTTATACAACCACACTTTTATCCGCATCTGCtatagcttcc                        >  1:471350/1‑132 (MQ=17)
                                                                                                                ttGTGCCTCTTGTGATGTATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGTTACAGCCACACTTTTATCCGCTTCTGCTTGTGCTTctct                      <  2:351266/141‑1 (MQ=255)
                                                                                                                ttGTGCCTCTTGTGATGTATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCTTGTGCTTctct                      <  2:195877/141‑1 (MQ=255)
                                                                                                                ttGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCTTGTGCTTctct                      >  2:1163131/1‑141 (MQ=255)
                                                                                                                  gtgCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCTTGTGCTTCTCttt                    >  2:94779/1‑141 (MQ=255)
                                                                                                                  gtgCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCTTGTGCTTCTCttt                    <  1:1142190/141‑1 (MQ=255)
                                                                                                                   tgCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGcttc                                                                                                                      >  1:40425/1‑42 (MQ=255)
                                                                                                                    gCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCTTGTGCTTCTCTTTTa                  <  1:617816/141‑1 (MQ=255)
                                                                                                                    gCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCTTGTGCTTCTCTTTTa                  <  2:765716/141‑1 (MQ=255)
                                                                                                                      ctctTGTGATGTATCACTCATATTATTTGCACTTGCTTCTTCTTTTTTCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATCCAGCCACACTTTTATCCGCTTCTGCTTGTGCTTCTCTTTTAGc                <  2:599906/141‑1 (MQ=255)
                                                                                                                         ttGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCTTGTGCTTCTCTTTTAGcttc             >  1:393988/1‑141 (MQ=255)
                                                                                                                           ttgATGTATCACTTAAATTATTTGCACTTGCTTCTTCTTTTTTCGCTTCTTGTTGTGCTTTCAATGCCACTGCTTTCCATTCTTCATTCCCTACAGCCACACTTTTATCCGCTTCTGCTTGTGCTTCTCTTTTAGcttctt           <  2:1083312/140‑1 (MQ=14)
                                                                                                                              aTGTATCACTCAAATTATTTGCACTTGCTTCTTGTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCTTGTGCTTCTCTTTTAGCTTCTTGaa        >  1:1072085/1‑141 (MQ=255)
                                                                                                                                atATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCTTGTGCTTCTCTTTTAGCTTCTTGAAtc      <  1:782291/140‑1 (MQ=255)
                                                                                                                                    cACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCTTGTGCTTCTCTTTTAGCTTCTTGAATCTCtt  <  1:659821/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                        
CTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCTTGTGCTTCTCTTTTAGCTTCTTGAATCTCTT  >  CP000730/1857731‑1858003

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