Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 114,825 G→A 14.6% S103F (TCC→TTC)  nadC ← quinolinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273114,8250GA14.6% 728.7 / 91.8 253S103F (TCC→TTC) nadCquinolinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (113/103);  new base A (24/13);  total (137/116)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 266 total aligned reads displayed.

GCGTATCCAACAACTGCGTGTTGGTGCCTTCCAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAGACGCCATTCTCGCGGGTGATCACCGTGGCATGAGAGCGAGA  >  CP009273/114676‑114973
                                                                                                                                                     |                                                                                                                                                    
gcgTATCCAACAACTGCGTGTTGGTGCCTTCCAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAgg                                                                                                                                                      >  1:806805/1‑150 (MQ=255)
  gTATCCAACAACTGCGTGTTGGTGCCTTCCAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGgat                                                                                                                                                    <  1:3500279/150‑1 (MQ=255)
  gTATCCAACAACTGCGTGTTGGTGCCTTCCAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGaat                                                                                                                                                    <  2:1749631/150‑1 (MQ=255)
  gTATCCAACAACTGCGTGTTGGTGCCTTCCAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGaat                                                                                                                                                    <  2:4080351/150‑1 (MQ=255)
    aTCCAACAACTGCGTGTTGGTGCCTTCCAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACACAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGaattg                                                                                                                                                  >  2:3886635/1‑150 (MQ=255)
    aTCCAACAACTGCGTGTTGGTGCCTTCCAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGgattg                                                                                                                                                  <  2:450912/150‑1 (MQ=255)
       caacaaCTGCGTGTTGGTGCCTTCCAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGgattgatt                                                                                                                                               >  1:3978403/1‑150 (MQ=255)
        aacaacTGCGTGTTGGTGCCTTCCAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGattgattg                                                                                                                                              >  1:225437/1‑150 (MQ=255)
             cTGCGTGTTGGTGCCTTCCAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCAtt                                                                                                                                         <  2:2935203/150‑1 (MQ=255)
             cTGCGTGTTGGTGCCTTCCAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCAtt                                                                                                                                         <  2:2935306/150‑1 (MQ=255)
              tGCGTGTTGGTGCCTTCCAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATtg                                                                                                                                        >  1:480155/1‑150 (MQ=255)
               gCGTGTTGGTGCCTTCCAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATtga                                                                                                                                       <  1:445235/150‑1 (MQ=255)
               gCGTGTTGGTGCCTTCCAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATtga                                                                                                                                       <  1:4255309/150‑1 (MQ=255)
                cGTGTTGGTGCCTTCCAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATtgat                                                                                                                                      <  2:3211319/150‑1 (MQ=255)
                 gtgtTGGTGCCTTCCAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATtgatg                                                                                                                                     >  2:828163/1‑150 (MQ=255)
                 gtgtTGGTGCCTTCCAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATtgatg                                                                                                                                     <  2:2287049/150‑1 (MQ=255)
                 gtgtTGGTGCCTTCCAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATtgatg                                                                                                                                     <  2:4516779/150‑1 (MQ=255)
                 gtgtTGGTGCCTTCCAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGAATTGATTGGCATtgatg                                                                                                                                     >  1:1688408/1‑150 (MQ=255)
                  tgtTGGTGCCTTCCAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGAATTGATTGGCATtgatga                                                                                                                                    >  1:4128527/1‑150 (MQ=255)
                    ttGGTGCCTTCCAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACa                                                                                                                                  >  1:743679/1‑150 (MQ=255)
                     tGGTGCCTTCCAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACAt                                                                                                                                 <  1:1477139/150‑1 (MQ=255)
                      ggTGCCTTCCAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATc                                                                                                                                <  1:3957005/150‑1 (MQ=255)
                        tGCCTTCCAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGc                                                                                                                              >  1:1746749/1‑150 (MQ=255)
                           cTTCCAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGt                                                                                                                           >  2:420368/1‑150 (MQ=255)
                               cAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTAATc                                                                                                                       <  2:763908/150‑1 (MQ=255)
                                aGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATcc                                                                                                                      >  1:3537352/1‑150 (MQ=255)
                                aGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGAATTGATTGGCATTGATGACATCGCCGTCATcc                                                                                                                      <  2:2460419/150‑1 (MQ=255)
                                  cAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGAATTGATTGGCATTGATGACATCGCCGTCATCcac                                                                                                                    <  1:3568262/150‑1 (MQ=255)
                                   aaTTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCcaca                                                                                                                   >  1:3178329/1‑150 (MQ=255)
                                    aTTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACAt                                                                                                                  <  2:600678/150‑1 (MQ=255)
                                     ttCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATg                                                                                                                 <  2:743679/150‑1 (MQ=255)
                                     ttCGACATAGTGGCGCACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATg                                                                                                                 >  2:3780753/1‑150 (MQ=255)
                                        gACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGAATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCa                                                                                                              >  1:2496012/1‑150 (MQ=255)
                                         aCATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAg                                                                                                             <  1:4529051/150‑1 (MQ=255)
                                         aCAGAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAg                                                                                                             <  2:3632568/150‑1 (MQ=255)
                                          cATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGa                                                                                                            >  2:2126803/1‑150 (MQ=255)
                                              gTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTAt                                                                                                        <  2:3272800/150‑1 (MQ=255)
                                                ggCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATgg                                                                                                      <  2:4571063/150‑1 (MQ=255)
                                                ggCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGAATTGATTGGCATTGATGACATCGCCGTCATCCACCTGCCAGATTATgg                                                                                                      <  2:568989/150‑1 (MQ=255)
                                                ggCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGAATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATgg                                                                                                      <  2:3755123/150‑1 (MQ=255)
                                                  cGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTg                                                                                                    <  2:2550578/150‑1 (MQ=255)
                                                   gTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGa                                                                                                   <  2:2630722/150‑1 (MQ=255)
                                                   gTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGAATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGa                                                                                                   >  1:280753/1‑150 (MQ=255)
                                                     aCCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACa                                                                                                 <  1:2011100/150‑1 (MQ=255)
                                                     aCCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGAGGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACa                                                                                                 >  2:2023232/1‑150 (MQ=255)
                                                          aCTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACAtcgtc                                                                                            <  1:2751461/150‑1 (MQ=255)
                                                          aCTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGAATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACAtcgtc                                                                                            <  1:2902999/150‑1 (MQ=255)
                                                             ggCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGcc                                                                                         <  1:3159831/150‑1 (MQ=255)
                                                                aaCTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCCAGTTCGAACAAGGAGTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGc                                                                                      >  2:1124378/1‑150 (MQ=255)
                                                                aaCTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGc                                                                                      >  2:3318364/1‑150 (MQ=255)
                                                                     cTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGtt                                                                                 <  1:4311783/150‑1 (MQ=255)
                                                                      tGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTg                                                                                >  1:3772926/1‑150 (MQ=255)
                                                                       gAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGa                                                                               >  1:4395320/1‑150 (MQ=255)
                                                                          aGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATa                                                                            >  1:599826/1‑150 (MQ=255)
                                                                          aGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGACTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATa                                                                            >  1:2918438/1‑150 (MQ=255)
                                                                            ggTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATaaa                                                                          >  2:2888245/1‑150 (MQ=255)
                                                                              tttGCACAAAATTAAGCTCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAaca                                                                        <  1:3323425/150‑1 (MQ=255)
                                                                               ttGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAacac                                                                       >  1:4074338/1‑150 (MQ=255)
                                                                                tGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACAcc                                                                      <  1:2439988/150‑1 (MQ=255)
                                                                                tGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACAcc                                                                      >  1:2614929/1‑150 (MQ=255)
                                                                                tGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACAcc                                                                      <  2:3976461/150‑1 (MQ=255)
                                                                                  cacaAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACctc                                                                    <  2:3554195/150‑1 (MQ=255)
                                                                                    caAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCtt                                                                  <  2:4027500/150‑1 (MQ=255)
                                                                                     aaaaTTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTc                                                                 >  1:1214300/1‑150 (MQ=255)
                                                                                       aaTTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCaa                                                               >  1:2631955/1‑150 (MQ=255)
                                                                                          tAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAAccc                                                            <  1:2687270/150‑1 (MQ=255)
                                                                                            aGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAg                                                          >  1:1548162/1‑150 (MQ=255)
                                                                                              cgcAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCg                                                        >  2:3403448/1‑150 (MQ=255)
                                                                                               gcAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGt                                                       >  1:2302619/1‑150 (MQ=255)
                                                                                                cAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGtt                                                      >  1:3448259/1‑150 (MQ=255)
                                                                                                 aGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGttt                                                     <  1:1015653/150‑1 (MQ=255)
                                                                                                      gTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGc                                                >  2:1330608/1‑150 (MQ=255)
                                                                                                          gCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCaaaa                                            <  1:2987497/150‑1 (MQ=255)
                                                                                                           cccGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAg                                           <  2:1587952/150‑1 (MQ=255)
                                                                                                           cccGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGAATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAg                                           >  2:4049867/1‑150 (MQ=255)
                                                                                                            ccGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAGa                                          >  2:3452973/1‑150 (MQ=255)
                                                                                                             cGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAGAc                                         <  2:487179/150‑1 (MQ=255)
                                                                                                             cGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAGAc                                         <  2:393857/150‑1 (MQ=255)
                                                                                                                tAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAGACGcc                                      >  2:1086011/1‑150 (MQ=255)
                                                                                                                 aaCAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAGACGCCa                                     >  1:966348/1‑150 (MQ=255)
                                                                                                                  aCAGCACGCGGGAGGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAGACGCCAt                                    >  2:296622/1‑150 (MQ=255)
                                                                                                                    aGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAGACGCCATtc                                  <  2:4020661/150‑1 (MQ=255)
                                                                                                                     gCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAGACGCCATtct                                 <  2:208567/150‑1 (MQ=255)
                                                                                                                        cgcgGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAGACGCCATTCTcgc                              >  1:1851781/1‑150 (MQ=255)
                                                                                                                           gggATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAGACGCCATTCTCGCggg                           >  2:3591835/1‑150 (MQ=255)
                                                                                                                            ggATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAGACGCCATTCTCGCGGGt                          >  2:3162036/1‑150 (MQ=255)
                                                                                                                               tGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAGACGCCATTCTCGCGGGTGAt                       >  2:119739/1‑150 (MQ=255)
                                                                                                                                tcgCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAGACGCCATTCTCGCGGGTGATc                      >  2:2407170/3‑150 (MQ=255)
                                                                                                                                gggCCTTCAAGTTCGAACAAGAATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAGACGCCATTCTCGCGGGTGATc                      <  2:1688408/150‑1 (MQ=255)
                                                                                                                                 ggCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAGACGCCATTCTCGCGGGTGATCa                     >  2:4288079/1‑150 (MQ=255)
                                                                                                                                    cTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAGACGCCATTCTCGCGGGTGATCACCg                  >  1:3473/1‑150 (MQ=255)
                                                                                                                                      tCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAGACGCCATTCTCGCGGGTGATCACCGTg                <  1:420368/150‑1 (MQ=255)
                                                                                                                                        aaGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAGACGCCATTCTCGCGGGTGATCACCGTGGc              >  1:1760343/1‑150 (MQ=255)
                                                                                                                                         aGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAGACGCCATTCTCGCGGGTGATCACCGTGGCa             >  1:2062998/1‑150 (MQ=255)
                                                                                                                                           ttCGAACAAGAATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAGACGCCATTCTCGCGGGTGATCACCGTGGCATg           >  2:2797659/1‑150 (MQ=255)
                                                                                                                                            tCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAGACGCCATTCTCGCGGGTGATCACCGTGGCATga          >  1:1984728/1‑150 (MQ=255)
                                                                                                                                              gAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACACAGCGTTTGCCGCAAAAGACGCCATTCTCGCGGGTGATCACCGTGGCATgaga        >  1:4369766/1‑150 (MQ=255)
                                                                                                                                                aCAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAGACGCCATTCTCGCGGGTGATCACCGTGGCATGAGAGc      <  2:4181860/150‑1 (MQ=255)
                                                                                                                                                aCAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAGACGCCATTCTCGCGGGTGATCACCGTGGCATGAGAGc      <  2:2847409/150‑1 (MQ=255)
                                                                                                                                                 cAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAGACGCCATTCTCGCGGGTGATCACCGTGGCATGAGAGCg     <  1:4573792/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCGTATCCAACAACTGCGTGTTGGTGCCTTCCAGCAATTCGACATAGTGGCGTACCTTACTGGCAACTCCTGAAAGGGTTTGCACAAAATTAAGCGCAGTGCGTTCGCCCGTTAACAGCACGCGGGATGGGCCTTCAAGTTCGAACAAGGATTGATTGGCATTGATGACATCGCCGTCATCCACATGCCAGATTATGGTGACATCGTCGCCTGCCAGTTGAATAAACACCTCTTCAACCCAGCGTTTGCCGCAAAAGACGCCATTCTCGCGGGTGATCACCGTGGCATGAGAGCGAGA  >  CP009273/114676‑114973

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: ‑