Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,912,184 T→G 32.2% intergenic (+1717/+366) USA300HOU_1773 → / ← USA300HOU_1775 camphor resistance protein/possible aldo/keto reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,912,1840TG32.2% 64.0 / 57.6 66intergenic (+1717/+366)USA300HOU_1773/USA300HOU_1775camphor resistance protein/possible aldo/keto reductase
Reads supporting (aligned to +/- strand):  ref base T (33/11);  new base G (10/11);  total (43/23)
Fisher's exact test for biased strand distribution p-value = 4.85e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

AAAGATATGCATACCGTAATGATAGATATGGTTTTAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGTTTTTGGATTTATGAAGGCTATTTTGGGTTTAACTCGAACGT  >  CP000730/1912045‑1912323
                                                                                                                                           |                                                                                                                                            
aaaGATATGCATACCGTAATGATAGATATGGTTTTAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAtt                                                                                                                                             <  1:1040713/141‑1 (MQ=255)
  aGATATGCATACCGTAATGATAGATATGGTTTTAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACCATGCATGAAACCAAATTCAAAATCCAATAGGAAAATCATGAAGAATtat                                                                                                                                           >  1:299601/1‑141 (MQ=255)
  aGATATGCATACCGTAATGATAGATATGGTTTTAAACGAGACTTCAAACTATATGAAGGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAtaca                                                                                                                                           >  2:35578/1‑138 (MQ=255)
   gATATGCATACCGTAATGATAGATATGGTTTTAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATtata                                                                                                                                          >  2:618516/1‑141 (MQ=255)
   gATATGCATACCGTAATGATAGATATGGTTTTAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAGtata                                                                                                                                          >  2:1045666/1‑141 (MQ=255)
     tatgctatacGTAATTATATATATGGTTTTAAACGAGACTTCAAACCACAGGAATGCGTGGAATGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAAt                                                                                                                                        <  2:895833/132‑1 (MQ=14)
      atGCATACCGTAATGATAGATATGGTTTTAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAGTATaagt                                                                                                                                       >  1:87847/1‑139 (MQ=255)
           tACCGTAATGATAGATATGGTTTTAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAGTATAATTGGGaa                                                                                                                                  <  1:497489/141‑1 (MQ=255)
           aaCCGTAATGATAGATATGGTTTTAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAGTATAATTGGGaa                                                                                                                                  <  1:662137/140‑1 (MQ=255)
             ccGTAATGATAGATATGGTTTTAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGAACAAATACATGAATCCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATa                                                                                                                                >  1:274205/1‑141 (MQ=14)
              cGTAATGATAGATATGGTTTTAAACGGGACTTCAAACTATATGAATGCGATGACCGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATAc                                                                                                                               <  2:641291/141‑1 (MQ=255)
                tAATGATAGATATGGTTTTAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAACTTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACtt                                                                                                                             >  1:83185/1‑141 (MQ=255)
                  aTGATAGATATGGTTTTAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTa                                                                                                                           >  1:958280/1‑141 (MQ=255)
                             ggTTTTAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAGTATAATTGGGAATACTTTAAAATCCAAAtt                                                                                                                >  1:1032117/1‑141 (MQ=255)
                                  tAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCa                                                                                                           >  1:853168/1‑141 (MQ=255)
                                  tAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAGTATAATTGGGAATACTTTAAAGTCCAAATTAATCa                                                                                                           <  1:1192838/141‑1 (MQ=255)
                                       gagaCTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGc                                                                                                     >  2:973206/1‑141 (MQ=255)
                                        cgaCTTCAAACTATATGAATGCGCTGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGACTCTTTTCATGTCCAAATTAATCAAAAGCt                                                                                                    >  2:243642/2‑141 (MQ=14)
                                         gaCTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCtt                                                                                                   >  1:743038/1‑141 (MQ=255)
                                             tCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAACTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTg                                                                                               <  2:429913/141‑1 (MQ=255)
                                             tCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTATAAGGCACCAATGAATGAAACAAAATACAAAATCCAATAAGAAAAACAGCAAAAATTAGAACCCGGAAAACATTAAAGTCCAAATAAATCAAAAGCTTTCta                                                                                               >  2:280988/1‑140 (MQ=255)
                                              cAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAGTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGa                                                                                              <  1:182517/141‑1 (MQ=255)
                                                     tatGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGaa                                                                                       >  1:503322/1‑141 (MQ=255)
                                                     tatGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGaa                                                                                       >  1:1106420/1‑141 (MQ=255)
                                                     tatGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTATAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGaa                                                                                       <  1:209885/141‑1 (MQ=255)
                                                       tGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATAAAAAGCTTTCTGAACAAGaaaa                                                                                     >  2:774313/1‑140 (MQ=255)
                                                            gCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGaaaa                                                                                <  1:293496/141‑1 (MQ=255)
                                                             cGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGaaaaa                                                                               >  2:998919/1‑141 (MQ=255)
                                                               aTGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGaaaaaaa                                                                             <  1:1185927/141‑1 (MQ=255)
                                                               aTGACTCGTCAGCATGTTCTTTGAGGCAACAATGCATGAAACAAAATTCAAAATCCAATAAGAAAATCATGAAAAATTATAATTGGGAATACTTTGAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAACCGaaaaaaa                                                                             >  2:305962/1‑141 (MQ=14)
                                                                 gACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCACTAAGAAAATCATGAAGAAGTATAATTGGGAATACTTTAAAGTACAAATAAAACAAAAGCTTTCTGAACCAGAAAAGaaaaagaac                                                                           >  2:841993/1‑137 (MQ=14)
                                                                 gACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATc                                                                           >  1:634779/1‑141 (MQ=255)
                                                                 gACTGTTCACCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCGAAATCCAATAAGAAAATCATGAAGAAGTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATc                                                                           >  2:770177/1‑141 (MQ=255)
                                                                  aCTGTTCAGCATGTTCTTTGAGGCAACACTGCTCGAAACCAAATTCAAAATCCAATAAGAAAATGATGAAGAAGTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCt                                                                          >  2:10585/1‑141 (MQ=14)
                                                                  aCTGTTCAGCATGTGCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGCGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCt                                                                          <  1:334648/141‑1 (MQ=255)
                                                                      ttCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCTATAg                                                                      >  1:714988/1‑141 (MQ=255)
                                                                          gCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAGTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCTATAGTCaa                                                                  <  1:430092/141‑1 (MQ=255)
                                                                             tGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAATCTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGa                                                               >  1:923410/1‑141 (MQ=255)
                                                                             tGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATGGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCTATAGTAAAAGa                                                               >  1:754390/1‑141 (MQ=255)
                                                                              gTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGaa                                                              >  2:1034240/1‑141 (MQ=39)
                                                                               ttCTTTGAGGCAACAATGCGTGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGGATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGaaa                                                             >  1:718472/1‑141 (MQ=255)
                                                                               ttCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGaaa                                                             >  1:110947/1‑141 (MQ=38)
                                                                                   ttGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAAtt                                                         >  1:971663/1‑141 (MQ=39)
                                                                                   ttGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAAtt                                                         >  2:564372/1‑141 (MQ=39)
                                                                                   ttGAGGCAACAATGCATGAAAACAAATTCAAAATACAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACTAAAAAAATCTATAGTCAAAGAAAAAtt                                                         <  2:1108149/141‑1 (MQ=255)
                                                                                       ggCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCCTGAAGACTTATAATTGGGAATACTTTCAAGTCCAAGCACTCCAGAGTCGTTCTATAGCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGAta                                                     >  2:843177/1‑140 (MQ=11)
                                                                                       ggCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATg                                                     >  1:422967/1‑141 (MQ=37)
                                                                                       ggCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAGTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATg                                                     >  2:957534/1‑141 (MQ=37)
                                                                                        gCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAGTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGt                                                    <  1:1095157/141‑1 (MQ=25)
                                                                                           acaaTGCATGAAACCAAATTCAAAATCCGATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAATTATCTATAGTCAAAGAAAAATTGATGTaga                                                 <  2:383047/141‑1 (MQ=12)
                                                                                             aaTGCATGAAACCAAATTCAAAATCAAAGAAGAAAATCATGAAGAAGTATAAGTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACAAAAAAAAAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGc                                               <  2:87847/141‑1 (MQ=11)
                                                                                               tGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAAAGGCt                                             >  2:1141559/1‑141 (MQ=17)
                                                                                                gCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAGTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTg                                            <  1:1045666/141‑1 (MQ=17)
                                                                                                 cATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGt                                           <  1:956947/141‑1 (MQ=32)
                                                                                                    gAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCACAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGtttt                                        <  1:398903/141‑1 (MQ=17)
                                                                                                    gAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAAGTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGtttt                                        <  1:1101268/141‑1 (MQ=17)
                                                                                                      aaCCCACTTCAAAACCCAATAAGGAAATCATGAAGAAGGATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGTTTTTg                                      <  1:1169909/141‑1 (MQ=11)
                                                                                                       aCCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGTTTTTgg                                     >  1:644331/1‑141 (MQ=32)
                                                                                                       aCCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAAAGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGTGTTTgg                                     >  2:150913/1‑141 (MQ=17)
                                                                                                         cAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGTTTTTGGAt                                   >  1:1023118/1‑141 (MQ=32)
                                                                                                            aTTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAACCTATAGTCAAAGAAAAACGGATGTAGAGCCTGTTTTGGGATTTa                                >  2:393580/1‑141 (MQ=255)
                                                                                                              tCAAAATCCAATAAGAAAATCATGAAGAAGTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGTTTTTGGATTTATg                              >  2:602978/1‑141 (MQ=17)
                                                                                                                                aTCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGATATCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAGATTGGTGTAGAGCCTGTTTTTGGATTTATGAAGGCTATTTTTGGTTTa            >  1:996914/1‑141 (MQ=11)
                                                                                                                                  cATGAAGAAGTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGTTTTTGGATTTATGAAGGCTATTTTGGGGTTAAc          >  2:1061251/1‑141 (MQ=17)
                                                                                                                                   aTGAAGAAGTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGTTTTTGGATTTATGAAAGCTATTTTTGGTTTAACt         >  1:400355/1‑141 (MQ=12)
                                                                                                                                     gaagaaTTATAATTGGAAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGTTTTTGGATTTATGAAGGCTATTTTGGGTTTAACTCg       >  1:592895/1‑141 (MQ=17)
                                                                                                                                       agaaGTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGTTTTTGGATTTATGAAGGCTATTTTGGGTTTAACTCGaa     <  1:108894/141‑1 (MQ=17)
                                                                                                                                          aTTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGTTTTTGGATTTATGAAGGCTATTTTGGGTTTAACTCGAACGt  <  2:379662/141‑1 (MQ=35)
                                                                                                                                           |                                                                                                                                            
AAAGATATGCATACCGTAATGATAGATATGGTTTTAAACGAGACTTCAAACTATATGAATGCGATGACTGTTCAGCATGTTCTTTGAGGCAACAATGCATGAAACCAAATTCAAAATCCAATAAGAAAATCATGAAGAATTATAATTGGGAATACTTTAAAGTCCAAATTAATCAAAAGCTTTCTGAACCAGAAACGAAAAAAATCTATAGTCAAAGAAAAATTGATGTAGAGCCTGTTTTTGGATTTATGAAGGCTATTTTGGGTTTAACTCGAACGT  >  CP000730/1912045‑1912323

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