Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,839,034 C→A 49.7% Q571K (CAG→AAG)  ade → cryptic adenine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,839,0340CA49.7% 119.7 / inf 271Q571K (CAG→AAG) adecryptic adenine deaminase
Reads supporting (aligned to +/- strand):  ref base C (73/63);  new base A (75/60);  total (148/123)
Fisher's exact test for biased strand distribution p-value = 8.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 276 total aligned reads displayed.

GATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCG  >  CP009273/3838886‑3839183
                                                                                                                                                    |                                                                                                                                                     
gatgaGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCa                                                                                                                                                      <  2:8527/150‑1 (MQ=255)
    aGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAgggg                                                                                                                                                  >  1:950728/1‑150 (MQ=255)
         cGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTAtt                                                                                                                                             >  1:3342233/1‑150 (MQ=255)
          gACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTAttt                                                                                                                                            >  1:3189118/1‑150 (MQ=255)
            cacGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGa                                                                                                                                          <  1:717654/150‑1 (MQ=255)
            cacGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGa                                                                                                                                          >  1:3408406/1‑150 (MQ=255)
             acGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGAt                                                                                                                                         <  1:4229855/150‑1 (MQ=255)
             acGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGAt                                                                                                                                         <  1:1628241/150‑1 (MQ=255)
              cgtcgcAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATg                                                                                                                                        <  2:3754548/147‑1 (MQ=255)
              cGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATg                                                                                                                                        <  2:1762420/150‑1 (MQ=255)
                      tCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAg                                                                                                                                <  2:2456860/150‑1 (MQ=255)
                      tCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAg                                                                                                                                <  2:4482328/150‑1 (MQ=255)
                         cTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGttt                                                                                                                             >  1:2422225/1‑150 (MQ=255)
                         cTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGttt                                                                                                                             <  1:3058830/150‑1 (MQ=255)
                           ggCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGc                                                                                                                           <  2:4302202/150‑1 (MQ=255)
                           ggCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGc                                                                                                                           <  2:452929/150‑1 (MQ=255)
                             cGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCt                                                                                                                         <  2:946942/150‑1 (MQ=255)
                             cGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCt                                                                                                                         >  2:1000469/1‑150 (MQ=255)
                              ggAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCtt                                                                                                                        >  2:2742285/1‑150 (MQ=255)
                                  cAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACt                                                                                                                    <  1:1074647/150‑1 (MQ=255)
                                    aaTTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTAc                                                                                                                  >  1:2250775/1‑150 (MQ=255)
                                     aTTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACg                                                                                                                 >  2:2369698/1‑150 (MQ=255)
                                         aCGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTgg                                                                                                             >  2:1402093/1‑150 (MQ=255)
                                              ttGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTc                                                                                                        >  1:873627/1‑150 (MQ=255)
                                               tGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCa                                                                                                       >  2:2246452/1‑150 (MQ=255)
                                               tGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCa                                                                                                       >  1:2718556/1‑150 (MQ=255)
                                               tGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCa                                                                                                       >  1:2718475/1‑150 (MQ=255)
                                                gAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCAc                                                                                                      >  1:1213905/1‑150 (MQ=255)
                                                gAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCAc                                                                                                      >  2:2607230/1‑150 (MQ=255)
                                                             cGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATaaaaa                                                                                         <  1:2033362/150‑1 (MQ=255)
                                                              gTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGACTAATaaaaaa                                                                                        >  2:2077857/1‑150 (MQ=255)
                                                               tGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATaaaaaaa                                                                                       >  2:2902870/1‑150 (MQ=255)
                                                                gAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAg                                                                                      >  1:3665902/1‑150 (MQ=255)
                                                                gAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAg                                                                                      >  1:3665910/1‑150 (MQ=255)
                                                                 aaTGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGc                                                                                     <  2:4428936/150‑1 (MQ=255)
                                                                 aaTGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGc                                                                                     <  2:523086/150‑1 (MQ=255)
                                                                  aTGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGcc                                                                                    <  2:2363114/150‑1 (MQ=255)
                                                                   tGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGccc                                                                                   <  2:3865642/150‑1 (MQ=255)
                                                                             tACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGc                                                                         >  1:76704/1‑150 (MQ=255)
                                                                                   aTGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCa                                                                   <  1:4269101/150‑1 (MQ=255)
                                                                                    tGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCaa                                                                  <  1:126280/150‑1 (MQ=255)
                                                                                    tGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCaa                                                                  <  1:125943/150‑1 (MQ=255)
                                                                                     gAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAATAAGCCCGGCGTCATGCCGGGCaca                                                                 >  2:2604818/1‑148 (MQ=255)
                                                                                     gAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCaaa                                                                 >  2:4290564/1‑150 (MQ=255)
                                                                                       gCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAATGCCCGGCGTCATGCCGGGCAAAAg                                                               >  2:3738218/1‑150 (MQ=255)
                                                                                       gCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAg                                                               >  2:1046108/1‑150 (MQ=255)
                                                                                        ccGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGt                                                              <  1:3703949/150‑1 (MQ=255)
                                                                                        ccGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGt                                                              <  1:2011285/150‑1 (MQ=255)
                                                                                          gTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCa                                                            >  2:573469/1‑150 (MQ=255)
                                                                                          gTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCa                                                            <  1:3734322/150‑1 (MQ=255)
                                                                                            ttATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCAcc                                                          >  1:1955856/1‑150 (MQ=255)
                                                                                            ttATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCAcc                                                          >  1:588959/1‑150 (MQ=255)
                                                                                            ttATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCAcc                                                          >  1:2841864/1‑150 (MQ=255)
                                                                                             tATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCa                                                         <  2:1162696/150‑1 (MQ=255)
                                                                                             tATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCa                                                         <  2:1163209/150‑1 (MQ=255)
                                                                                             tATTCAGATGGCGTTTCTTTCACTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCa                                                         >  1:4451455/1‑150 (MQ=255)
                                                                                                tCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGtt                                                      <  1:4064314/150‑1 (MQ=255)
                                                                                                  aGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTAc                                                    >  1:1247800/1‑150 (MQ=255)
                                                                                                    aTGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTCACCAGACAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGt                                                  <  1:1869287/150‑1 (MQ=255)
                                                                                                     tGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGtt                                                 <  1:1713286/150‑1 (MQ=255)
                                                                                                      ggCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGttt                                                >  1:612356/1‑150 (MQ=255)
                                                                                                       gCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTa                                               >  1:1017165/1‑150 (MQ=255)
                                                                                                             ctttctCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCCGTTACGTTTATGCCAc                                         >  2:1801271/1‑150 (MQ=255)
                                                                                                             ctttctCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCAc                                         >  2:4053964/1‑150 (MQ=255)
                                                                                                              tttctCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACt                                        >  1:2604739/1‑150 (MQ=255)
                                                                                                              tttctCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACt                                        >  2:1736100/1‑150 (MQ=255)
                                                                                                               ttctCTGCCAGTGATCCCCGCGCTAAAACTACCCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTg                                       <  2:4434979/150‑1 (MQ=255)
                                                                                                               ttctCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTg                                       <  2:2400063/150‑1 (MQ=255)
                                                                                                                 ctctGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTc                                     >  1:3094142/1‑150 (MQ=255)
                                                                                                                   ctGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCaa                                   >  2:4145902/1‑150 (MQ=255)
                                                                                                                   ctGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCaa                                   >  2:4145656/1‑150 (MQ=255)
                                                                                                                   ctGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCaa                                   >  2:3381409/1‑150 (MQ=255)
                                                                                                                     gCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAAct                                 <  2:1954593/150‑1 (MQ=255)
                                                                                                                       cAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAActgc                               <  1:2344152/150‑1 (MQ=255)
                                                                                                                       cAGTGATCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAActgc                               >  2:3943345/1‑150 (MQ=255)
                                                                                                                         gTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAActgctg                             >  2:1458994/1‑150 (MQ=255)
                                                                                                                          tGATCCCCGCGCTAAAACTAACCCGTCAGGGGCTATTTGAGGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGa                            <  2:3565971/150‑1 (MQ=255)
                                                                                                                          tGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGa                            >  1:3225209/1‑150 (MQ=255)
                                                                                                                            aTCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAt                          <  1:4314722/150‑1 (MQ=255)
                                                                                                                            aTCCCCGCGCTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAt                          <  1:3546929/150‑1 (MQ=255)
                                                                                                                              ccccGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAttt                        >  1:1904219/1‑150 (MQ=255)
                                                                                                                               cccGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACTCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCACAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAtttt                       >  1:1342674/1‑150 (MQ=255)
                                                                                                                               cccGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAtttt                       >  2:1578723/1‑150 (MQ=255)
                                                                                                                                  gcgcTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAttttttt                    <  2:3234121/150‑1 (MQ=255)
                                                                                                                                    gcTAAAACTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTtct                  >  2:2134161/1‑150 (MQ=255)
                                                                                                                                     cTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTactc                 >  2:584645/1‑150 (MQ=255)
                                                                                                                                        aaaCTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTcgcg              >  1:327631/1‑150 (MQ=255)
                                                                                                                                        aaaCTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTcgcg              <  1:2349173/150‑1 (MQ=255)
                                                                                                                                         caCTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCgcgg             <  2:2124905/149‑1 (MQ=255)
                                                                                                                                         aaCTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCgcgg             <  2:4370363/150‑1 (MQ=255)
                                                                                                                                          aCTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCgcggc            <  1:3923027/150‑1 (MQ=255)
                                                                                                                                          aCTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCgcggc            <  1:3020012/150‑1 (MQ=255)
                                                                                                                                           cTAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTCCGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCgcggcg           <  2:3265361/150‑1 (MQ=255)
                                                                                                                                            tAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGcggcgg          >  2:3404497/1‑150 (MQ=255)
                                                                                                                                            tAACCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGcggcgg          >  2:514748/1‑150 (MQ=255)
                                                                                                                                             aaCCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGa         >  1:1418820/1‑150 (MQ=255)
                                                                                                                                             aaCCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGa         >  1:665306/1‑150 (MQ=255)
                                                                                                                                              aCCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGAt        >  1:264177/1‑150 (MQ=255)
                                                                                                                                              aCCAGTAAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGAt        >  2:263329/1‑150 (MQ=255)
                                                                                                                                                    cAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCg  >  2:2721415/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCG  >  CP009273/3838886‑3839183

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: