Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,807,163 C→T 6.6% G487R (GGG→AGG)  USA300HOU_2651 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,807,1630CT 6.6% 200.8 / 2.7 75G487R (GGG→AGG) USA300HOU_2651hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (34/34);  new base T (2/3);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.43e-01

GTTTGTGAAATATAATTCCAAGTAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCTATGGCAAAAACAGTTTGTGACAATGATTTCTTTTGA  >  CP000730/2807023‑2807299
                                                                                                                                            |                                                                                                                                         
gTTTGTGAAATATAATTCCAAGTAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCTTCTTTATTTGGccc                                                                                                                                           >  1:355164/1‑141 (MQ=255)
   tgtcAAATATAATTCCAAGGAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAgg                                                                                                                                        <  1:406778/141‑1 (MQ=255)
   tactaaaTAAGATTCCAAGTAACTTGAAATGTTTTGTATTTTAGCAAACTGAACTCTCCCAAACGTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAgg                                                                                                                                        <  1:567997/137‑1 (MQ=255)
    gtgAAATATAATTCCAAGTAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTCTGTCGCCCAGGa                                                                                                                                       >  1:813091/1‑141 (MQ=255)
    gtgAAATATAATTCCAAGTAACTTTAAATGTTTTGTTTTTCAGCATGTTAAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACAAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGa                                                                                                                                       >  2:620152/1‑141 (MQ=255)
     tgAAATATAATTCCAAGTAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCGACGaa                                                                                                                                      >  1:61517/1‑141 (MQ=255)
     tgAAATATAATTCCAAGTAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCAACATTACGTCGCCCaagca                                                                                                                                      >  2:118645/1‑137 (MQ=255)
      gAAATATAATTCCAAGTAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCACCATTATGTCGCCCAGGAAt                                                                                                                                     <  1:545677/141‑1 (MQ=255)
      gAAATATAATCCCAAGTAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTAAATACTTTTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAAt                                                                                                                                     >  1:607963/1‑141 (MQ=255)
          tataATTCCAAGTAACTTTAAATGTTTTGTTTTTCTGCATGTTGAACTCTCCCAAACTTGTCTTCCAAATTAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTg                                                                                                                                 >  2:234033/1‑141 (MQ=255)
             aaTTCCAAGTAACTTTAAATGTTTTGTTTTTCAGCATATTGAACGCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAAGGGGAGGTGAATCATCATTATGTCGCCACGAAATACTGCGa                                                                                                                              >  2:457934/1‑141 (MQ=255)
              tttCTAAGTAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGAt                                                                                                                             <  2:247592/140‑1 (MQ=255)
                tCCAAGTAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCGTTATGTCGCCCAGGACTACTGCGATtc                                                                                                                           >  1:392619/1‑141 (MQ=255)
                  cAAGTACCTTTAATTGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAATATAATGTTGTAAAAATTAACAAACCAACTTGCTATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATtcat                                                                                                                         <  2:27677/141‑1 (MQ=255)
                   aaGTAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATtcatc                                                                                                                        <  2:443037/141‑1 (MQ=255)
                     gtgaccttAAATGTTTTCTATTTCAGCATGTTGAACTCTCCCACACTTGTCTTCCAAAAAAAAAAAATAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCAc                                                                                                                      <  2:744907/135‑1 (MQ=255)
                     gTAACTTTAATTGTTTTGTTTTTCAGCATGATGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCAc                                                                                                                      <  2:668072/141‑1 (MQ=255)
                     gTAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCACAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCAc                                                                                                                      >  2:139636/1‑141 (MQ=255)
                     gTAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCAc                                                                                                                      <  1:693890/141‑1 (MQ=255)
                          tttAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTcttg                                                                                                                 >  1:580358/1‑141 (MQ=255)
                              aaTGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAACACTGCGATTCATCACTcttgcttg                                                                                                             <  1:533739/141‑1 (MQ=255)
                                tgttttgtttTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGgt                                                                                                           >  1:545428/1‑141 (MQ=255)
                                   tttgttttTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGtgtg                                                                                                        <  2:179232/141‑1 (MQ=255)
                                      gttttTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTg                                                                                                     <  1:624318/141‑1 (MQ=255)
                                      gttttTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTg                                                                                                     <  2:418316/141‑1 (MQ=255)
                                       tttttCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCTAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGt                                                                                                    >  2:177457/1‑141 (MQ=255)
                                       tttttCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCTAGGAATACTGCGATCCATCACTCTTGCTTGGTGTGCTGt                                                                                                    >  2:320870/1‑141 (MQ=255)
                                        ttttCAGCATTTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAGAATTAACAAACCAACTTGCAATGGTAGGTGAATAATCATTATGTCGCCCAGATATACTTCGATTTATCACTTTTGCTTGGTTTGCTGTc                                                                                                   >  2:91826/1‑141 (MQ=255)
                                        ttttCAGCATGTTTAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTc                                                                                                   <  2:265290/141‑1 (MQ=255)
                                        ttttCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTACTGTc                                                                                                   <  2:731612/141‑1 (MQ=255)
                                        tttgCATCATGTTGATCTCCCCCAACCTCGTCTGACAACATAAGGTTGTACAAATTAACAAACCAACTTGCAATGGTATGGGACTCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTc                                                                                                   <  2:91629/141‑1 (MQ=255)
                                          ttCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACTAACTTGCAATGGTTGGTGGATCATTATTATGTTGAGCAGGGATTGTGGGATTCATCACTTTTGCTTGGTTTGCTGTCat                                                                                                 >  2:281556/1‑140 (MQ=255)
                                           tCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAAt                                                                                                >  2:220216/1‑141 (MQ=255)
                                           tCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGCCCTGTCAAt                                                                                                >  2:571213/1‑141 (MQ=255)
                                                  gTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTa                                                                                         <  1:821332/141‑1 (MQ=255)
                                                      aaCTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAg                                                                                     <  1:95429/141‑1 (MQ=255)
                                                        ctctCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCCGGAATACTGCGAGTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTaaaaact                                                                                   >  1:784046/1‑136 (MQ=255)
                                                        ctctCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑ct                                                                                  >  2:130079/1‑141 (MQ=255)
                                                         tctcCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAAAAAACCAACGTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑ctc                                                                                 <  1:782104/141‑1 (MQ=255)
                                                                   ttGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCa                                                                       <  2:437225/141‑1 (MQ=255)
                                                                    tgttatCCAAAATAATGTGGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCTATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCAt                                                                      <  2:463598/136‑1 (MQ=255)
                                                                        ttCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCATGTgg                                                                  >  1:357852/1‑141 (MQ=255)
                                                                        ttCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCATGTgg                                                                  <  2:740293/141‑1 (MQ=255)
                                                                         tCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTAGGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAAGACAGGTAATAG‑CTCTTGAAATGCATGTGGa                                                                 >  1:379303/1‑141 (MQ=255)
                                                                            aaaaTAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCATGTGGATCa                                                              >  2:373376/1‑141 (MQ=255)
                                                                            aaaaTAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCATGTGGATCa                                                              <  1:837952/141‑1 (MQ=255)
                                                                                 aatGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATAATTATGCAGCTCACGAATACTGCGATTCATCACCCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑ACCTTGCACTGCAGGTGGATCATaacc                                                         >  2:346339/1‑139 (MQ=255)
                                                                                  atGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCATGTGGATCATAatca                                                        >  2:54068/1‑141 (MQ=255)
                                                                                  atGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATCTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCATGTGGATCCTAatca                                                        <  2:526287/141‑1 (MQ=255)
                                                                                   tgttgtAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAATTGCATGTGGATCATAatcgt                                                       >  2:638333/1‑139 (MQ=255)
                                                                                   tgttgtAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCATGTGGATCATAatcat                                                       <  1:599849/141‑1 (MQ=255)
                                                                                     ttgtAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCATGTGGATCATAAtcatca                                                     <  1:93331/141‑1 (MQ=255)
                                                                                     ttgtAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCATGTGGATCATAAtcatca                                                     >  1:380586/1‑141 (MQ=255)
                                                                                     ttgtAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCATGTGGATCATAAtcatca                                                     >  2:192151/1‑141 (MQ=255)
                                                                                     ttgtAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGAAATAG‑CTCTTGAAATGCATGTGGATCATAAtcatca                                                     >  2:65104/1‑141 (MQ=255)
                                                                                          aaaaTTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCATGTGGATCATAATCATCATGTTg                                                <  2:399648/141‑1 (MQ=255)
                                                                                           aaaTTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTCAAATTACAGGTAATAG‑CTCTTGAAATGCATGTGGATCATAATCATCATGTTGc                                               >  1:649372/1‑141 (MQ=255)
                                                                                            aaatcaCAAACCAACCTGCAATGGAAAGTGAATCATCATTATGTCGCCCAGGAACACCGCCATTCATCACTCCTGCTTGGTGTGCTGTCAATACACGTAATAG‑CTCTTGAAATGCATGTGGATCATAATCATCATgaagca                                              <  2:13011/136‑6 (MQ=255)
                                                                                            aaTTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATATAGGTAATAG‑CTCTTGAAATGCATGTGGATCATAATCATCATGTTGCa                                              <  2:256373/141‑1 (MQ=255)
                                                                                             aTTAACAAACCAACATGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTCGTGTGCTGTCAAGACAGGTAATAG‑CTCTTGAAATGCATGTGGATCATAATCATCATGTTGCat                                             >  2:572695/1‑141 (MQ=255)
                                                                                                    aaCCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCt                                      >  1:572277/1‑141 (MQ=255)
                                                                                                     actaaCTTGCAATGGTAGGTGAATCATCATGATGTCGCACAGGAATACTGCGGTTCATTAGTCTTGGTTGGTTTGCTGGCAATAGGGGTAATAG‑CTATTTCCATGCATGTGGCGCATACTCATCATGTTGccggtctgcga                                     >  2:234758/4‑131 (MQ=255)
                                                                                                           ttGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCATGTGGGTCATAATCGACATGTTGCATATATGCTTTTGCaa                               >  1:189433/1‑141 (MQ=255)
                                                                                                            tGCAATGGTAGGTGAATCATCATTATGTCGCCTAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCTATGGCaaa                              <  1:177457/141‑1 (MQ=255)
                                                                                                            tGCAATGGTAGGTGAATCATCATTATGTCGCCTAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGGCAATACAGGTAATAG‑CTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCTATGGCaaa                              <  1:320870/141‑1 (MQ=255)
                                                                                                               aaTGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGATTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCTATGGCAAAAAc                           <  2:297948/141‑1 (MQ=255)
                                                                                                                atcgTAGGTGAATAATTATTCTGTCGCCCAGGACTATTGCGATTTCTCGCTTTTCCTTGGTGGTCGGGCAATAAAGGTCTCAGACTGTT‑TCTTGCATTTGGATCATAATCATCATGTTGCATATATGCTATGGCAAAAACa                          <  2:804179/138‑1 (MQ=255)
                                                                                                                 tGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGCCACGCATGCGGATCATAATAATCATGTTGCATCTATGCTACGGCAAAAACAg                         >  1:798979/1‑141 (MQ=255)
                                                                                                                 tGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCTATGGCAAAAACAg                         >  1:268742/1‑141 (MQ=255)
                                                                                                                     aGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCTATGGCAAAAACAGATtg                     >  1:558620/1‑141 (MQ=255)
                                                                                                                          aaTCATCATTATGTCGCCCAGGAACACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCTATGGCAAATACAGATTGTGACa                >  1:611474/1‑141 (MQ=255)
                                                                                                                           atcatcATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCTATGGCAAAAACAGTTTGTGACaa               <  1:14303/141‑1 (MQ=255)
                                                                                                                              atcatTATGTGGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCTATGGCAAAAACAGTTTGTGACAATGa            <  2:255320/141‑1 (MQ=255)
                                                                                                                                 atTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCTATGGCAAAAACAGTTTGTGACAATGAttt         >  2:194817/1‑141 (MQ=255)
                                                                                                                                        cGCCTAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCTATGGCAAAAACAGTTTGTGACAATGATTTCTTTTGa  <  1:384780/141‑1 (MQ=255)
                                                                                                                                        cGCACAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCTATGGCAAAAACAGTTTGTGACAATGATTTCTTTTGa  <  1:255370/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                         
GTTTGTGAAATATAATTCCAAGTAACTTTAAATGTTTTGTTTTTCAGCATGTTGAACTCTCCCAAACTTGTCTTCCAAAATAATGTTGTAAAAATTAACAAACCAACTTGCAATGGTAGGTGAATCATCATTATGTCGCCCAGGAATACTGCGATTCATCACTCTTGCTTGGTGTGCTGTCAATACAGGTAATAG‑CTCTTGAAATGCATGTGGATCATAATCATCATGTTGCATATATGCTATGGCAAAAACAGTTTGTGACAATGATTTCTTTTGA  >  CP000730/2807023‑2807299

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