Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,838,672 T→A 9.4% L450Q (CTG→CAG)  ade → cryptic adenine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,838,6720TA 9.4% 1006.5 / 69.1 304L450Q (CTG→CAG) adecryptic adenine deaminase
Reads supporting (aligned to +/- strand):  ref base T (148/127);  new base A (16/13);  total (164/140)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 322 total aligned reads displayed.

GCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGCCGAAGAGATGGCGCTGGCGGTCAATCAGGTGATTCAGGATGGCG  >  CP009273/3838523‑3838821
                                                                                                                                                     |                                                                                                                                                     
   aTTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACggctggc                                                                                                                                                    <  1:1083636/150‑1 (MQ=255)
        cGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCgg                                                                                                                                               >  1:1317394/1‑150 (MQ=255)
         gCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACCTGAGCGTTACGGGCAACGGCTGGCTCCGGc                                                                                                                                              >  2:2665262/1‑150 (MQ=255)
           ccGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCtt                                                                                                                                            <  2:3804701/150‑1 (MQ=255)
            cGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTtg                                                                                                                                           >  2:1865135/1‑150 (MQ=255)
             ggAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTtgt                                                                                                                                          >  1:1967367/1‑150 (MQ=255)
              gAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTtgtg                                                                                                                                         <  2:4870865/150‑1 (MQ=255)
                aaaCGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGt                                                                                                                                       >  1:3131516/1‑150 (MQ=255)
                 aaCGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGtt                                                                                                                                      <  1:2277991/150‑1 (MQ=255)
                  aCGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGttt                                                                                                                                     >  1:350742/1‑150 (MQ=255)
                    gCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGc                                                                                                                                   >  1:5112443/1‑150 (MQ=255)
                     cTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCt                                                                                                                                  >  2:4871704/1‑150 (MQ=255)
                     cTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCt                                                                                                                                  >  1:225700/1‑150 (MQ=255)
                      tATCGGGTCATTGACGTCATCCATAACGAAGTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGCTCGCGATGATGTGAGCTTTAGTGCCGTACTTGAGCGTTACGGGCAACGGCAGGCTCCGGCTTGTGGTTTGCtt                                                                                                                                 <  2:3986364/150‑1 (MQ=255)
                       aTCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTg                                                                                                                                >  2:2339788/1‑150 (MQ=255)
                       aTCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTg                                                                                                                                <  1:1150410/150‑1 (MQ=255)
                         cGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCAGGCTCCGGCTTGTGGTTTGCTTggc                                                                                                                              >  2:1505086/1‑150 (MQ=255)
                          gggTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTggcg                                                                                                                             >  1:891509/1‑150 (MQ=255)
                           ggTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTggcgg                                                                                                                            >  2:4950002/1‑150 (MQ=255)
                            gTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTggcggc                                                                                                                           >  1:1782756/1‑150 (MQ=255)
                               aTTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCAGGCTCCGGCGTGTGGTTTGCTTGGCGGCttt                                                                                                                        >  1:3325426/1‑150 (MQ=255)
                               aTTGACGTCATCCAGAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCttt                                                                                                                        <  2:1182730/150‑1 (MQ=255)
                                  gACGTCATCCATAACGAATTGAGTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGCGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGa                                                                                                                     <  1:5148175/150‑1 (MQ=255)
                                    cGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACt                                                                                                                   <  2:2052096/150‑1 (MQ=255)
                                       cATCCATAACGAATTGATTACGCACTCACACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGACGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGTCGGCTTTGGACtgaa                                                                                                                >  1:249494/1‑150 (MQ=255)
                                         tCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACtgaatg                                                                                                              >  1:1568821/1‑150 (MQ=255)
                                         tCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACtgaatg                                                                                                              >  1:4129742/1‑150 (MQ=255)
                                           cATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACtgaatgaa                                                                                                            >  1:3010369/1‑150 (MQ=255)
                                              aaCGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGt                                                                                                         >  2:3405649/1‑150 (MQ=255)
                                              aaCGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATAGCGATGATGTGAGCTTTATTGCCTTACTTTAGCTTTACGTGCAACGGCTGGCTCCGTCTTGTGGTTTGCTTGTCGGCTTTTTACTGAATGAATGt                                                                                                         <  2:1383718/150‑1 (MQ=255)
                                                cGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTgc                                                                                                       <  2:849508/150‑1 (MQ=255)
                                                cGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTgc                                                                                                       <  1:1694565/150‑1 (MQ=255)
                                                cGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTgc                                                                                                       >  1:4322283/1‑150 (MQ=255)
                                                cGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTgc                                                                                                       <  2:3699640/150‑1 (MQ=255)
                                                     tGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCgctgg                                                                                                  >  2:2595961/1‑150 (MQ=255)
                                                      gATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCgctggc                                                                                                 <  1:1074747/150‑1 (MQ=255)
                                                      gATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCgctggc                                                                                                 <  1:2159535/150‑1 (MQ=255)
                                                       aTTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCAGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCgctggct                                                                                                <  2:621812/150‑1 (MQ=255)
                                                       aTTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCAGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCgctggct                                                                                                <  2:4267927/150‑1 (MQ=255)
                                                         tACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGc                                                                                              >  2:3125911/1‑150 (MQ=255)
                                                            gCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGAc                                                                                           >  1:4999809/1‑150 (MQ=255)
                                                            gCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGAc                                                                                           >  1:3696991/1‑150 (MQ=255)
                                                             cactcccactccAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACg                                                                                          <  2:108851/150‑1 (MQ=255)
                                                              actcccactccAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACgg                                                                                         <  1:2638045/150‑1 (MQ=255)
                                                              actcccactccAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCAGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACgg                                                                                         >  2:4629624/1‑150 (MQ=255)
                                                                 cccactccAGCGTCTACATCGTAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGAGTTACGGGCAACGGCTGGCTCCGGATTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCa                                                                                      <  1:344364/150‑1 (MQ=255)
                                                                 cccactccAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCa                                                                                      <  2:279507/150‑1 (MQ=255)
                                                                  ccactccAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAg                                                                                     <  1:2401357/150‑1 (MQ=255)
                                                                  ccactccAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTGGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTGGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAg                                                                                     <  2:3898949/150‑1 (MQ=255)
                                                                     ctccAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCa                                                                                  <  2:4950254/150‑1 (MQ=255)
                                                                       cccgCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATg                                                                                <  2:4725442/147‑1 (MQ=255)
                                                                       ccAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATg                                                                                <  1:4694182/150‑1 (MQ=255)
                                                                       ccAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATg                                                                                >  2:1161372/1‑150 (MQ=255)
                                                                         aGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGAc                                                                              >  2:3349428/1‑150 (MQ=255)
                                                                                aCAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCata                                                                       <  2:891509/150‑1 (MQ=255)
                                                                                 cAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCataa                                                                      >  1:1528816/1‑150 (MQ=255)
                                                                                 cAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCataa                                                                      >  1:4089474/1‑150 (MQ=255)
                                                                                 cAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCataa                                                                      <  2:2905653/150‑1 (MQ=255)
                                                                                  aGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCataat                                                                     >  1:595082/1‑150 (MQ=255)
                                                                                    cGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAatat                                                                   >  2:2469924/1‑150 (MQ=255)
                                                                                       aaaTGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATtgt                                                                >  1:738492/1‑150 (MQ=255)
                                                                                           ggTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTgtggtg                                                            >  2:3026144/1‑150 (MQ=255)
                                                                                            gTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGa                                                           >  1:3159046/1‑150 (MQ=255)
                                                                                            gTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGa                                                           <  1:3553988/150‑1 (MQ=255)
                                                                                             ttttGATCGCGATGATGTGAGCTTTATTGCCTTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGAt                                                          >  1:2681682/1‑150 (MQ=255)
                                                                                              tttGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATc                                                         <  1:1188289/150‑1 (MQ=255)
                                                                                              tttGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATc                                                         <  2:579064/150‑1 (MQ=255)
                                                                                                 gATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGt                                                      >  2:1029355/1‑150 (MQ=255)
                                                                                                 gATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGGTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGt                                                      <  2:3638315/150‑1 (MQ=255)
                                                                                                   tCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCg                                                    >  1:4493424/1‑150 (MQ=255)
                                                                                                    cgcgATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGc                                                   >  2:2226193/1‑150 (MQ=255)
                                                                                                        atgatgTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTg                                               <  2:2412200/150‑1 (MQ=255)
                                                                                                        atgatgTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTg                                               <  1:161534/150‑1 (MQ=255)
                                                                                                         tgatgTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGc                                              >  1:1932239/1‑150 (MQ=255)
                                                                                                          gatgTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGcc                                             <  2:738492/150‑1 (MQ=255)
                                                                                                          gatgTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGcc                                             <  1:2543421/150‑1 (MQ=255)
                                                                                                             gtgAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGATGGTGCGCTGGCTGCGACGGTCAGCCATTACAGCCCTAGTAGTGTGGTGATCGGTCGCTGTGCCGaa                                          >  1:4162896/1‑150 (MQ=255)
                                                                                                               gAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGCCGAaga                                        <  2:4085581/150‑1 (MQ=255)
                                                                                                                aGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGCCGAagag                                       >  1:3266159/1‑150 (MQ=255)
                                                                                                                aGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGCCGAagag                                       >  2:1595506/1‑150 (MQ=255)
                                                                                                                     tATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGCCGAAGAGATGgc                                  <  1:4567393/150‑1 (MQ=255)
                                                                                                                     tATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGCCGAAGAGATGgc                                  <  2:413589/150‑1 (MQ=255)
                                                                                                                       ttGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGCCGAAGAGATGgcgc                                >  2:1678796/1‑150 (MQ=255)
                                                                                                                        tGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGCCGAAGAGATGGCGCt                               <  1:3618128/150‑1 (MQ=255)
                                                                                                                         gCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGCCGAAGAGATGGCGCTg                              >  2:520555/1‑150 (MQ=255)
                                                                                                                           cGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGCCGAAGAGATGGCGCTGGc                            >  1:4837735/1‑150 (MQ=255)
                                                                                                                           cGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGCCGAAGAGATGGCGCTGGc                            >  2:3228773/1‑150 (MQ=255)
                                                                                                                            gTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGCCGAAGAGATGGCGCTGGCg                           >  1:98763/1‑150 (MQ=255)
                                                                                                                            gTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGCCGAAGAGATGGCGCTGGCg                           <  1:4871704/150‑1 (MQ=255)
                                                                                                                               cTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGCCGAAGAGATGGCGCTGGCGGTc                        <  2:4431457/150‑1 (MQ=255)
                                                                                                                                     cGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGCCGAAGAGATGGCGCTGGCGGTCAATCAg                  >  1:3679952/1‑150 (MQ=255)
                                                                                                                                      gTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGCCGAAGAGATGGCGCTGGCGGTCAATCAgg                 >  2:2574832/1‑150 (MQ=255)
                                                                                                                                        tACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGCCGAAGAGATGGCGCTGGCGGTCAATCAGGTg               >  2:1633105/1‑150 (MQ=255)
                                                                                                                                        tACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGCCGAAGAGATGGCGCTGGCGGTCAATCAGGTg               >  1:681937/1‑150 (MQ=255)
                                                                                                                                               aaCGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGCCGAAGAGATGGCGCTGGCGGTCAATCAGGTGATTCAgg        >  1:3904914/1‑150 (MQ=255)
                                                                                                                                                  ggctggctCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGCCGAAGAGATGGCGCTGGCGGTCAATCAGGTGATTCAGGATg     <  1:2899036/150‑1 (MQ=255)
                                                                                                                                                  ggctggctCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGCCGAAGAGATGGCGCTGGCGGTCAATCAGGTGATTCAGGATg     >  1:199832/1‑150 (MQ=255)
                                                                                                                                                    ctggctCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGCCGAAGAGATGGCGCTGGCGGTCAATCAGGTGATTCAGGATggc   >  1:4769143/1‑150 (MQ=255)
                                                                                                                                                     tggctCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGCCGAAGAGATGGCGCTGGCGGTCAATCAGGTGATTCAGGATggcg  >  1:3512612/1‑150 (MQ=255)
                                                                                                                                                     tggctCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGCCGAAGAGATGGCGCTGGCGGTCAATCAGGTGATTCAGGATggcg  >  2:5065682/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAATTTACGCCCGGAAAACGCTATCGGGTCATTGACGTCATCCATAACGAATTGATTACGCACTCCCACTCCAGCGTCTACAGCGAAAATGGTTTTGATCGCGATGATGTGAGCTTTATTGCCGTACTTGAGCGTTACGGGCAACGGCTGGCTCCGGCTTGTGGTTTGCTTGGCGGCTTTGGACTGAATGAAGGTGCGCTGGCTGCGACGGTCAGCCATGACAGCCATAATATTGTGGTGATCGGTCGCAGTGCCGAAGAGATGGCGCTGGCGGTCAATCAGGTGATTCAGGATGGCG  >  CP009273/3838523‑3838821

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: