Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,506,430 A→T 19.0% L22Q (CTG→CAG)  fecA ← ferric citrate outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,506,4300AT19.0% 688.2 / 150.7 258L22Q (CTG→CAG) fecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (88/121);  new base T (18/31);  total (106/152)
Fisher's exact test for biased strand distribution p-value = 5.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 268 total aligned reads displayed.

CGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATT  >  CP009273/4506281‑4506576
                                                                                                                                                     |                                                                                                                                                  
cGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGgca                                                                                                                                                    >  1:3518010/1‑150 (MQ=255)
 gTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGgcag                                                                                                                                                   >  1:1999500/1‑150 (MQ=255)
   tGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGgcagca                                                                                                                                                 <  2:4478441/150‑1 (MQ=255)
    gCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGcagcag                                                                                                                                                <  1:163734/150‑1 (MQ=255)
      tctGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCTGCAGGc                                                                                                                                              <  1:735259/150‑1 (MQ=255)
         gCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCa                                                                                                                                           <  2:4797787/150‑1 (MQ=255)
         gCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCa                                                                                                                                           <  1:3608871/150‑1 (MQ=255)
          cTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAg                                                                                                                                          <  2:752722/150‑1 (MQ=255)
            tGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGc                                                                                                                                        <  1:1407723/150‑1 (MQ=255)
                gcgcgTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAAt                                                                                                                                    <  1:109155/150‑1 (MQ=255)
                gcgcgTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAAt                                                                                                                                    <  1:2561162/150‑1 (MQ=255)
                gcgcgTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCGGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCTGCAGGCTCAGGCGAAt                                                                                                                                    <  1:2865436/150‑1 (MQ=255)
                 cgcgTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCTGCAGGCTCAGGCGAATg                                                                                                                                   <  2:1665088/150‑1 (MQ=255)
                   cgTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGgt                                                                                                                                 >  1:2389158/1‑150 (MQ=255)
                    gTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGgtg                                                                                                                                <  1:4845626/150‑1 (MQ=255)
                     tCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGgtgt                                                                                                                               <  2:4295294/150‑1 (MQ=255)
                      cAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGtt                                                                                                                              <  1:1317869/150‑1 (MQ=255)
                       aGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTa                                                                                                                             <  2:4369626/150‑1 (MQ=255)
                         gCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAAc                                                                                                                           <  1:3967684/150‑1 (MQ=255)
                          cTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAAcc                                                                                                                          >  2:3010421/1‑150 (MQ=255)
                             gCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCTGCAGGCTCAGGCGAATGGTGTTAACCaaa                                                                                                                       <  2:637470/150‑1 (MQ=255)
                             gCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCaaa                                                                                                                       <  2:4053173/150‑1 (MQ=255)
                                 cAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCTGCAGGCTCAGGCGAATGGTGTTAACCAAAGgtg                                                                                                                   <  2:3552387/150‑1 (MQ=255)
                                   aCCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGtgtt                                                                                                                 <  2:267485/150‑1 (MQ=255)
                                      gagaGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTgttgtt                                                                                                              >  1:4414715/1‑150 (MQ=255)
                                        gagGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGtttt                                                                                                            <  2:3136364/150‑1 (MQ=255)
                                         agGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTa                                                                                                           >  1:242799/1‑150 (MQ=255)
                                           ggTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCTGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACg                                                                                                         <  1:4731285/150‑1 (MQ=255)
                                           ggTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCTGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACg                                                                                                         <  1:2280586/150‑1 (MQ=255)
                                              aaaTCCGCTGTGTGCTGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGaaa                                                                                                      <  2:3500235/150‑1 (MQ=255)
                                              aaaTCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGaaa                                                                                                      <  1:2015464/150‑1 (MQ=255)
                                               aaTCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGaaaa                                                                                                     >  2:4243029/1‑150 (MQ=255)
                                                aTCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGaaaaa                                                                                                    <  1:4531899/150‑1 (MQ=255)
                                                aTCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGaaaaa                                                                                                    <  1:1012466/150‑1 (MQ=255)
                                                aTCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGaaaaa                                                                                                    >  2:882077/1‑150 (MQ=255)
                                                aTCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGaaaaa                                                                                                    >  1:1526273/1‑150 (MQ=255)
                                                  ccGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCTGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAAcg                                                                                                  >  2:2446580/1‑150 (MQ=255)
                                                     cTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGt                                                                                               >  1:1737237/1‑150 (MQ=255)
                                                     cTGTGTGCGGCATACTGAGTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGt                                                                                               <  1:4474593/150‑1 (MQ=255)
                                                      tgtgtgCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTa                                                                                              >  1:397640/1‑150 (MQ=255)
                                                      tgtgtgCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTa                                                                                              >  1:1726048/1‑150 (MQ=255)
                                                       gtgtgCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCTGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTaa                                                                                             >  1:2993349/1‑150 (MQ=255)
                                                        tgtgCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAAc                                                                                            >  1:347448/1‑150 (MQ=255)
                                                         gtgCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACg                                                                                           >  1:2496909/1‑150 (MQ=255)
                                                             ggCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGt                                                                                       >  2:2855868/1‑150 (MQ=255)
                                                               cATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCTGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCa                                                                                     <  1:3581128/150‑1 (MQ=255)
                                                               cATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCa                                                                                     <  1:3343745/150‑1 (MQ=255)
                                                               cATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCa                                                                                     <  1:903866/150‑1 (MQ=255)
                                                                aTACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCAt                                                                                    >  2:4872163/1‑150 (MQ=255)
                                                                aTACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCAt                                                                                    <  2:1907559/150‑1 (MQ=255)
                                                                   cTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCTGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATAcc                                                                                 <  2:2785836/150‑1 (MQ=255)
                                                                    tGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCTGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCt                                                                                >  1:1233524/1‑150 (MQ=255)
                                                                    tGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCt                                                                                >  1:4697869/1‑150 (MQ=255)
                                                                     gATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCtt                                                                               >  1:4288583/1‑150 (MQ=255)
                                                                      aTTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCTGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTc                                                                              <  1:756017/150‑1 (MQ=255)
                                                                        tGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTccc                                                                            <  2:3071284/150‑1 (MQ=255)
                                                                          aGCGCTTTATCGAGCGATCCCGGTGCGATATTAGCCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCca                                                                          >  2:4575723/1‑150 (MQ=255)
                                                                           gcgcTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCcat                                                                         >  2:3889450/1‑150 (MQ=255)
                                                                                 tATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCcatcatcat                                                                   <  1:3836150/150‑1 (MQ=255)
                                                                                  aTCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCAtt                                                                  <  2:1231057/150‑1 (MQ=255)
                                                                                  aTCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCAtt                                                                  >  1:4177205/1‑150 (MQ=255)
                                                                                    cGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCAtttt                                                                >  2:2216221/1‑150 (MQ=255)
                                                                                     gAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCAttttt                                                               >  2:3481592/1‑150 (MQ=255)
                                                                                       gCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTttgt                                                             >  2:3978823/1‑150 (MQ=255)
                                                                                         gATCCCGGTGCGATATTCACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAACACGCGTAACGGCGTCATACCTTCCCCATCATCATTTttgttg                                                           >  2:2907839/1‑150 (MQ=255)
                                                                                             ccGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGtt                                                       >  1:3109515/1‑150 (MQ=255)
                                                                                               ggTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCa                                                     >  1:693443/1‑150 (MQ=255)
                                                                                               ggTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCa                                                     <  1:1911815/150‑1 (MQ=255)
                                                                                                        ttAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGt                                            <  1:4670956/150‑1 (MQ=255)
                                                                                                         tAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCAAATCATCATTTTTGTTGTGTTCAGATATGAGTc                                           >  2:3976889/1‑150 (MQ=255)
                                                                                                          aaCCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCTGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCg                                          >  1:965847/1‑150 (MQ=255)
                                                                                                           aCCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGa                                         <  2:1257430/150‑1 (MQ=255)
                                                                                                            ccTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGaa                                        <  1:1037510/150‑1 (MQ=255)
                                                                                                             cTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAAc                                       >  2:2847963/1‑150 (MQ=255)
                                                                                                             cTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAAc                                       <  2:672549/150‑1 (MQ=255)
                                                                                                              tgtgCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGaaag                                      <  1:2611907/150‑3 (MQ=255)
                                                                                                               gtgCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACga                                     <  2:58174/150‑1 (MQ=255)
                                                                                                                 gcagcaAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACgaga                                   <  1:2547501/150‑1 (MQ=255)
                                                                                                                  cagcaACAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGaa                                  <  2:190282/150‑1 (MQ=255)
                                                                                                                  cagcaAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGaa                                  >  2:14052/1‑150 (MQ=255)
                                                                                                                   agcaAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAg                                 <  2:3932705/150‑1 (MQ=255)
                                                                                                                     caAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAggg                               >  2:150951/1‑150 (MQ=255)
                                                                                                                      aaaaGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGt                              >  1:522674/1‑150 (MQ=255)
                                                                                                                       aaaGCGGAAAACGAGAGACCGGCCAGCGGCTGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTa                             <  2:3785385/150‑1 (MQ=255)
                                                                                                                         aGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTaaa                           >  1:1768978/1‑150 (MQ=255)
                                                                                                                            ggAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAg                        >  2:2414709/1‑150 (MQ=255)
                                                                                                                             gAAAACGAGAGACCGGCCAGCGGCTGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAgg                       >  1:4002894/1‑150 (MQ=255)
                                                                                                                             gAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAgg                       >  1:4777379/1‑150 (MQ=255)
                                                                                                                             gAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAgg                       >  1:3701961/1‑150 (MQ=255)
                                                                                                                               aaaCGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGAc                     <  2:2023609/150‑1 (MQ=255)
                                                                                                                               aaaCGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGAc                     <  2:3155374/150‑1 (MQ=255)
                                                                                                                                aaCGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACa                    <  2:155089/150‑1 (MQ=255)
                                                                                                                                aaCGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACa                    <  2:3223308/150‑1 (MQ=255)
                                                                                                                                   gagagaCCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATc                 >  1:649196/1‑150 (MQ=255)
                                                                                                                                    agagaCCGGCCAGCGGCAGCAGGCTCAGTCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAATGGTAAAAAGGACAATCg                <  2:3393572/150‑1 (MQ=255)
                                                                                                                                            gCCAGCGGCTGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGa        >  2:472925/1‑150 (MQ=255)
                                                                                                                                            gCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGa        <  1:722274/150‑1 (MQ=255)
                                                                                                                                               aGCGGCTGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAatt     <  2:3210264/150‑1 (MQ=255)
                                                                                                                                                 cGGCTGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAattat   <  1:599878/150‑1 (MQ=255)
                                                                                                                                                 cGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAattat   >  2:230576/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATT  >  CP009273/4506281‑4506576

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: