Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,020,228 G→T 23.1% A368A (GCG→GCT guaD → guanine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,020,2280GT23.1% 840.6 / 216.9 333A368A (GCG→GCTguaDguanine deaminase
Reads supporting (aligned to +/- strand):  ref base G (111/145);  new base T (26/51);  total (137/196)
Fisher's exact test for biased strand distribution p-value = 1.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 338 total aligned reads displayed.

AAAGTTAAAGTGGGCATGGGAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGAC  >  CP009273/3020079‑3020377
                                                                                                                                                     |                                                                                                                                                     
aaaGTTAAAGTGGGCATGGGAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCg                                                                                                                                                       >  1:718665/1‑150 (MQ=255)
 aaGTTAAAGTGGGCATGGGAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGa                                                                                                                                                      >  2:756837/1‑150 (MQ=255)
 aaGTTAAAGTGGGCATGGGAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGa                                                                                                                                                      >  2:2944108/1‑150 (MQ=255)
  aGTTAAAGTGGGCATGGGAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGaa                                                                                                                                                     >  1:2443371/1‑150 (MQ=255)
   gTTAAAGTGGGCATGGGAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGaaa                                                                                                                                                    <  1:406739/150‑1 (MQ=255)
   gTTAAAGTGGGCATGGGAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGaaa                                                                                                                                                    >  2:1218525/1‑150 (MQ=255)
    ttAAAGTGGGCATGGGAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAAt                                                                                                                                                   <  1:615557/150‑1 (MQ=255)
    ttAAAGTGGGCATGGGAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAAt                                                                                                                                                   <  1:3045800/150‑1 (MQ=255)
    ttAAAGTGGGCATGGGAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAAt                                                                                                                                                   <  2:1061619/150‑1 (MQ=255)
      aaaGTGGGCATGGGAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAAtct                                                                                                                                                 >  2:4320969/1‑150 (MQ=255)
      aaaGTGGGCATGGGAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAAtct                                                                                                                                                 >  2:1867804/1‑150 (MQ=255)
      aaaGTGGGCATGGGAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAAtct                                                                                                                                                 <  2:2029385/150‑1 (MQ=255)
        aGTGGGCATGGGAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCTAAATctct                                                                                                                                               <  1:4646220/150‑1 (MQ=255)
        aGTGGGCATGGGAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCTAAATctct                                                                                                                                               >  2:3581289/1‑150 (MQ=255)
         gTGGGCATGGGAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTg                                                                                                                                              >  2:644494/1‑150 (MQ=255)
               aTGGGAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCtt                                                                                                                                        <  2:2252259/150‑1 (MQ=255)
               aTGGGAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCtt                                                                                                                                        >  2:3039645/1‑150 (MQ=255)
               aTGGGAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCtt                                                                                                                                        >  2:4807247/1‑150 (MQ=255)
                tGGGAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCTAAATCTCTGGGCCTTg                                                                                                                                       <  1:3144864/150‑1 (MQ=255)
                 gggAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGa                                                                                                                                      <  1:4802256/150‑1 (MQ=255)
                  ggAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCTAAATCTCTGGGCCTTGAc                                                                                                                                     <  1:2036082/150‑1 (MQ=255)
                  ggAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGAc                                                                                                                                     <  1:885618/150‑1 (MQ=255)
                  ggAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCATACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGAc                                                                                                                                     <  2:4512602/150‑1 (MQ=255)
                   gAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCTAAATCTCTGGGCCTTGACg                                                                                                                                    <  1:3220522/150‑1 (MQ=255)
                     aCGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGAt                                                                                                                                  >  1:4515183/1‑150 (MQ=255)
                      cGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCTAAATCTCTGGGCCTTGACGAtt                                                                                                                                 <  1:4690239/150‑1 (MQ=255)
                       ggATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGAttt                                                                                                                                <  2:3974997/150‑1 (MQ=255)
                          tatCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGAt                                                                                                                             <  1:8440/150‑1 (MQ=255)
                           atCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGAtt                                                                                                                            <  1:2478855/150‑1 (MQ=255)
                            tCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTg                                                                                                                           <  1:1474879/150‑1 (MQ=255)
                             cGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCTAAATCTCTGGGCCTTGACGATTTGATTgg                                                                                                                          >  1:4756094/1‑150 (MQ=255)
                               gTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCa                                                                                                                        <  1:1244561/150‑1 (MQ=255)
                                  ccGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCTAAATCTCTGGGCCTTGACGATTTGATTGGCAACt                                                                                                                     <  2:871192/150‑1 (MQ=255)
                                  ccGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCTAAATCTCTGGGCCTTGACGATTTGATTGGCAACt                                                                                                                     <  1:1771265/150‑1 (MQ=255)
                                   cGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACtt                                                                                                                    >  2:1440157/1‑150 (MQ=255)
                                     gAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACtttt                                                                                                                  <  1:1156022/150‑1 (MQ=255)
                                      aaCCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACttttt                                                                                                                 >  1:4290032/1‑150 (MQ=255)
                                          aCTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCt                                                                                                             <  1:216467/150‑1 (MQ=255)
                                              tCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCa                                                                                                         >  1:1517334/1‑150 (MQ=255)
                                              tCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCa                                                                                                         <  2:4009491/150‑1 (MQ=255)
                                                   aTGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCTAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAagag                                                                                                    <  2:515058/150‑1 (MQ=255)
                                                   aTGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAagag                                                                                                    <  1:2305073/150‑1 (MQ=255)
                                                      ctgcAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCt                                                                                                 <  1:1866185/150‑1 (MQ=255)
                                                       tgcAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCTAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTg                                                                                                <  1:4744905/150‑1 (MQ=255)
                                                        gcAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGa                                                                                               >  2:39492/1‑150 (MQ=255)
                                                         cAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCTAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGAt                                                                                              >  2:1608547/1‑150 (MQ=255)
                                                           aaCGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGAttt                                                                                            >  1:4147788/1‑150 (MQ=255)
                                                           aaCGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGAttt                                                                                            <  1:1914256/150‑1 (MQ=255)
                                                            aCGCTGAACGAAGCCTACACAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCTAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTc                                                                                           <  2:213584/150‑1 (MQ=255)
                                                            aCGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCTAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTc                                                                                           <  2:4439793/150‑1 (MQ=255)
                                                             cGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCg                                                                                          <  2:1973112/150‑1 (MQ=255)
                                                                   aCGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGa                                                                                    <  2:4653909/150‑1 (MQ=255)
                                                                          cTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAAcc                                                                             <  2:2113363/150‑1 (MQ=255)
                                                                              aaaGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCAcc                                                                         >  2:4219706/1‑150 (MQ=255)
                                                                                aGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGc                                                                       >  1:3953593/1‑150 (MQ=255)
                                                                                   aTTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCAc                                                                    >  1:3441232/1‑150 (MQ=255)
                                                                                     tGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTc                                                                  <  1:1545320/150‑1 (MQ=255)
                                                                                       cAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCg                                                                >  2:3178848/1‑150 (MQ=255)
                                                                                       cAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCg                                                                >  1:2332701/1‑150 (MQ=255)
                                                                                         aTTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCt                                                              <  1:3397603/150‑1 (MQ=255)
                                                                                         aTTACAAGGCAATCGCCTCTCGGCATATGAAGAGTAGAACCTGGCCACTCTCCGAGGAGCGAAATCTCTGGCCCTTGACGATATGATTGGCAACATATTACCTGGCAAACAGGCTGATTTAGTGGTGATGCAACCCACCGCCACTCCGCt                                                              >  1:2992745/1‑150 (MQ=255)
                                                                                          ttACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCTAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTa                                                             <  1:4481313/150‑1 (MQ=255)
                                                                                          ttACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCTAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTa                                                             <  1:3275475/150‑1 (MQ=255)
                                                                                            aCAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCTAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTAca                                                           <  2:3450998/150‑1 (MQ=255)
                                                                                            aCAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTAca                                                           >  2:3334841/1‑150 (MQ=255)
                                                                                             cAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCTAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTAcag                                                          <  2:1463492/150‑1 (MQ=255)
                                                                                               aGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTAcagca                                                        >  2:1861241/1‑150 (MQ=255)
                                                                                                ggCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCTAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTAcagcag                                                       >  1:3671931/1‑150 (MQ=255)
                                                                                                 gCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACagcagc                                                      <  1:2629156/150‑1 (MQ=255)
                                                                                                   tATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTg                                                    >  2:3432625/1‑150 (MQ=255)
                                                                                                      cGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTgcgc                                                 >  1:4655873/1‑150 (MQ=255)
                                                                                                       gCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCt                                                <  1:2886021/150‑1 (MQ=255)
                                                                                                         ctctCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTAt                                              >  2:2501242/1‑150 (MQ=255)
                                                                                                         ctctCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTAt                                              >  1:3937539/1‑150 (MQ=255)
                                                                                                         ctctCGGCAGATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCTAAATCTCGGGGCCTTGACGATTTGATTGGCAACTTTGTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTAt                                              <  2:4096102/150‑1 (MQ=255)
                                                                                                          tctcGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCTAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATg                                             <  1:697936/150‑1 (MQ=255)
                                                                                                          tctcGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCTAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATg                                             <  2:833264/150‑1 (MQ=255)
                                                                                                              ggCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGTGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACaa                                         >  2:2728488/1‑150 (MQ=255)
                                                                                                                 atatGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAActc                                      >  2:3533793/1‑150 (MQ=255)
                                                                                                                   atGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTg                                    <  2:2385532/150‑1 (MQ=255)
                                                                                                                     gAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGtt                                  <  1:79816/150‑1 (MQ=255)
                                                                                                                      aaGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGttt                                 >  2:4301653/1‑150 (MQ=255)
                                                                                                                      aaGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGttt                                 <  1:3336226/150‑1 (MQ=255)
                                                                                                                       aGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTc                                >  1:4119765/1‑150 (MQ=255)
                                                                                                                            tttACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAg                           <  1:2651749/150‑1 (MQ=255)
                                                                                                                             ttACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGt                          <  1:4347828/150‑1 (MQ=255)
                                                                                                                               aCCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCg                        <  1:3247949/150‑1 (MQ=255)
                                                                                                                                  tGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACa                     <  1:3858424/150‑1 (MQ=255)
                                                                                                                                  tGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACa                     <  1:1046927/150‑1 (MQ=255)
                                                                                                                                          cTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTc             <  1:357820/150‑1 (MQ=255)
                                                                                                                                             ggcggAGCTAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtg          >  1:3984957/1‑150 (MQ=255)
                                                                                                                                             ggcggAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtg          <  1:1633553/150‑1 (MQ=255)
                                                                                                                                               cggAGCTAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgat        >  2:2560794/1‑150 (MQ=255)
                                                                                                                                               cggAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgat        <  1:3614636/150‑1 (MQ=255)
                                                                                                                                                 gAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatga      >  2:4497517/1‑150 (MQ=255)
                                                                                                                                                  aGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatgat     >  1:2956179/1‑150 (MQ=255)
                                                                                                                                                   gCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatgatg    <  2:3693807/150‑1 (MQ=255)
                                                                                                                                                   gCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatgatg    <  2:1231453/150‑1 (MQ=255)
                                                                                                                                                    cGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtgatgatga   >  2:4791274/1‑150 (MQ=255)
                                                                                                                                                     gAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGAc  >  1:512910/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAAGTTAAAGTGGGCATGGGAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGAC  >  CP009273/3020079‑3020377

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: