Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,018,834 G→C 11.2% P29A (CCA→GCA)  USA300HOU_1869 ← transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,018,8340GC11.2% 166.8 / 10.8 64P29A (CCA→GCA) USA300HOU_1869transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (26/28);  new base C (3/4);  total (30/34)
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.58e-01

TTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAACTCACTTAGTATATACCTTGAGTTTACCAGTAC  >  CP000730/2018696‑2018973
                                                                                                                                           |                                                                                                                                           
ttCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGtt                                                                                                                                           <  1:292279/141‑1 (MQ=255)
 tCTTAGTTCATCCACATTATCGCTTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTg                                                                                                                                          <  1:240691/141‑1 (MQ=255)
 tCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGttt                                                                                                                                          >  2:536535/1‑140 (MQ=255)
 tCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTg                                                                                                                                          <  2:12332/141‑1 (MQ=255)
  cTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGCTAATGATTAGAATGGTTAATAATATTTTTAATAGTCTTCGTTTAATTTTTTTCTTTTTTTGTGGTTTGCAAATTGGGTGa                                                                                                                                         >  2:371671/1‑141 (MQ=255)
  cTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAATTGGTTgt                                                                                                                                         >  1:503001/1‑140 (MQ=255)
  cTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACGGGTTGa                                                                                                                                         >  1:91052/1‑141 (MQ=255)
     aGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGataa                                                                                                                                      <  2:669544/141‑1 (MQ=255)
     aGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGAGATAATGATTAGAATGGTTAATAATATTTTTAATAGTCTTCGTTAACTTTTTTTCTTTTTGGGCGTTTTGGCACGTGGTTGataa                                                                                                                                      >  1:230204/1‑141 (MQ=255)
         cATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATCATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGTTGGTTTGCCAACTGGTTGTTAATACg                                                                                                                                  >  1:203313/1‑141 (MQ=255)
         cATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACg                                                                                                                                  >  2:132082/1‑141 (MQ=255)
            ccACATTATCGCGTGTAGCTAAAAAATACATGATACCAATAAACAATGCGATAACGATTAGAGTGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGTCGGTTTGTCACTTGGTGGCTAATACGtat                                                                                                                               >  1:74138/1‑141 (MQ=255)
            ccACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGtat                                                                                                                               >  2:23707/1‑141 (MQ=255)
                   aTCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATAAGTATTATAGTg                                                                                                                        >  1:332824/1‑141 (MQ=255)
                       cgTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGt                                                                                                                    <  1:357872/141‑1 (MQ=255)
                       cgTGTAGATAAAAAATACATGATACAAATTAACAATGCGATGCTGCTTTGAATGCTTAATAATATTTTAATCGTTTGTCGCTTAACCACGCGCTGTGTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGt                                                                                                                    <  1:389519/141‑1 (MQ=255)
                               tAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGttt                                                                                                            <  1:196314/141‑1 (MQ=255)
                                   aaaTACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTTTTTTTCTTCTATTTGCCAACTGCTTTATAATACGTATTATAGTGAGGTTGGTGTTTCatt                                                                                                        >  2:405223/1‑140 (MQ=255)
                                   aaaTACATGATACCAATAAACAATGCGATAATGATTACAATGGTTAATAATATTTTTAATAGAATTCGTTTACTTTTTTTGTTATTCGGCGGTTTTCCAACAGCTTGCTAATACGTCTTATCGTTCGGTTCGTGTGTGata                                                                                                        >  2:223650/1‑141 (MQ=255)
                                      taaaTGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGc                                                                                                     <  1:709109/138‑1 (MQ=255)
                                        cATGATACCAATATACAATGCGATAATGATTCGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGGTGATAATACGTATTATAGTGAGGTTCGTGTTTCATTTGCTc                                                                                                   <  1:310326/141‑1 (MQ=255)
                                        cATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTc                                                                                                   >  1:43962/1‑141 (MQ=255)
                                        cATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTc                                                                                                   <  1:151107/141‑1 (MQ=255)
                                        cATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGGTTGCCAACTGGTTGATAATACGCATTATAGTGAGGTTCGTGTTTCATTTGCTc                                                                                                   <  1:310653/141‑1 (MQ=255)
                                        cATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTGCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTc                                                                                                   <  1:384248/141‑1 (MQ=255)
                                         aTGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATGTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCa                                                                                                  >  1:426881/1‑141 (MQ=255)
                                           taTACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCaaa                                                                                                <  1:472037/140‑1 (MQ=255)
                                            aTACCAAGAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCaaaa                                                                                               <  1:452340/141‑1 (MQ=255)
                                                  aTAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCa                                                                                         >  1:354711/1‑141 (MQ=255)
                                                  aTAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCATTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCa                                                                                         <  1:550128/141‑1 (MQ=255)
                                                  aTAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGGATTCGTTTACTTTTTTTCTTTTTAGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGAGCTTCAGAGGGTCAAAATGTTCa                                                                                         >  2:524106/1‑141 (MQ=255)
                                                    aaaCAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTTTTTTTCGTCGGTTTGCCAACTGTTTGATAATATGTATTATTGTGAGGTTCGTGTTTCATTTGCTCAATATGTTCAtt                                                                                       >  2:333627/1‑141 (MQ=255)
                                                    aaaCAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGCTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCAtt                                                                                       >  1:404803/1‑141 (MQ=255)
                                                        attgagATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCAtttgag                                                                                   <  1:78227/136‑1 (MQ=255)
                                                          tGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGCTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCAtttgagtt                                                                                 <  1:331454/141‑1 (MQ=255)
                                                          cgCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGCTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCAtttgagtt                                                                                 <  1:41752/140‑1 (MQ=255)
                                                             gATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTGGTGTTTCATATGCTGAGACTGTGCAtttgagtttgg                                                                              >  1:69607/1‑140 (MQ=255)
                                                             gATAATGATTAGAATGGTTAATAATATTATTAATAGTAATCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTACATATGCTCAAAATGTTCAtttgagtttga                                                                              >  2:739820/1‑141 (MQ=255)
                                                              aTAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGGGGTTCGGGTTTCATATGCTCAAAATGTTCAtttgagtttgag                                                                             >  1:29401/1‑141 (MQ=255)
                                                              aTAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCAtttgagtttgag                                                                             >  1:733734/1‑141 (MQ=255)
                                                                aaTGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTCATGCGTCGGTTTTCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTa                                                                           <  1:587446/141‑1 (MQ=255)
                                                                       aGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCg                                                                    <  1:401177/141‑1 (MQ=255)
                                                                         aaTGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGct                                                                  <  1:112142/141‑1 (MQ=255)
                                                                                  taataTTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCTTCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCCAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCAt                                                         >  1:266854/1‑141 (MQ=255)
                                                                                    atatTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGc                                                       <  1:302441/141‑1 (MQ=255)
                                                                                     tatTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCg                                                      >  1:371559/1‑141 (MQ=255)
                                                                                     tatTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCg                                                      >  1:163451/1‑141 (MQ=255)
                                                                                      atTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGt                                                     >  2:684295/1‑141 (MQ=255)
                                                                                      atTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGGTTGAGTACCTATCGCTTCTTTTCATGCGt                                                     <  1:392214/141‑1 (MQ=255)
                                                                                       tttttAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTCAGTTTGAGTACCTATCGCTTCTTTTCATGCAtt                                                    >  2:239009/1‑141 (MQ=255)
                                                                                                       tttttCTTTTTGTCTTTTTCGGCGGTTTGTCAACTTGTTGATAATACGTATTATAGTGCGGTTCGTGTTTCATTTGCTCATCATGTTCATTTGAGTTTGAGTATCTTCCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTaa                                     <  2:303432/140‑1 (MQ=255)
                                                                                                           aCTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAACt                                 <  1:172816/141‑1 (MQ=255)
                                                                                                           aCTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAAAACGTGTTATAGTAAGATTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTAACTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAAca                                 >  1:277609/1‑140 (MQ=255)
                                                                                                           aCTTGTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAACt                                 >  2:560451/1‑141 (MQ=255)
                                                                                                                ttttctttttcGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAACTCACtt                            <  1:634353/141‑1 (MQ=255)
                                                                                                                    ctttttcGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGATCCTTCTTTTAAAACTCACTTAGta                        >  1:660088/1‑141 (MQ=255)
                                                                                                                           gggggTTTGGCAACTGTTTGATAATACGTATTATAGTGGGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTTCCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAACTCACTTAGTATATACCt                 <  2:481920/138‑1 (MQ=255)
                                                                                                                            gctgTTTGCCAATTGGTTAATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAACTCACTTAGTATATACCtt                <  2:203313/138‑1 (MQ=255)
                                                                                                                            gcggTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAACTCACTTAGTATATACCtt                <  1:132082/141‑1 (MQ=255)
                                                                                                                             cggTTTGCCAATTGTTTGATAATACGTATTATAGGGAGGTTTGTTTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAACCTCACTTAGTATATACCTTg               <  2:282603/141‑1 (MQ=255)
                                                                                                                             cggTTTGCCAACTGGTTGATAAGACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAACTCACTTAGTATATACCTTg               <  1:303866/141‑1 (MQ=255)
                                                                                                                                tttGCCAACTGCTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAACTCACTTAGTATATACCTTGAGt            <  2:404803/141‑1 (MQ=255)
                                                                                                                                tctGCCAACGGCTTGATAAAACGGATTAAAGTGAGGTGCGTGTTTCATATGCTCAAAATTTTCATTTGAGTTTGAGCACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAACTCACTTAGTATATACCTTGAGt            <  1:333627/139‑1 (MQ=255)
                                                                                                                                   gCCAACTGGTTTATAATACGTATTATAGTGAGGTTCTTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTAACGCTTCTTTTCATGCGTTTGCGTCTTCTTTTAAAACTCACTTAGTATATACCTTGAGTTTa         <  2:332824/141‑1 (MQ=255)
                                                                                                                                          ggTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAACATGTTCGTTTCAATCTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAACTCACTTAGTATATACCTTGAGTTTACCCGTAc  <  1:551433/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                           
TTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAATGCGATAATGATTAGAATGGTTAATAATATTTTTAATAGTATTCGTTTACTTTTTTTCTTTTTCGGCGGTTTGCCAACTGGTTGATAATACGTATTATAGTGAGGTTCGTGTTTCATATGCTCAAAATGTTCATTTGAGTTTGAGTACCTATCGCTTCTTTTCATGCGTTTGCTCCTTCTTTTAAAACTCACTTAGTATATACCTTGAGTTTACCAGTAC  >  CP000730/2018696‑2018973

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: