Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,381,462 G→A 11.6% E14E (GAG→GAA)  orn → oligoribonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,381,4620GA11.6% 1020.4 / 90.2 326E14E (GAG→GAA) ornoligoribonuclease
Reads supporting (aligned to +/- strand):  ref base G (165/123);  new base A (24/14);  total (189/137)
Fisher's exact test for biased strand distribution p-value = 6.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 327 total aligned reads displayed.

TATTTTTACTCAATCGTACAGACTCCTGGTCGCCCCTGATGGGCAAAACATCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTCTGATGAACAGCTGGCGCTGATGGATGACTGGAACG  >  CP009273/4381314‑4381607
                                                                                                                                                    |                                                                                                                                                 
tATTTTTACTCAATCGTACAGACTCCTGGTCGCCCCTGATGGGCAAAACATCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTgaga                                                                                                                                                  >  1:743833/1‑150 (MQ=255)
   ttttACTCAATCGTACAGACTCCTGGTCGCCCCTGATGGGCAAAACATCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGa                                                                                                                                               <  2:3273075/150‑1 (MQ=255)
    tttACTCAATCGTACAGACTCCTGGTCGCCCCTGATGGGCAAAACATCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAAATGAc                                                                                                                                              >  2:3016906/1‑150 (MQ=255)
     ttACTCAATCGTACAGACTCCTGGTCGCCCCTGATGGGCAAAACATCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGAcc                                                                                                                                             <  1:2776806/150‑1 (MQ=255)
        cTCAATCGTACAGACTCCTGGTCGCCCCTGATGGGCAAAACATCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGt                                                                                                                                          >  2:1707644/1‑150 (MQ=255)
        cTCAATCGTACAGACTCCTGGTCGCCCCTGATGGGCAAAACATCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGt                                                                                                                                          >  1:3812472/1‑150 (MQ=255)
            aTCGTACAGACTCCTGGTCGCCCCTGATGGGCAAAACATCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAAATGACCGGTCTgg                                                                                                                                      <  2:1592445/150‑1 (MQ=255)
               gTACAGACTCCTGGTCGCCCCTGATGGGCAAAACATCTATGATACACGCAATTTTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATc                                                                                                                                   <  1:708425/150‑1 (MQ=255)
                tACAGACTCCTGGTCGCCCCTGATGGGCAAAACATCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATcc                                                                                                                                  >  2:2898870/1‑150 (MQ=255)
                 aCAGACTCCTGGTCGCGCCTGATGGGCAAAACATCTATGATACCCGCAATTGTGGATCAATATTGTCACTGTGAATTGGTGGAAAATAGCATGAGTGCCATTGAATACAACCTGATTTGGATCGATCTTGAGATGACCGGTATGGATccc                                                                                                                                 <  1:17210/150‑1 (MQ=255)
                 aCAGACTCCTGGTCGCCCCTGATGGGCAAAACATCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATccc                                                                                                                                 >  1:2841963/1‑150 (MQ=255)
                 aCAGACTCCTGGTCGCCCCTGATGGGCAAAACATCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCGGATTTGGATCGATCTTGAGATGACCGGTCTGGATccc                                                                                                                                 >  2:1083314/1‑150 (MQ=255)
                   aGACTCCTGGTCGCCCCTGATGGGCAAAACATCTATTATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGa                                                                                                                               <  1:1281354/150‑1 (MQ=255)
                   aGACTCCTGGTCGCCCCTGATGGGCAAAACATCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGa                                                                                                                               >  1:3190546/1‑150 (MQ=255)
                    gACTCCTGGTCGCCCCTGATGGGCAAAACATCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAg                                                                                                                              >  1:3483935/1‑150 (MQ=255)
                      cTCCTGGTCGCCCCTGATGGGCAAAACATCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAgcg                                                                                                                            <  2:514911/150‑1 (MQ=255)
                               gCCCCTGATGGGCAAAACATCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCat                                                                                                                   >  2:3028755/1‑150 (MQ=255)
                                ccccTGATGGGCAAAACATCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCatt                                                                                                                  >  2:3597509/1‑150 (MQ=255)
                                 cccTGATGGGCAAAACATCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCatta                                                                                                                 >  1:628881/1‑150 (MQ=255)
                                  ccTGATGGGCAAAACATCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCattat                                                                                                                <  2:2762402/150‑1 (MQ=255)
                                    tGATGGGCAAAACATCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAAATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTg                                                                                                              >  2:3808981/1‑150 (MQ=255)
                                       tGGGCAAAACATCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAAATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTgaga                                                                                                           <  1:1086428/150‑1 (MQ=255)
                                        gggCAAAACATCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGAt                                                                                                          <  2:2654801/150‑1 (MQ=255)
                                           cAAAACATCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGc                                                                                                       <  2:589894/150‑1 (MQ=255)
                                            aaaaCATCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGcc                                                                                                      <  2:3628486/150‑1 (MQ=255)
                                              aaCATCTATGATACACGCAATTGTGGATCAGTATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCAc                                                                                                    >  1:4274795/1‑150 (MQ=255)
                                               aCATCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAAATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACg                                                                                                   <  2:3152577/150‑1 (MQ=255)
                                                 aTCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAAATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCt                                                                                                 >  2:4059789/1‑150 (MQ=255)
                                                  tCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTg                                                                                                <  2:64235/150‑1 (MQ=255)
                                                     aTGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTg                                                                                             >  1:4163135/1‑150 (MQ=255)
                                                     aTGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTg                                                                                             <  1:4368383/150‑1 (MQ=255)
                                                     aTGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTg                                                                                             <  1:1408348/150‑1 (MQ=255)
                                                     aTGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTg                                                                                             >  1:3171670/1‑150 (MQ=255)
                                                       gATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGAc                                                                                           <  1:281109/150‑1 (MQ=255)
                                                       gATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAAATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGAc                                                                                           >  2:1883875/1‑150 (MQ=255)
                                                        aTACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGAcc                                                                                          >  2:348944/1‑150 (MQ=255)
                                                        aTACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGAcc                                                                                          <  1:3263552/150‑1 (MQ=255)
                                                        aTACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGAcc                                                                                          >  1:1793517/1‑150 (MQ=255)
                                                         tACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAAATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCg                                                                                         <  2:3550878/150‑1 (MQ=255)
                                                           cacGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGAt                                                                                       <  2:280686/150‑1 (MQ=255)
                                                            acGCAATTGTGGATAAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTTCCAATGAAAACAAACTGATTTGGATCGATCTTGAAATGACAGGTCTGGATCCAGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATg                                                                                      >  2:3976747/1‑150 (MQ=255)
                                                             cGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGc                                                                                     >  1:2902515/1‑150 (MQ=255)
                                                                aaTTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCaa                                                                                  <  2:2357774/150‑1 (MQ=255)
                                                                 aTTGTGGATCAATATAGTCACTGTGAATGTGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTAGGATCGATCTTGAGATGACCGGTCTGGATCCCTAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGAAGCCAAc                                                                                 >  2:2294374/1‑150 (MQ=255)
                                                                  ttGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAAcc                                                                                >  1:2236343/1‑150 (MQ=255)
                                                                    gtgGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTg                                                                              <  2:4601240/150‑1 (MQ=255)
                                                                      ggATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGaa                                                                            <  1:4334848/150‑1 (MQ=255)
                                                                       gATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACACCCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAat                                                                           >  2:1916274/1‑150 (MQ=255)
                                                                        aTCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAata                                                                          >  2:1246783/1‑150 (MQ=255)
                                                                             tataGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTg                                                                     <  1:3793589/150‑1 (MQ=255)
                                                                             tataGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTg                                                                     >  2:3201199/1‑150 (MQ=255)
                                                                              ataGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTgg                                                                    <  1:1859204/150‑1 (MQ=255)
                                                                               taGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGc                                                                   >  2:3883122/1‑150 (MQ=255)
                                                                                 gTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAg                                                                 >  1:3455127/1‑150 (MQ=255)
                                                                                     cTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAgg                                                             >  2:1829918/1‑150 (MQ=255)
                                                                                       gtgAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGGTGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGc                                                           <  2:566859/150‑1 (MQ=255)
                                                                                       gtgAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGc                                                           <  1:1904557/150‑1 (MQ=255)
                                                                                       gtgAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAAATGACCGGTCTGGATCCCGAGCGCGTTCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGc                                                           <  2:2045319/150‑1 (MQ=255)
                                                                                        tgAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGcc                                                          <  1:1560446/150‑1 (MQ=255)
                                                                                         gAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCg                                                         >  2:2320327/1‑150 (MQ=255)
                                                                                         gAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCg                                                         >  1:1688853/1‑150 (MQ=255)
                                                                                         gAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCg                                                         >  2:1054879/1‑150 (MQ=255)
                                                                                         gAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCg                                                         <  1:3990994/150‑1 (MQ=255)
                                                                                          aaTGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGa                                                        >  1:868918/1‑150 (MQ=255)
                                                                                          aaTGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGa                                                        >  1:3508304/1‑150 (MQ=255)
                                                                                            tGGGTGGAAAATAGCATGAGTGCCAATGCACACACCCTGATTTGGATCGCTCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGAcc                                                      <  1:4513457/150‑1 (MQ=255)
                                                                                            tGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGAcc                                                      <  1:4395193/150‑1 (MQ=255)
                                                                                               gTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCAtt                                                   <  1:1045485/150‑1 (MQ=255)
                                                                                                 ggAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGc                                                 <  1:1972786/150‑1 (MQ=255)
                                                                                                  gAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCa                                                <  2:1813665/150‑1 (MQ=255)
                                                                                                  gAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCa                                                <  1:2641861/150‑1 (MQ=255)
                                                                                                   aaaaTAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAg                                               <  1:4121791/150‑1 (MQ=255)
                                                                                                          cATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCa                                        <  2:1926035/150‑1 (MQ=255)
                                                                                                            tGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGt                                      >  2:4334571/1‑150 (MQ=255)
                                                                                                             gAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTc                                     <  1:3595290/150‑1 (MQ=255)
                                                                                                                tGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTCtga                                  <  2:847236/150‑1 (MQ=255)
                                                                                                                  ccAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTCtgatg                                >  1:2748202/1‑150 (MQ=255)
                                                                                                                  ccAATGAAAACAACCTGATTTGGATCGATCTTGAAATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCAGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTCtgatg                                >  2:2777048/1‑150 (MQ=255)
                                                                                                                   cAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTCtgatga                               >  2:2837593/1‑150 (MQ=255)
                                                                                                                    aaTGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTCTGATGaa                              >  2:2548069/1‑150 (MQ=255)
                                                                                                                           acaacCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTCTGATGAACAGCTgg                       >  1:3058631/1‑150 (MQ=255)
                                                                                                                            caacCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTCTGATGAACAGCTGgc                      >  1:1423387/1‑150 (MQ=255)
                                                                                                                              acCTGATTTGGATCGATCTTGAAATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTCTGATGAACAGCTGgcgc                    >  1:2098496/1‑150 (MQ=255)
                                                                                                                                  gATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTCTGATGAACAGCTGGCGCTgat                >  1:449161/1‑150 (MQ=255)
                                                                                                                                   aTTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTCTGATGAACAGCTGGCGCTgatg               >  1:3893948/1‑150 (MQ=255)
                                                                                                                                    tttGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTCTGATGAACAGCTGGCGCTgatgg              >  1:193226/1‑150 (MQ=255)
                                                                                                                                    tttGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTCTGATGAACAGCTGGCGCTgatgg              >  1:1263546/1‑150 (MQ=255)
                                                                                                                                     ttGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTCTGATGAACAGCTGGCGCTgatgga             >  1:309553/1‑150 (MQ=255)
                                                                                                                                      tGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTCTGATGAACAGCTGGCGCTgatggat            <  2:1549064/150‑1 (MQ=255)
                                                                                                                                       ggATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTCTGATGAACAGCTGGCGCTgatggatg           >  2:1661046/1‑150 (MQ=255)
                                                                                                                                        gatcgatcTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTCTGATGAACAGCTGGCGCTGATGGATGa          >  2:3684877/1‑150 (MQ=255)
                                                                                                                                            gatcTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTCTGATGAACAGCTGGCGCTGATGGATGACTgg      <  1:2216343/150‑1 (MQ=255)
                                                                                                                                             atcTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGAGTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTCTGATGAACAGCTGGCGCTGATGGATGACTGGa     <  2:2172918/150‑1 (MQ=255)
                                                                                                                                             atcTTGAGATGACACTACTGGATCCCGAGCGCGATCGCATTATTTAGATTGCCACGCTGGTGACAGATGCCAACCTGAATATTCTTGCAGAAGGGCCGACCATTTCAGTACACCAGTCTGATGAACAGCTGGCGCTGATGCATGACTGGa     >  2:1843317/1‑150 (MQ=255)
                                                                                                                                              tcTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTCTGATGAACAGCTGGCGCTGATGGATGACTGGaa    >  2:2208452/1‑150 (MQ=255)
                                                                                                                                              tcTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTCTGATGAACAGCTGGCGCTGATGGATGACTGGaa    >  2:3352409/1‑150 (MQ=255)
                                                                                                                                                ttGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTCTGATGAACAGCTGGCGCTGATGGATGACTGGAACg  >  1:537092/1‑150 (MQ=255)
                                                                                                                                                ttGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTCTGATGAACAGCTGGCGCTGATGGATGACTGGAACg  >  1:3909932/1‑150 (MQ=255)
                                                                                                                                                ttGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTCTGATGAACAGCTGGCGCTGATGGATGACTGGAACg  >  1:1456006/1‑150 (MQ=255)
                                                                                                                                                ttGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGACAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTCTGATGAACAGCTGGCGCTGATGGATGACTGGAACg  >  1:2466713/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TATTTTTACTCAATCGTACAGACTCCTGGTCGCCCCTGATGGGCAAAACATCTATGATACACGCAATTGTGGATCAATATAGTCACTGTGAATGGGTGGAAAATAGCATGAGTGCCAATGAAAACAACCTGATTTGGATCGATCTTGAGATGACCGGTCTGGATCCCGAGCGCGATCGCATTATTGAGATTGCCACGCTGGTGACCGATGCCAACCTGAATATTCTGGCAGAAGGGCCGACCATTGCAGTACACCAGTCTGATGAACAGCTGGCGCTGATGGATGACTGGAACG  >  CP009273/4381314‑4381607

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