Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,618,167 G→A 6.2% E346K (GAG→AAG)  nadR → trifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,618,1670GA 6.2% 1083.0 / 39.7 305E346K (GAG→AAG) nadRtrifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (145/141);  new base A (13/6);  total (158/147)
Fisher's exact test for biased strand distribution p-value = 1.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 311 total aligned reads displayed.

TTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTC  >  CP009273/4618020‑4618316
                                                                                                                                                   |                                                                                                                                                     
tttGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGgag                                                                                                                                                     <  1:1614747/150‑1 (MQ=255)
    cGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGaaca                                                                                                                                                 >  1:1148799/1‑150 (MQ=255)
    cGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGaaca                                                                                                                                                 >  2:377041/1‑150 (MQ=255)
        gAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAacac                                                                                                                                             <  2:1543968/150‑1 (MQ=255)
         aaaTATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACg                                                                                                                                            <  2:2119011/150‑1 (MQ=255)
         aaaTATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACg                                                                                                                                            <  2:2672179/150‑1 (MQ=255)
          aaTATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGc                                                                                                                                           >  2:2868048/1‑150 (MQ=255)
             atGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCgt                                                                                                                                        >  2:3497442/1‑150 (MQ=255)
             atGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCgt                                                                                                                                        >  1:4368239/1‑150 (MQ=255)
              tGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCgtg                                                                                                                                       <  2:4757291/150‑1 (MQ=255)
                ccAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACTCGCCgtggg                                                                                                                                     >  2:1022077/1‑150 (MQ=255)
                 cAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTAGATCTGGTGATCCTGCTGGAGAACAACACGCCgtgggt                                                                                                                                    >  1:3998959/1‑150 (MQ=255)
                  aaTAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCgtgggtg                                                                                                                                   >  1:689231/1‑150 (MQ=255)
                  aaTAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCgtgggtg                                                                                                                                   >  1:5070098/1‑150 (MQ=255)
                  aaTAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCgtgggtg                                                                                                                                   <  1:4273543/150‑1 (MQ=255)
                   aTAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCgtgggtgg                                                                                                                                  >  2:1215023/1‑150 (MQ=255)
                   aTAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCgtgggtgg                                                                                                                                  >  2:4828419/1‑150 (MQ=255)
                      aaGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCgg                                                                                                                               >  1:1313222/1‑150 (MQ=255)
                        gTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGAt                                                                                                                             <  1:2578744/150‑1 (MQ=255)
                        gTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGAt                                                                                                                             >  1:2874675/1‑150 (MQ=255)
                        gTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGAt                                                                                                                             >  1:406787/1‑150 (MQ=255)
                        gTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGAt                                                                                                                             >  2:4631291/1‑150 (MQ=255)
                         tGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATg                                                                                                                            <  2:2483222/150‑1 (MQ=255)
                            cATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGtt                                                                                                                         >  1:3189645/1‑150 (MQ=255)
                             ttttATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGttt                                                                                                                        <  2:3159775/149‑1 (MQ=255)
                               ttATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTAc                                                                                                                      <  1:4389813/150‑1 (MQ=255)
                                tATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGAGGGTTTACg                                                                                                                     >  2:1191499/1‑150 (MQ=255)
                                 aTCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGc                                                                                                                    >  1:795288/1‑150 (MQ=255)
                                 aTCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGc                                                                                                                    >  1:3409064/1‑150 (MQ=255)
                                 aTCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGc                                                                                                                    >  2:3465365/1‑150 (MQ=255)
                                      tACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCt                                                                                                               <  1:1501056/150‑1 (MQ=255)
                                       aCCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTc                                                                                                              <  2:4233147/150‑1 (MQ=255)
                                         cGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCgg                                                                                                            <  2:764089/150‑1 (MQ=255)
                                            ttttGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAg                                                                                                         <  2:232467/150‑1 (MQ=255)
                                              ttGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGtt                                                                                                       <  2:2991316/150‑1 (MQ=255)
                                              ttGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGtt                                                                                                       <  2:353469/150‑1 (MQ=255)
                                              ttGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGtt                                                                                                       <  2:1593232/150‑1 (MQ=255)
                                               tGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTc                                                                                                      <  1:1837964/150‑1 (MQ=255)
                                                gTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCg                                                                                                     >  2:1970056/1‑150 (MQ=255)
                                                gTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCg                                                                                                     >  1:2547501/1‑150 (MQ=255)
                                                 tCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCgg                                                                                                    <  1:3015447/150‑1 (MQ=255)
                                                   accacTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTg                                                                                                  >  2:568602/1‑150 (MQ=255)
                                                      acTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGAAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGAt                                                                                               <  2:1700195/150‑1 (MQ=255)
                                                          aGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCa                                                                                           >  2:4142456/1‑150 (MQ=255)
                                                              gTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCTGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAaga                                                                                       <  1:2378953/150‑1 (MQ=255)
                                                              gTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAaga                                                                                       <  1:3588832/150‑1 (MQ=255)
                                                               ttCTGCAAAAAGTACGCAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGACTACCGTTTCGATCTGGTGATCCTGCTGGAGAACACCACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAagag                                                                                      <  2:4560870/150‑1 (MQ=255)
                                                                 cTGCAAAAAGTACGACGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGtt                                                                                    <  1:540465/150‑1 (MQ=255)
                                                                  tGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGAAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTc                                                                                   <  1:1033894/150‑1 (MQ=255)
                                                                   gCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAAAACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTcc                                                                                  >  2:802438/1‑150 (MQ=255)
                                                                   gCAAAAAGGACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCTGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTcc                                                                                  <  2:4430978/150‑1 (MQ=255)
                                                                    cAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCa                                                                                 >  2:206199/1‑150 (MQ=255)
                                                                     aaaaaGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAg                                                                                >  2:1099685/1‑150 (MQ=255)
                                                                     aaaaaGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAg                                                                                >  1:2098771/1‑150 (MQ=255)
                                                                     aaaaaGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAg                                                                                <  2:1295328/150‑1 (MQ=255)
                                                                      aaaaGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGa                                                                               >  2:5169113/1‑150 (MQ=255)
                                                                      aaaaGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGAGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGa                                                                               >  1:1993989/1‑150 (MQ=255)
                                                                        aaGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAAc                                                                             <  2:3373307/150‑1 (MQ=255)
                                                                          gTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACtt                                                                           >  2:3676487/1‑150 (MQ=255)
                                                                           tACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTg                                                                          >  2:4687892/1‑150 (MQ=255)
                                                                               aaGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCtgg                                                                      <  1:2456518/150‑1 (MQ=255)
                                                                                 gggCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCtggtg                                                                    >  2:3211053/1‑150 (MQ=255)
                                                                                      ggAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGAt                                                               >  2:50998/1‑150 (MQ=255)
                                                                                             ccGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGaa                                                        >  1:2726494/1‑150 (MQ=255)
                                                                                              cGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAag                                                       <  1:3054692/150‑1 (MQ=255)
                                                                                                 tCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAgaga                                                    <  1:4510747/150‑1 (MQ=255)
                                                                                                    tGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGAAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACa                                                 >  1:4313048/1‑150 (MQ=255)
                                                                                                     gCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACaa                                                <  1:3938889/150‑1 (MQ=255)
                                                                                                     gCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACaa                                                <  2:1662107/150‑1 (MQ=255)
                                                                                                       aGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAata                                              <  2:3860146/150‑1 (MQ=255)
                                                                                                          cgcTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCg                                           >  1:4905855/1‑150 (MQ=255)
                                                                                                          cgcTGATTGATGAAGACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCg                                           <  2:4321882/150‑1 (MQ=255)
                                                                                                            cTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGaa                                         >  1:5133269/1‑150 (MQ=255)
                                                                                                            cTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGaa                                         >  1:2797916/1‑150 (MQ=255)
                                                                                                             tgattgatGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAAt                                        <  2:2509592/150‑1 (MQ=255)
                                                                                                             tgattgatGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAAt                                        <  2:3675124/150‑1 (MQ=255)
                                                                                                               attgatGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAAttt                                      >  2:1234966/1‑150 (MQ=255)
                                                                                                                ttgatGAATACCTCTTCGATCTGGTGATCCTGCTGGATAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGTACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTtg                                     <  1:4784881/150‑1 (MQ=255)
                                                                                                                ttgatGAATACCGTTTAGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGACCAATATCGAATTtg                                     >  1:2218080/1‑150 (MQ=255)
                                                                                                                     gaATACCGTTTCGATCTGGTTATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCgg                                >  1:3703029/1‑150 (MQ=255)
                                                                                                                       aTACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGAAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGt                              >  2:71470/1‑150 (MQ=255)
                                                                                                                        tACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGtt                             <  2:819370/150‑1 (MQ=255)
                                                                                                                        tACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGtt                             <  1:3793447/150‑1 (MQ=255)
                                                                                                                            gTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAag                         >  2:4027174/1‑150 (MQ=255)
                                                                                                                            gTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAag                         <  1:2712382/150‑1 (MQ=255)
                                                                                                                             tttCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAaga                        >  1:2413152/1‑150 (MQ=255)
                                                                                                                             tttCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAaga                        <  1:3822206/150‑1 (MQ=255)
                                                                                                                                     tGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGAtt                >  2:1354182/1‑150 (MQ=255)
                                                                                                                                        tGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACg             <  1:4655940/150‑1 (MQ=255)
                                                                                                                                         gATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGACTTTGTGCGGGTTGAAGAGGAAGATTACGa            >  1:4997022/1‑150 (MQ=255)
                                                                                                                                         gATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGAGGGTTGACGAGGACGATTACGa            >  1:3146839/1‑150 (MQ=255)
                                                                                                                                          aTCCTGCTGGAGAACAACCCGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGAc           <  2:2443866/150‑1 (MQ=255)
                                                                                                                                          aTCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGAc           <  2:5035656/150‑1 (MQ=255)
                                                                                                                                           tCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACa          <  2:1787307/150‑1 (MQ=255)
                                                                                                                                            ccTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAg         <  1:2185255/150‑1 (MQ=255)
                                                                                                                                              tgctgGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTc       >  1:4581025/1‑150 (MQ=255)
                                                                                                                                              tgctgGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTc       <  2:689231/150‑1 (MQ=255)
                                                                                                                                                ctgGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGt     >  2:1474071/1‑150 (MQ=255)
                                                                                                                                                 tgGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGtt    >  1:4449997/1‑150 (MQ=255)
                                                                                                                                                 tgGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGtt    <  2:1128278/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTC  >  CP009273/4618020‑4618316

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: