Predicted mutation
evidence seq id position mutation annotation gene description
RA CP043523 1,867,833 T→G L103L (CTT→CTG FTN78_chr1914 → Site‑specific recombinase, DNA invertase Pin related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0435231,867,8330TG99.8% 1537.1 / ‑4.9 435L103L (CTT→CTGFTN78_chr1914Site‑specific recombinase, DNA invertase Pin related protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (257/177);  total (258/177)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 435 total aligned reads displayed.

GGCGATACTTTAATTGTGACAAAGTTAGATCGATTAGCCAGGAACACTAAAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTTATATATTAAATATTGGACTTATCGATAATACACCGACAGGTCAACTTATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAACTTTGCTAAAATCCATGATCCGAACTTTAAAGAAGGACGTCCTCAGAAATTTACTGAAGAACAAATTC  >  CP043523/1867687‑1867975
                                                                                                                                                  |                                                                                                                                              
ggtgacaCTTTAATTGTGACAAAGTTAGACCGATTAGCCAGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGAtt                                                                                                                                             >  1:1596278/7‑150 (MQ=255)
  tgacaCTTTAATTGTGACAAAGTTAGACCGATTAGCCAGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGAtttt                                                                                                                                           <  2:1425214/146‑1 (MQ=255)
  tgacaCTTTAATTGTGACAAAGTTAGACCGATTAGCCAGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGAtttt                                                                                                                                           <  1:805684/146‑1 (MQ=255)
  tgacaCTTTAATTGTGACAAAGTTAGACCGATTAGCCAGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGAtttt                                                                                                                                           >  1:198695/5‑150 (MQ=255)
    acaCTTTAATTGTGACAAAGTTAGACCGATTAGCCAGGAACACTAGAGAAGCTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTa                                                                                                                                         >  2:4607835/3‑150 (MQ=255)
         ttAATTGTGACAAAGTTAGACCGATTAGCCAGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATATAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATc                                                                                                                                    >  2:4244257/1‑150 (MQ=255)
         ttAATTGTGACAAAGTTAGACCGATTAGCCAGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATc                                                                                                                                    <  1:2791374/150‑1 (MQ=255)
               gtgACAAAGTTAGACCGATTAGCCAGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGt                                                                                                                              >  2:4428403/1‑150 (MQ=255)
               gtgACAAAGTTAGACCGATTAGCCAGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGt                                                                                                                              <  1:1316043/150‑1 (MQ=255)
               gtgACAAAGTTAGACCGATTAGCCAGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGt                                                                                                                              <  1:2197499/150‑1 (MQ=255)
                 gACAAAGTTAGACCGATTAGCCAGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGc                                                                                                                            <  2:3390027/150‑1 (MQ=255)
                  aCAAAGTTAGACCGATTAGCCAGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCt                                                                                                                           <  1:1735493/150‑1 (MQ=255)
                   cAAAGTTAGACCGATTAGCCAGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATTCTAAATATGGGTATTATTGATAACACTCCCTCAGGTCAGCTTATCTTTACGATTTTTAGTGcat                                                                                                                          >  2:2065535/1‑148 (MQ=255)
                     aaGTTAGACCGATTAGCCAGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCtttt                                                                                                                        <  1:530508/150‑1 (MQ=255)
                      aGTTAGACCGATTAGCCAGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTg                                                                                                                       >  2:1319063/1‑150 (MQ=255)
                      aGTTAGACCGATTAGCCAGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTg                                                                                                                       >  2:1702992/1‑150 (MQ=255)
                       gTTAGACCGATTAGCCAGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTgc                                                                                                                      >  1:2948916/1‑150 (MQ=255)
                         tAGACCGATTAGCCAGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTgcgc                                                                                                                    >  2:2850373/1‑150 (MQ=255)
                          aGACCGATTAGCCAGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCa                                                                                                                   >  1:398438/1‑150 (MQ=255)
                               gATTAGCCAGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTg                                                                                                              >  2:1481263/1‑150 (MQ=255)
                                 ttAGCCAGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATACTACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGaa                                                                                                            <  1:3411536/150‑1 (MQ=255)
                                  tAGCCAGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAAc                                                                                                           <  2:392391/150‑1 (MQ=255)
                                     ccAGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACgag                                                                                                        <  2:3115261/150‑1 (MQ=255)
                                      cAGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACgaga                                                                                                       >  1:1671831/1‑150 (MQ=255)
                                      cAGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACgaga                                                                                                       >  2:2201174/1‑150 (MQ=255)
                                       aGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGCATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGAc                                                                                                      >  2:3764570/1‑150 (MQ=255)
                                       aGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGAc                                                                                                      >  2:2616501/1‑150 (MQ=255)
                                       aGGAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGAc                                                                                                      <  2:2266637/150‑1 (MQ=255)
                                        ggAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGAcc                                                                                                     <  2:216041/150‑1 (MQ=255)
                                         gAACACTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCt                                                                                                    >  1:3001625/1‑150 (MQ=255)
                                             acTAGAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTAtt                                                                                                >  1:2411332/1‑150 (MQ=255)
                                                  agaagTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAAc                                                                                           >  2:3251044/1‑150 (MQ=255)
                                                   gaagTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAAcc                                                                                          >  1:2705473/1‑150 (MQ=255)
                                                    aagTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATATTCAGTGCTTTTGCGCAGTTTGGACGAGAGCTGATTGTAACCa                                                                                         >  1:4366510/1‑150 (MQ=255)
                                                     agTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAg                                                                                        >  1:2861112/1‑150 (MQ=255)
                                                      gTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAgg                                                                                       <  1:3128440/150‑1 (MQ=255)
                                                      gTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAgg                                                                                       >  2:404290/1‑150 (MQ=255)
                                                       ttttAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGa                                                                                      <  1:152241/150‑1 (MQ=255)
                                                        tttAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGAc                                                                                     >  1:1852170/1‑150 (MQ=255)
                                                        tttAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGAc                                                                                     <  2:774946/150‑1 (MQ=255)
                                                          tAGAAATTGTTCAATCGTTATTTAATAGAGGTATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAAACAGGACTc                                                                                   >  1:3839789/1‑150 (MQ=255)
                                                             aaaTTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCagg                                                                                <  2:576693/150‑1 (MQ=255)
                                                               aTTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCaggaa                                                                              >  2:3515917/1‑150 (MQ=255)
                                                                ttGTTCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCaggaag                                                                             >  1:2648889/1‑150 (MQ=255)
                                                                   ttCAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAggaaggaa                                                                          >  1:3274935/1‑150 (MQ=255)
                                                                     cAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGaaaa                                                                        <  1:2460155/150‑1 (MQ=255)
                                                                     cAATCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGaaaa                                                                        <  1:1006156/150‑1 (MQ=255)
                                                                     cAATCGTTATTTAATAGAGGAATTAAGGTTCAGATATTAAATACTGGTCTTATTGATAATAGACCCACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGaaaa                                                                        <  2:4123122/150‑1 (MQ=255)
                                                                       aTCGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGaaaaaa                                                                      <  1:3102302/150‑1 (MQ=255)
                                                                         cGTTATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGaaaaaatt                                                                    >  2:3680644/1‑148 (MQ=255)
                                                                           ttATTTAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAAtttt                                                                  >  2:503918/1‑150 (MQ=255)
                                                                               ttAATAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTa                                                              >  1:1236814/1‑150 (MQ=255)
                                                                                   tAGAGGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTaagct                                                          <  2:874654/150‑4 (MQ=255)
                                                                                    agagGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTaagctt                                                         <  2:196351/150‑5 (MQ=255)
                                                                                    agagGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTaagctt                                                         <  2:2897239/150‑5 (MQ=255)
                                                                                     gagGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCt                                                               >  2:2292620/1‑143 (MQ=255)
                                                                                      agGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTaagcttca                                                       >  1:1619709/1‑144 (MQ=255)
                                                                                      agGAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTaagcttca                                                       >  2:2426555/1‑144 (MQ=255)
                                                                                       ggAATTAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTaagcttcan                                                      <  1:2530922/150‑8 (MQ=255)
                                                                                            tAAGGTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATcc                                                 >  2:1003242/1‑150 (MQ=255)
                                                                                               ggTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATccaag                                              >  2:1459826/1‑147 (MQ=255)
                                                                                               ggTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATccaag                                              <  1:2185976/150‑4 (MQ=255)
                                                                                                gTTCATATATTAAATATTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATccaagt                                             >  1:3458240/1‑146 (MQ=255)
                                                                                                             atatTGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTaaggaagggc                                >  1:2279320/1‑148 (MQ=255)
                                                                                                                ttGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTc                             >  2:1906786/1‑150 (MQ=255)
                                                                                                                ttGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTc                             >  1:3620262/1‑150 (MQ=255)
                                                                                                                 tGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTcc                            >  1:3519243/1‑150 (MQ=255)
                                                                                                                 tGGTCTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTcc                            >  2:1219777/1‑150 (MQ=255)
                                                                                                                  ggtcTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCt                           >  1:639493/4‑150 (MQ=255)
                                                                                                                   gtcTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTc                          <  2:3578598/148‑1 (MQ=255)
                                                                                                                   gtcTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAAGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTc                          <  2:712051/148‑1 (MQ=255)
                                                                                                                    tcTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCa                         >  2:841495/2‑150 (MQ=255)
                                                                                                                    tcTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCa                         <  1:4607835/149‑1 (MQ=255)
                                                                                                                     cTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCaa                        >  1:882729/1‑149 (MQ=255)
                                                                                                                     cTTATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCaa                        >  2:1549814/1‑149 (MQ=255)
                                                                                                                      ttATTGATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCaaa                       <  1:4317181/150‑3 (MQ=255)
                                                                                                                        attgATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCAaaaa                     >  1:687001/4‑150 (MQ=255)
                                                                                                                        attgATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCAaaaa                     >  1:1807140/4‑150 (MQ=255)
                                                                                                                          tgATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCAAAAAtt                   <  1:2074386/149‑1 (MQ=255)
                                                                                                                           gATAATACACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCAAAAAttt                  <  1:3779407/150‑1 (MQ=255)
                                                                                                                            ataataCACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCAAAAATTTa                 <  2:4068768/150‑1 (MQ=255)
                                                                                                                            ataataCACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTGTTAAGGAAGGGCGTCCTCAAAACTTTa                 >  2:2953860/1‑150 (MQ=255)
                                                                                                                              aataCACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCAAAAATTTAcc               <  1:1341525/150‑2 (MQ=255)
                                                                                                                                taCACCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCAAAAATTTACCga             <  1:3680644/150‑1 (MQ=255)
                                                                                                                                 acacCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCAAAAATTTACCgaa            <  1:4244257/150‑1 (MQ=255)
                                                                                                                                 acacCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCAAAAATTTACCgaa            >  1:584/1‑150 (MQ=255)
                                                                                                                                 acacCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCAAAAATTTACCgaa            >  2:1135991/1‑150 (MQ=255)
                                                                                                                                 acacCTACTGGTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAACTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCAAAAATTTACCgaa            >  1:2336264/1‑150 (MQ=255)
                                                                                                                                   acctactggTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCAAAAATTTACCgaaga          >  2:2725279/8‑150 (MQ=255)
                                                                                                                                   acctactggTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCAAAAATTTACCgaaga          <  1:1457459/143‑1 (MQ=255)
                                                                                                                                    cctactggTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCAAAAATTTACCgaagaa         >  2:2387750/7‑150 (MQ=255)
                                                                                                                                      tactggTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCAAAAATTTACCGAAGAACa       >  1:756307/5‑150 (MQ=255)
                                                                                                                                      tactggTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCAAAAATTTACCGAAGAACa       <  1:1155068/146‑1 (MQ=255)
                                                                                                                                       actggTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCAAAAATTTACCGAAGAACaa      >  2:1984733/4‑150 (MQ=255)
                                                                                                                                       actggTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCAAAAATTTACCGAAGAACaa      >  1:1126306/4‑150 (MQ=255)
                                                                                                                                        ctggTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCAAAAATTTACCGAAGAACaaa     >  1:4531499/3‑150 (MQ=255)
                                                                                                                                         tggTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCAAAAATTTACCGAAGAACAAAt    >  2:1633679/2‑150 (MQ=255)
                                                                                                                                          ggTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCAAAAATTTACCGAAGAACAAAtt   <  2:4336209/150‑1 (MQ=255)
                                                                                                                                          ggTCAGCTGATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCAAAAATTTACCGAAGAACAAAtt   >  1:4199862/1‑150 (MQ=255)
                                                                                                                                           gTCAGCTGATTTTTACCATCTTAAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACAAAGAAAAAAAAAAATTTTGCTAAGCTTCATGATCCAAGTTTTAAGGAAGGGCGTCCTCAAAAATTTACCGAAGAACAAATTc  <  2:644462/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
GGCGATACTTTAATTGTGACAAAGTTAGATCGATTAGCCAGGAACACTAAAGAAGTTTTAGAAATTGTTCAATCGTTATTTAATAGAGGAATTAAGGTTTATATATTAAATATTGGACTTATCGATAATACACCGACAGGTCAACTTATTTTTACCATCTTCAGTGCTTTTGCGCAGTTTGAACGAGACCTTATTGTAACCAGGACTCAGGAAGGAAAAAACTTTGCTAAAATCCATGATCCGAACTTTAAAGAAGGACGTCCTCAGAAATTTACTGAAGAACAAATTC  >  CP043523/1867687‑1867975

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: