Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,477,840 G→T 6.0% D175Y (GAC→TAC)  hrpA → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,477,8400GT 6.0% 958.8 / 28.2 249D175Y (GAC→TAC) hrpAputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (116/118);  new base T (11/4);  total (127/122)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 255 total aligned reads displayed.

TCGGCCATACCCAGCCGCGTCGTCTGGCGGCAAGAACAGTGGCGAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAA  >  CP009273/1477691‑1477989
                                                                                                                                                     |                                                                                                                                                     
tCGGCCATACCCAGCCGCGTCGTCTGGCGGCAAGAACAGTGGCGAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCg                                                                                                                                                       <  2:3342739/150‑1 (MQ=255)
  ggCCATACCCAGCCGCGTCGTCTGGCGGCAAGAACAGTGGCGAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGAc                                                                                                                                                     <  2:4041997/150‑1 (MQ=255)
   gCCATACCCAGCCGCGTCGTCTGGCGGCAAGAACAGTGGCGAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCTACg                                                                                                                                                    <  1:662823/150‑1 (MQ=255)
      aTACCCAGCCGCGTCGTCTGGCGGCAAGAACAGTGGCGAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTa                                                                                                                                                 >  1:3803500/1‑150 (MQ=255)
      aTACCCAGCCGCGTCGTCTGGCGGCAAGAACAGTGGCGAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTa                                                                                                                                                 >  2:4319276/1‑150 (MQ=255)
      aTACCCAGCCGCGTCGTCTGGCGGCAAGAACAGTGGCGAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGAACACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTa                                                                                                                                                 >  1:1505427/1‑150 (MQ=255)
            aGCCGCGTCGTCTGGCGGCAAGAACAGTGGCGAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCctgc                                                                                                                                           >  1:2249956/1‑150 (MQ=255)
              ccGCGTCGTCTGGCGGCAAGAACAGTGGCGAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCctgctg                                                                                                                                         >  1:4112038/1‑150 (MQ=255)
              ccGCGTCGTCTGGCGGCAAGAACAGTGGCGAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCctgctg                                                                                                                                         <  1:585819/150‑1 (MQ=255)
                 cgtcgtCTGGCGGCAAGAACAGTGGCGAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCg                                                                                                                                      <  1:1096947/150‑1 (MQ=255)
                   tcgtcTGGCGGCAAGAACAGTGGCGAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGga                                                                                                                                    <  2:4426796/150‑1 (MQ=255)
                   tcgtcTGGCGGCAAGAACAGTGGCGAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGga                                                                                                                                    <  2:4076775/150‑1 (MQ=255)
                    cgtcTGGCGGCAAGAACAGTGGCGAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGgag                                                                                                                                   >  1:163298/1‑150 (MQ=255)
                     gtcTGGCGGCAAGAACAGTGGCGAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGAGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGgaga                                                                                                                                  >  1:4013391/1‑150 (MQ=255)
                        tGGCGGCAAGAACAGTGGCGAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATcc                                                                                                                               <  1:2561122/150‑1 (MQ=255)
                         ggcggcAAGAACAGTGGCGAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCa                                                                                                                              >  1:262601/1‑150 (MQ=255)
                              cAAGAACAGTGGCGAACCGAATTGCGGACGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAg                                                                                                                         <  1:3376665/150‑1 (MQ=255)
                                aGAACAGTGGCGAACCTTATTGCGTAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGAc                                                                                                                       <  1:1134998/150‑1 (MQ=255)
                                 gAACAGTGGCGAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGAcc                                                                                                                      <  2:3292640/150‑1 (MQ=255)
                                 gAACAGTGGCGAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGAcc                                                                                                                      <  1:3613763/150‑1 (MQ=255)
                                  aaCAGTGGCGAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCg                                                                                                                     >  1:3092254/1‑150 (MQ=255)
                                    cAGTGGCGAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGcc                                                                                                                   >  2:2419047/1‑150 (MQ=255)
                                        ggCGAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCctgct                                                                                                               >  1:3245220/1‑150 (MQ=255)
                                           gAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCTACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGAt                                                                                                            >  1:480691/1‑150 (MQ=255)
                                           gAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGAt                                                                                                            >  2:2085927/1‑150 (MQ=255)
                                             aCCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGc                                                                                                          >  1:3002527/1‑150 (MQ=255)
                                              ccGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCa                                                                                                         <  1:1812419/150‑1 (MQ=255)
                                               cGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAg                                                                                                        >  1:1847913/1‑150 (MQ=255)
                                               cGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAg                                                                                                        <  2:158910/150‑1 (MQ=255)
                                                gTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGt                                                                                                       >  1:283812/1‑150 (MQ=255)
                                                gTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGt                                                                                                       >  1:284531/1‑150 (MQ=255)
                                                 tATTGCGGAAGAGATGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTa                                                                                                      <  1:1413418/150‑1 (MQ=255)
                                                        gAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCTACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACt                                                                                               >  1:3499174/1‑150 (MQ=255)
                                                         aaGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTa                                                                                              <  2:2055882/150‑1 (MQ=255)
                                                             gCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCat                                                                                          <  1:537483/150‑1 (MQ=255)
                                                              cTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCatt                                                                                         <  2:2603385/150‑1 (MQ=255)
                                                               tGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCatta                                                                                        <  1:1353145/150‑1 (MQ=255)
                                                                gAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCATCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCattat                                                                                       >  1:3664209/1‑150 (MQ=255)
                                                                gAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCattat                                                                                       <  2:2237383/150‑1 (MQ=255)
                                                                gAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCattat                                                                                       <  2:897331/150‑1 (MQ=255)
                                                                  aaaCGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTg                                                                                     >  2:839300/1‑150 (MQ=255)
                                                                           cGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAgcgc                                                                            >  1:4524461/1‑150 (MQ=255)
                                                                           cGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAgcgc                                                                            <  1:2633021/150‑1 (MQ=255)
                                                                             ggCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCAc                                                                          <  1:100064/150‑1 (MQ=255)
                                                                             ggCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCAc                                                                          <  1:2008385/150‑1 (MQ=255)
                                                                              gCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACg                                                                         <  2:2848257/150‑1 (MQ=255)
                                                                                 gTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAAc                                                                      <  1:381634/150‑1 (MQ=255)
                                                                                   tGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGc                                                                    >  1:1054037/1‑150 (MQ=255)
                                                                                    gCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCTACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCa                                                                   >  1:977310/1‑150 (MQ=255)
                                                                                      aTCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGc                                                                 >  2:411448/1‑150 (MQ=255)
                                                                                       tCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGcc                                                                <  2:2725388/150‑1 (MQ=255)
                                                                                         ggTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTg                                                              <  2:305788/150‑1 (MQ=255)
                                                                                         ggTTACAAAGTGCGTTTCAGAGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTg                                                              <  1:3835401/150‑1 (MQ=255)
                                                                                          gTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCTACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGa                                                             >  2:3403439/1‑150 (MQ=255)
                                                                                             aCAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAata                                                          >  1:3571205/1‑150 (MQ=255)
                                                                                               aaaGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATc                                                        >  1:843911/1‑150 (MQ=255)
                                                                                                aaGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCg                                                       <  2:3857283/150‑1 (MQ=255)
                                                                                                aaGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCg                                                       <  1:3858440/150‑1 (MQ=255)
                                                                                                aaGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCg                                                       <  2:1380923/150‑1 (MQ=255)
                                                                                                   tGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGAtt                                                    <  2:1814304/150‑1 (MQ=255)
                                                                                                     cGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGAtttt                                                  >  2:2840052/1‑150 (MQ=255)
                                                                                                     cGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGAtttt                                                  <  2:711683/150‑1 (MQ=255)
                                                                                                      gTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGAttttt                                                 >  2:631291/1‑150 (MQ=255)
                                                                                                       tttCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGAtttttt                                                <  1:3920136/150‑1 (MQ=255)
                                                                                                         tCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGc                                              >  2:3590375/1‑150 (MQ=255)
                                                                                                         tCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGc                                              >  2:1897079/1‑150 (MQ=255)
                                                                                                           aGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTc                                            >  2:4352088/1‑150 (MQ=255)
                                                                                                           aGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTc                                            <  2:2421797/150‑1 (MQ=255)
                                                                                                            gCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCg                                           <  2:123079/150‑1 (MQ=255)
                                                                                                            gCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCg                                           <  2:3201463/150‑1 (MQ=255)
                                                                                                              gATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGc                                         >  2:1245942/1‑150 (MQ=255)
                                                                                                               aTCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCt                                        >  1:1382705/1‑150 (MQ=255)
                                                                                                               aTCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCt                                        >  2:1976293/1‑150 (MQ=255)
                                                                                                               aTCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCt                                        <  2:3121108/150‑1 (MQ=255)
                                                                                                                   cGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTAttt                                    >  1:440769/1‑150 (MQ=255)
                                                                                                                   cGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTAttt                                    >  1:1775511/1‑150 (MQ=255)
                                                                                                                   cGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTAttt                                    <  2:2593100/150‑1 (MQ=255)
                                                                                                                   cGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTAttt                                    <  1:796643/150‑1 (MQ=255)
                                                                                                                         tGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAaga                              >  1:3551775/1‑150 (MQ=255)
                                                                                                                          gATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAagag                             <  2:2748648/150‑1 (MQ=255)
                                                                                                                          gATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAagag                             >  1:243831/1‑150 (MQ=255)
                                                                                                                           aTAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGt                            >  1:374/1‑150 (MQ=255)
                                                                                                                           aTAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGt                            <  1:3313376/150‑1 (MQ=255)
                                                                                                                              acacGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTtgc                         <  1:2709292/150‑1 (MQ=255)
                                                                                                                               cacGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTtgct                        <  1:2238148/150‑1 (MQ=255)
                                                                                                                                acGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTtgctg                       <  2:4193222/150‑1 (MQ=255)
                                                                                                                                 cGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTtgctgc                      >  1:423551/1‑150 (MQ=255)
                                                                                                                                      gTCAAGCTGATGACCTACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCgcgg                 >  1:2699205/1‑150 (MQ=255)
                                                                                                                                      gTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCgcgg                 >  1:2278773/1‑150 (MQ=255)
                                                                                                                                      gTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCgcgg                 >  1:13098/1‑150 (MQ=255)
                                                                                                                                        cAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCgcggcg               >  1:3836620/1‑150 (MQ=255)
                                                                                                                                        cAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCgcggcg               >  1:2170185/1‑150 (MQ=255)
                                                                                                                                        cAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCgcggcg               >  1:632378/1‑150 (MQ=255)
                                                                                                                                         aaGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGt              >  2:3184919/1‑150 (MQ=255)
                                                                                                                                          aGCTGATGACCTACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTc             >  1:1550134/1‑150 (MQ=255)
                                                                                                                                          aGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTc             >  1:1912519/1‑150 (MQ=255)
                                                                                                                                          aGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTc             >  1:4189919/1‑150 (MQ=255)
                                                                                                                                             tgatgaCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTg          <  2:1671385/150‑1 (MQ=255)
                                                                                                                                               atgaCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGAc        <  2:3635568/150‑1 (MQ=255)
                                                                                                                                                     gACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTaaaa  >  1:4530687/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGGCCATACCCAGCCGCGTCGTCTGGCGGCAAGAACAGTGGCGAACCGTATTGCGGAAGAGCTGAAAACGGAGCCGGGCGGTTGCATCGGTTACAAAGTGCGTTTCAGCGATCACGTAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAA  >  CP009273/1477691‑1477989

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