Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 196,832 C→A 100% I806I (ATC→ATA)  bamA → BamABCDE complex OM biogenesis outer membrane pore‑forming assembly factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273196,8320CA99.4% 1058.3 / ‑2.4 329I806I (ATC→ATA) bamABamABCDE complex OM biogenesis outer membrane pore‑forming assembly factor
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (148/178);  total (149/180)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 329 total aligned reads displayed.

CCGGACTATAGTGATCCAAGCAATATCCGTATGTCTGCGGGTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATCGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGCAAATGGGATGGTAAGGAGTTTATTGTGAAAAAGT  >  CP009273/196686‑196978
                                                                                                                                                  |                                                                                                                                                  
   gACTATAGTGATCCAAGCAATATCCGTATGTCTGCGGGTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAATGCAGAACAGTTCCAGTTTAACATAGGTaaa                                                                                                                                              >  1:2551610/1‑150 (MQ=255)
     cTATAGTGATCCAAGCAATATCCGTATGTCTGCGGGTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAAc                                                                                                                                            <  2:1444641/150‑1 (MQ=255)
      tataGTGATCCAAGCAATATCCGTATGTCTGCGGGTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACACGGCAGAACAGTTCCAGTTTAACATAGTTAAAAcc                                                                                                                                           >  2:2892896/1‑150 (MQ=255)
      tataGTGATCCAAGCAATATCCGTATGTCTGCGGGTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAAcc                                                                                                                                           <  1:1715340/150‑1 (MQ=255)
         aGTGATCCAAGCAATATCCGTATGTCTGCGGGTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTgg                                                                                                                                        <  1:499552/150‑1 (MQ=255)
          gTGATCCAAGCAATATCCGTATGTCTGCGGGTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGt                                                                                                                                       <  2:1466561/150‑1 (MQ=255)
           tGATCCAAGCAATATCCGTATGTCTGCGGGTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTa                                                                                                                                      <  1:4936467/150‑1 (MQ=255)
            gATCCAAGCAATATCCGTATGTCTGCGGGTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTaa                                                                                                                                     <  1:4080388/150‑1 (MQ=255)
            gATCCAAGCAATATCCGTATGTCTGCGGGTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTaa                                                                                                                                     <  1:2968931/150‑1 (MQ=255)
             aTCCAAGCAATATCCGTATGTCTGCGGGTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAg                                                                                                                                    <  1:3406109/150‑1 (MQ=255)
             aTCCAAGCAATATCCGTATGTCTGCGGGTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAg                                                                                                                                    <  2:267903/150‑1 (MQ=255)
              tCCAAGCAATATCCGTATGTCTGCGGGTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAgt                                                                                                                                   >  1:886094/1‑150 (MQ=255)
                 aaGCAATATCCGTATGTCTGCGGGTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGtt                                                                                                                                >  2:2608952/1‑150 (MQ=255)
                  aGCAATATCCGTATGTCTGCGGGTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTtc                                                                                                                               <  2:1626974/150‑1 (MQ=255)
                  aGCAATATCCGTATGTCTGCGGGTATCACATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTtc                                                                                                                               >  2:294135/1‑150 (MQ=255)
                      atatCCGTATGTCTGCGGGTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCca                                                                                                                           >  2:2706506/1‑150 (MQ=255)
                       tatCCGTATGTCTGCGGGTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGTGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCcac                                                                                                                          <  1:4450492/150‑1 (MQ=255)
                       tatCCGTATGTCTGCGGGTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCcac                                                                                                                          <  2:1991728/150‑1 (MQ=255)
                        atCCGTATGTCTGCTGTTATCGCATTACAATGGATGTACCCATTGGGGACGTTGGTTTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTTTTCTCcaca                                                                                                                         <  2:1581110/150‑1 (MQ=255)
                        atCCGTATGTCTGCGGGTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGGTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCGGGGAAGTGTTCTCcaca                                                                                                                         >  1:2955534/1‑150 (MQ=255)
                           cGTATGTCTGCGGGTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAg                                                                                                                      <  2:2275903/150‑1 (MQ=255)
                              aTGTCTGCGGGTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGaa                                                                                                                   <  1:3076010/150‑1 (MQ=255)
                               tGTCTGCGGGTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAAt                                                                                                                  <  2:4668426/150‑1 (MQ=255)
                               tGTCTGCGGGTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAAt                                                                                                                  <  2:2638560/150‑1 (MQ=255)
                                gTCTGCGGGTATCGATTTTCATTGGATGTCCCCATTGGGGCAGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAAGAGTTCCAGTTTTACATAGGTAAAACCTGGTTAGTGTTCTCCAGAATGGATTg                                                                                                                 <  2:5036903/150‑1 (MQ=255)
                                     cGGGTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATgtagtg                                                                                                            <  1:2041797/150‑1 (MQ=255)
                                       ggTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATgtagtggt                                                                                                          <  1:4483937/150‑1 (MQ=255)
                                        gTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATCGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATgtagtggta                                                                                                         <  1:1373973/150‑1 (MQ=255)
                                        gTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATgtagtggta                                                                                                         <  2:397231/150‑1 (MQ=255)
                                        gTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATgtagtggta                                                                                                         >  2:843475/1‑150 (MQ=255)
                                        gTATCGCATTACAATGGATGTCCCAATTGGTGCCGTTGGTGTTATCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTAACACTTTAACATAGTTAAAAACTGGTAAGTGTTCTACACCCAGGAATgtagaggta                                                                                                         >  2:3756426/1‑150 (MQ=255)
                                         tATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATgtagtggtag                                                                                                        <  2:3877561/150‑1 (MQ=255)
                                              cATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAAACTGGTAAGTGTTCTCCACAAAGGAATGTAGGGGTAGTGTAg                                                                                                   >  1:2004367/1‑150 (MQ=255)
                                              cATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGGTCAAAAAGTACGGTGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCGGGTAAGTGTTATCCACAAAGGAATGTAGTGGTAGTGTAg                                                                                                   >  1:4345354/1‑150 (MQ=255)
                                                  aCAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGAt                                                                                               <  1:2404971/150‑1 (MQ=255)
                                                   cAATGGAGGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATg                                                                                              <  2:5100418/150‑1 (MQ=255)
                                                    aaTGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGa                                                                                             >  1:4316258/1‑150 (MQ=255)
                                                          tGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAg                                                                                       >  1:1319979/1‑150 (MQ=255)
                                                           gTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAgg                                                                                      <  2:886094/150‑1 (MQ=255)
                                                           gTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAgg                                                                                      <  1:5160961/150‑1 (MQ=255)
                                                             ccccATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCg                                                                                    >  1:2671894/1‑150 (MQ=255)
                                                              cccATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGa                                                                                   >  2:1591128/1‑150 (MQ=255)
                                                               ccATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGAt                                                                                  <  1:466618/150‑1 (MQ=255)
                                                                cATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATc                                                                                 >  2:2092513/1‑150 (MQ=255)
                                                                 aTTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCa                                                                                >  2:523121/1‑150 (MQ=255)
                                                                 aTTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCa                                                                                <  1:5272001/150‑1 (MQ=255)
                                                                  ttGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAATGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACCAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCaa                                                                               >  1:1788615/1‑150 (MQ=255)
                                                                  ttGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCaa                                                                               <  2:4178575/150‑1 (MQ=255)
                                                                  ttGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCaa                                                                               <  2:5064995/150‑1 (MQ=255)
                                                                        ccGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCTATATAAg                                                                         <  1:752697/150‑1 (MQ=255)
                                                                        ccGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAg                                                                         <  2:1501723/150‑1 (MQ=255)
                                                                          gTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGAt                                                                       >  2:1158336/1‑150 (MQ=255)
                                                                           ttGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATc                                                                      <  2:2823210/150‑1 (MQ=255)
                                                                           ttGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATc                                                                      <  1:2442857/150‑1 (MQ=255)
                                                                              gtgtTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAACCCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGcc                                                                   <  1:2912009/150‑1 (MQ=255)
                                                                               tgtTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCg                                                                  >  2:4490197/1‑150 (MQ=255)
                                                                                    tcCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCa                                                             <  2:109025/150‑1 (MQ=255)
                                                                                          gCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCaaa                                                       >  2:48356/1‑150 (MQ=255)
                                                                                            ccAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGa                                                     >  1:2648652/1‑150 (MQ=255)
                                                                                              aGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTTTTCTCCACAAAGGAATGTAGTGGTAGTGTACCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAAc                                                   <  1:4230412/150‑1 (MQ=255)
                                                                                              aGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAAc                                                   <  2:549439/150‑1 (MQ=255)
                                                                                              aGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAAc                                                   <  1:3854038/150‑1 (MQ=255)
                                                                                               gCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACt                                                  <  1:3116298/150‑1 (MQ=255)
                                                                                                ccGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCATTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTg                                                 <  1:2301749/150‑1 (MQ=255)
                                                                                                 cGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCACAGAACTGc                                                >  2:578916/1‑150 (MQ=255)
                                                                                                     cAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCAccc                                            >  2:4117868/1‑150 (MQ=255)
                                                                                                      aaaaaGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCt                                           >  2:3287381/1‑150 (MQ=255)
                                                                                                        aaaGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTcc                                         <  1:3804087/150‑1 (MQ=255)
                                                                                                         aaGTACGCTGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTTGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCg                                        <  1:1897964/150‑1 (MQ=255)
                                                                                                         aaGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCg                                        <  1:1878602/150‑1 (MQ=255)
                                                                                                          aGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCgg                                       >  2:2399541/1‑150 (MQ=255)
                                                                                                             aCGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAACCCTGGTAAGTGTTCTCCACACAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGc                                    <  1:3626660/150‑1 (MQ=255)
                                                                                                             aCGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGc                                    <  1:4788312/150‑1 (MQ=255)
                                                                                                               gATGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGCaa                                  <  1:1167905/150‑1 (MQ=255)
                                                                                                                aTGGAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGCaaa                                 <  2:4315749/150‑1 (MQ=255)
                                                                                                                  ggAGACAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGCAAATg                               >  1:2993261/1‑150 (MQ=255)
                                                                                                                   gagaCAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGCAAATgg                              <  1:3899115/150‑1 (MQ=255)
                                                                                                                     gaCAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGCAAATGGGa                            >  1:1589570/1‑150 (MQ=255)
                                                                                                                     gaCAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGCAAATGGGa                            <  1:5057435/150‑1 (MQ=255)
                                                                                                                      aCAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCCAAGACCTGCACCCTCCGGTGCACATGGGAt                           >  1:46240/1‑150 (MQ=255)
                                                                                                                       cAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTTCACCCTCCGGTGCAAATGGGATg                          <  2:4738000/150‑1 (MQ=255)
                                                                                                                       cAAGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGCAAATGGGATg                          <  1:5174231/150‑1 (MQ=255)
                                                                                                                        aaGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGCAAATGGGATgg                         <  1:663384/150‑1 (MQ=255)
                                                                                                                        aaGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGCAAATGGGATgg                         <  1:1264913/150‑1 (MQ=255)
                                                                                                                         aGGCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGCAAATGGGATGGt                        >  1:3339834/1‑150 (MQ=255)
                                                                                                                          ggCAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGCAAATGGGATGGTa                       >  2:1237585/1‑150 (MQ=255)
                                                                                                                            cAGAACAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGCAAATGGGATGGTAAg                     <  2:2502353/150‑1 (MQ=255)
                                                                                                                               aaCAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGCAAATGGGATGGTAAGTAg                  <  2:1313709/150‑1 (MQ=255)
                                                                                                                                aCAGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGCAAATGGGATGGTAAGGAGt                 >  1:2374097/1‑150 (MQ=255)
                                                                                                                                  aGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACACAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGCAAATGGGATGGTAAGGAGttt               <  1:2798143/150‑1 (MQ=255)
                                                                                                                                  aGTTCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGCAAATGGGATGGTAAGGAGttt               <  2:3069929/150‑1 (MQ=255)
                                                                                                                                     tacAGTTTAACATAGGTAAAACCTGGTTAGTGTTCTCCACATAGGATTGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATACGATCGCCGGGCCTCGCAAAGAACTGCACCCTCCGGTGCAAATGGGATGGTAAGGAGTTTAtt            <  1:197391/148‑1 (MQ=255)
                                                                                                                                     tCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGCAAATGGGATGGTAAGGAGTTTAtt            <  2:443318/150‑1 (MQ=255)
                                                                                                                                     tCCAGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGCAAATGGGATGGTAAGGAGTTTAtt            >  2:3160766/1‑150 (MQ=255)
                                                                                                                                        aGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGCAAATGGGATGGTAAGGAGTTTATtgtg         >  2:2197505/1‑150 (MQ=255)
                                                                                                                                        aGTTTAACATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGCAAATGGGATGGTAAGGAGTTTATtgtg         >  2:1198447/1‑150 (MQ=255)
                                                                                                                                             aaCATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGCAAATGGGATGGTAAGGAGTTTATTGTGaaaaa    <  2:770520/150‑1 (MQ=255)
                                                                                                                                               cATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATTTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGCAAATGGGATGGTAAGGAGTTTATTGTGAAAAAGt  <  2:3475738/150‑1 (MQ=255)
                                                                                                                                               cATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGCAAATGGGATGGTAAGGAGTTTATTGTGAAAAAGt  >  1:1510210/1‑150 (MQ=255)
                                                                                                                                               cATAGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGCAAATGGGATGGTAAGGAGTTTATTGTGAAAAAGt  >  1:1413946/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CCGGACTATAGTGATCCAAGCAATATCCGTATGTCTGCGGGTATCGCATTACAATGGATGTCCCCATTGGGGCCGTTGGTGTTCTCCTACGCCCAGCCGTTCAAAAAGTACGATGGAGACAAGGCAGAACAGTTCCAGTTTAACATCGGTAAAACCTGGTAAGTGTTCTCCACAAAGGAATGTAGTGGTAGTGTAGCGATGACTTTAGGCGATCAATATAAGATCGCCGGGCCACGCAAAGAACTGCACCCTCCGGTGCAAATGGGATGGTAAGGAGTTTATTGTGAAAAAGT  >  CP009273/196686‑196978

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