Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,825,763 C→A 26.2% A159A (GCG→GCT)  astC ← succinylornithine transaminase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,825,7630CA26.2% 473.3 / 158.7 220A159A (GCG→GCT) astCsuccinylornithine transaminase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (107/55);  new base A (30/26);  total (139/81)
Fisher's exact test for biased strand distribution p-value = 1.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 225 total aligned reads displayed.

GACCTTGTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTT  >  CP009273/1825615‑1825911
                                                                                                                                                    |                                                                                                                                                    
gACCTTGTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGgcg                                                                                                                                                     >  2:2974242/1‑150 (MQ=255)
  ccTTGTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCa                                                                                                                                                   >  1:4835596/1‑150 (MQ=255)
     tGTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCaaaa                                                                                                                                                <  1:2424573/150‑1 (MQ=255)
       tAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATc                                                                                                                                              <  1:413851/150‑1 (MQ=255)
        aaaaaCGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATcc                                                                                                                                             >  2:3023010/1‑150 (MQ=255)
        aaaaaCGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATcc                                                                                                                                             <  1:2488628/150‑1 (MQ=255)
         aaaaCGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCt                                                                                                                                            >  2:3349610/1‑150 (MQ=255)
         aaaaCGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCt                                                                                                                                            >  1:2850456/1‑150 (MQ=255)
            aCGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCtgtg                                                                                                                                         >  2:5116036/1‑150 (MQ=255)
                 ttGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAg                                                                                                                                    <  1:3275743/150‑1 (MQ=255)
                  tGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGATCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCAGCCGGCAGTGGCGAAAAATCCTGTGAATAgg                                                                                                                                   >  1:2994982/1‑150 (MQ=255)
                    ctggctggCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAggct                                                                                                                                 >  2:1532376/1‑150 (MQ=255)
                     tggctggcACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAggctg                                                                                                                                >  2:2244078/1‑150 (MQ=255)
                        ctggcACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAggctggct                                                                                                                             <  1:2551291/150‑1 (MQ=255)
                        ctggcACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAggctggct                                                                                                                             <  1:182803/150‑1 (MQ=255)
                             accacACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCa                                                                                                                        >  2:3873169/1‑150 (MQ=255)
                               cacaCCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCAcc                                                                                                                      >  1:955695/1‑150 (MQ=255)
                                  acCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGc                                                                                                                   >  2:3746867/1‑150 (MQ=255)
                                  acCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGc                                                                                                                   >  2:4405966/1‑150 (MQ=255)
                                    cGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCAc                                                                                                                 >  2:4522145/1‑150 (MQ=255)
                                     gCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACt                                                                                                                <  2:3574293/150‑1 (MQ=255)
                                     gCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACt                                                                                                                >  2:5194183/1‑150 (MQ=255)
                                       nttCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGa                                                                                                              >  1:5182160/2‑150 (MQ=255)
                                       cTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGa                                                                                                              >  1:1335085/1‑150 (MQ=255)
                                          cccccTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAg                                                                                                           <  1:2294116/150‑1 (MQ=255)
                                           ccccTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGt                                                                                                          >  2:2360855/1‑150 (MQ=255)
                                            cccTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTa                                                                                                         <  2:1584757/150‑1 (MQ=255)
                                             ccTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTaa                                                                                                        <  1:4904613/150‑1 (MQ=255)
                                              cTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTaaa                                                                                                       >  2:4342184/1‑150 (MQ=255)
                                                   tGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCg                                                                                                  <  1:3228743/150‑1 (MQ=255)
                                                   tGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCg                                                                                                  >  1:700788/1‑150 (MQ=255)
                                                   tGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCg                                                                                                  <  2:579461/150‑1 (MQ=255)
                                                    gggTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGt                                                                                                 >  1:1686413/1‑150 (MQ=255)
                                                     ggTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTg                                                                                                >  1:1879118/1‑150 (MQ=255)
                                                     ggTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTg                                                                                                <  2:1229037/150‑1 (MQ=255)
                                                      gTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGc                                                                                               >  2:3615453/1‑150 (MQ=255)
                                                       ttCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCg                                                                                              >  2:4647923/1‑150 (MQ=255)
                                                        tCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGa                                                                                             >  1:3660551/1‑150 (MQ=255)
                                                         cGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGAc                                                                                            >  2:4665286/1‑150 (MQ=255)
                                                          gACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGAcc                                                                                           >  1:1300731/1‑150 (MQ=255)
                                                            cAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCAt                                                                                         <  2:3998709/150‑1 (MQ=255)
                                                             aaTCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATg                                                                                        <  2:2328404/150‑1 (MQ=255)
                                                                 aCCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAAc                                                                                    <  2:5097044/150‑1 (MQ=255)
                                                                   cGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGc                                                                                  >  1:4268588/1‑150 (MQ=255)
                                                                   cGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGc                                                                                  >  1:3614315/1‑150 (MQ=255)
                                                                   cGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGc                                                                                  >  1:3178634/1‑150 (MQ=255)
                                                                    gCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCa                                                                                 <  2:1381190/150‑1 (MQ=255)
                                                                        aGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCAtttt                                                                             <  2:193631/150‑1 (MQ=255)
                                                                           tAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGa                                                                          >  2:4813149/1‑150 (MQ=255)
                                                                              agTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACg                                                                       >  2:371954/1‑150 (MQ=255)
                                                                              agTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAAca                                                                       >  1:3876916/1‑149 (MQ=255)
                                                                                 cgtcAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCa                                                                    <  1:2431703/150‑1 (MQ=255)
                                                                                  gtcAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCAc                                                                   >  2:2884217/1‑150 (MQ=255)
                                                                                    cAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGa                                                                 >  1:2054215/1‑150 (MQ=255)
                                                                                     aaTCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGAt                                                                >  2:2736768/1‑150 (MQ=255)
                                                                                        cAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGcc                                                             >  2:1495861/1‑150 (MQ=255)
                                                                                        cAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGcc                                                             <  2:4871209/150‑1 (MQ=255)
                                                                                        cAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGcc                                                             >  2:555073/1‑150 (MQ=255)
                                                                                        cAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGcc                                                             <  2:2564685/150‑1 (MQ=255)
                                                                                         aGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCg                                                            <  1:1013237/150‑1 (MQ=255)
                                                                                              gcTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAACATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCtt                                                       <  2:2330458/150‑1 (MQ=255)
                                                                                                tGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTAt                                                     <  1:1982098/150‑1 (MQ=255)
                                                                                                 ggCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATg                                                    <  1:3154416/150‑1 (MQ=255)
                                                                                                   cAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGc                                                  <  2:782057/150‑1 (MQ=255)
                                                                                                   cAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGc                                                  >  1:2480520/1‑150 (MQ=255)
                                                                                                       gTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGcc                                              >  1:1102938/1‑150 (MQ=255)
                                                                                                       gTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGcc                                              >  1:3213503/1‑150 (MQ=255)
                                                                                                        ttAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCg                                             >  2:688535/1‑150 (MQ=255)
                                                                                                         tAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGt                                            >  1:2167409/1‑150 (MQ=255)
                                                                                                          aaTATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTa                                           <  1:1324177/150‑1 (MQ=255)
                                                                                                           atatCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAg                                          >  1:2661234/1‑150 (MQ=255)
                                                                                                           atatCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAg                                          <  1:3436685/150‑1 (MQ=255)
                                                                                                             atCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCg                                        >  1:3104018/1‑150 (MQ=255)
                                                                                                               cGTTATATGCAGCCTGACGAATATCCGCCGGCCGTGGCGCAAAAGCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGt                                      <  2:10266/150‑1 (MQ=255)
                                                                                                                 ttATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCg                                    <  2:44581/150‑1 (MQ=255)
                                                                                                                  tataTGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGt                                   >  1:1851600/1‑150 (MQ=255)
                                                                                                                   atatGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTg                                  >  2:483745/1‑150 (MQ=255)
                                                                                                                     atGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGTCTGCCGTAGCGGTCGTGAg                                >  2:4488252/1‑150 (MQ=255)
                                                                                                                      tGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGc                               >  2:1338704/1‑150 (MQ=255)
                                                                                                                        cagcaTGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGa                             >  1:4666814/1‑150 (MQ=255)
                                                                                                                         agcaTGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGaa                            >  2:1792301/1‑150 (MQ=255)
                                                                                                                          gcaTGACGAATATCCGCCGGCAGTGGAGCATAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAAt                           <  1:4595907/150‑1 (MQ=255)
                                                                                                                           caTGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAAtt                          >  1:57448/1‑150 (MQ=255)
                                                                                                                           caTGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAAtt                          >  1:3835376/1‑150 (MQ=255)
                                                                                                                           caTGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAAtt                          <  1:4520428/150‑1 (MQ=255)
                                                                                                                           caTGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAAtt                          <  1:1933214/150‑1 (MQ=255)
                                                                                                                            aTGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAAttt                         <  1:1797970/150‑1 (MQ=255)
                                                                                                                              gACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTAc                       >  2:3512675/1‑150 (MQ=255)
                                                                                                                                cGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTAcgc                     >  2:4814029/1‑150 (MQ=255)
                                                                                                                                   atatCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGcc                  <  2:4232086/150‑1 (MQ=255)
                                                                                                                                   atatCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGcc                  >  2:1715272/1‑150 (MQ=255)
                                                                                                                                   atatCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGcc                  >  2:518928/1‑150 (MQ=255)
                                                                                                                                    tatCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCa                 >  2:252740/1‑150 (MQ=255)
                                                                                                                                      tCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTACACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGt               >  2:3135563/1‑150 (MQ=255)
                                                                                                                                        cgccgGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGttt             >  1:115903/1‑150 (MQ=255)
                                                                                                                                         gccgGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGtttt            >  2:2832395/1‑150 (MQ=255)
                                                                                                                                         gccgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGtttt            >  2:2154742/1‑150 (MQ=255)
                                                                                                                                            ggCAGTGGAGCAAAATCCTGTGAATAGGCTGGTTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAgc         >  1:233997/1‑150 (MQ=255)
                                                                                                                                                gTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGcgccg     >  1:4387354/1‑150 (MQ=255)
                                                                                                                                                 tGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGcgccgc    >  1:1081877/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GACCTTGTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTT  >  CP009273/1825615‑1825911

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