Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 1,915,165 T→C intergenic (+54/+14) RPPX_08535 → / ← RPPX_08540 ferrioxamine B receptor/5S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099741,915,1650TC99.5% 539.0 / ‑6.0 193intergenic (+54/+14)RPPX_08535/RPPX_08540ferrioxamine B receptor/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (75/114);  total (76/117)
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.98e-01

Only 100 of 198 total aligned reads displayed.

CCAATGAAAGCTACATTACCTCATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTGTT  >  CP009974/1915017‑1915304
                                                                                                                                                    |                                                                                                                                           
ccAATGAAAGCTACATTACCTCATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGctt                                                                                                                                           >  2:3031789/1‑151 (MQ=255)
 cAATGAAAGCTACATTACCTCATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGcttt                                                                                                                                          >  1:3945235/1‑151 (MQ=255)
 cAATGAAAGCTACATTACCTCATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGcttt                                                                                                                                          >  1:2591835/1‑151 (MQ=255)
   aTGAAAGCTACATTACCTCATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACTGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTgg                                                                                                                                        <  1:640685/151‑2 (MQ=255)
    tGAAAGCTACATTACCTCATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTggg                                                                                                                                       <  1:3282142/151‑3 (MQ=255)
     gAAAGCTACATTACCTCATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTgggt                                                                                                                                      <  1:3582779/151‑4 (MQ=255)
     gAAAGCTACATTACCTCATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTgggt                                                                                                                                      >  1:2491683/1‑148 (MQ=255)
                tACCTCATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGa                                                                                                                          <  2:1961528/151‑1 (MQ=255)
                   cTCATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGa                                                                                                                       >  1:3800264/1‑151 (MQ=255)
                     cATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATg                                                                                                                     <  2:2361254/151‑1 (MQ=255)
                      aTGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGa                                                                                                                    >  2:4010525/1‑151 (MQ=255)
                       tgtgCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGAc                                                                                                                   >  1:663541/1‑151 (MQ=255)
                       tgtgCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGAc                                                                                                                   <  1:3361328/151‑1 (MQ=255)
                       tgtgCCAGCCTGAACTACTGCTACATGAGTGAAGAGCGCAATGTAAGCGCTACAGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGAc                                                                                                                   <  2:707473/151‑1 (MQ=255)
                          gCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTa                                                                                                                <  2:3901909/151‑1 (MQ=255)
                           ccAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTAc                                                                                                               <  1:1315934/151‑1 (MQ=255)
                            cAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTAct                                                                                                              <  2:1721056/151‑1 (MQ=255)
                                 tGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTcac                                                                                                         >  1:3566073/1‑151 (MQ=255)
                                     ctactGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATgg                                                                                                     <  1:305430/151‑1 (MQ=255)
                                      tactGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATggg                                                                                                    >  2:248800/1‑151 (MQ=39)
                                      tactGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATggg                                                                                                    <  1:2996415/151‑1 (MQ=255)
                                       actGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATgggg                                                                                                   >  1:519225/1‑151 (MQ=39)
                                         tGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGaa                                                                                                 >  2:2234193/1‑151 (MQ=38)
                                         tGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGaa                                                                                                 <  1:1016248/151‑1 (MQ=38)
                                         tGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGaa                                                                                                 <  1:1016973/151‑1 (MQ=38)
                                          gCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGaaa                                                                                                >  1:3263438/1‑151 (MQ=38)
                                          gCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGaaa                                                                                                >  2:1086315/1‑151 (MQ=38)
                                          gCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGaaa                                                                                                >  1:1672424/1‑151 (MQ=38)
                                              cATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAAcccc                                                                                            >  1:3411296/1‑151 (MQ=38)
                                                   gTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTACTAACAGATATTGAAACCAGTAAAAAACAGAGCACCTATCTGAATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACt                                                                                       <  1:339994/151‑1 (MQ=25)
                                                     gAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTAc                                                                                     >  1:2644760/1‑151 (MQ=38)
                                                      aaGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTAcc                                                                                    <  2:3153126/151‑1 (MQ=37)
                                                       agagCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCa                                                                                   <  1:363333/151‑1 (MQ=37)
                                                       agagCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCa                                                                                   <  1:1223660/151‑1 (MQ=37)
                                                         agCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATc                                                                                 <  2:2117424/151‑1 (MQ=37)
                                                           cgcAATGTAAGCGCTAGGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAATCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCgg                                                                               <  2:4202889/151‑1 (MQ=25)
                                                            gcAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGc                                                                              <  1:4895409/151‑1 (MQ=37)
                                                             cAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATTACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCg                                                                             >  1:5434430/1‑151 (MQ=37)
                                                              aaTGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCCCACAGCGGCATTGGGTCTGAATCTTGACGATGACCTACTCTCCCATAGGGAAACCCACCACTACCATCGGCGa                                                                            >  2:3526277/1‑151 (MQ=25)
                                                              aaTGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCACACATGGGGAACCCCCACACTACCATCGGCGa                                                                            <  2:64592/151‑1 (MQ=25)
                                                              aaTGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTAACGATGACCTACTCTCACATGTGGAAACCCCACACTACCATCGGCGa                                                                            >  2:3676279/1‑151 (MQ=25)
                                                               aTGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGAt                                                                           <  1:4257200/151‑1 (MQ=37)
                                                               aTGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGAt                                                                           >  1:2948016/1‑151 (MQ=37)
                                                               aTGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATAGTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGAt                                                                           >  1:3436478/1‑151 (MQ=25)
                                                                tGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATg                                                                          <  2:4635148/151‑1 (MQ=25)
                                                                 gTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGc                                                                         >  2:1400466/1‑151 (MQ=25)
                                                                 gTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGc                                                                         <  1:4589033/151‑1 (MQ=25)
                                                                 gTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGc                                                                         <  1:2155052/151‑1 (MQ=25)
                                                                 gTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCGTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGc                                                                         >  2:341622/1‑151 (MQ=25)
                                                                 gTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGATGTGGGTTTGAAGCTTGACGATGACCTACTCTCACATGGGGAAACCCCAAACTACCATCGGCGATGc                                                                         >  1:532612/1‑151 (MQ=25)
                                                                  tAAGCGCTAAGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCa                                                                        <  1:4967785/151‑1 (MQ=25)
                                                                   aaccgcTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCAt                                                                       <  1:3301559/148‑1 (MQ=25)
                                                                    aGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATc                                                                      <  1:1209381/151‑1 (MQ=25)
                                                                            ggTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCact                                                              >  1:2667094/1‑151 (MQ=25)
                                                                             gTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCactg                                                             >  1:2713593/1‑150 (MQ=21)
                                                                              tCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCactgc                                                            <  1:4048039/151‑3 (MQ=21)
                                                                              tCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGAGTCaatac                                                            >  1:14760/1‑151 (MQ=21)
                                                                                cgCTATCAGTTCTAACCGATATTGAGGCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCActgctg                                                          <  2:4661931/150‑1 (MQ=18)
                                                                                  cTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAg                                                        <  1:2082809/151‑1 (MQ=21)
                                                                                    ctCAGTTCTAACCGATATTGAATCCCGTGAAAGACAACGCCCCTGTCCGTATCGCAGACGGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGtt                                                      <  2:532612/150‑1 (MQ=18)
                                                                                     tCAGTTCTAACTGTTATTGAAGCCCTTGAAAGACAACGCCCCTGTCTGTAACGAAGTCGGGGGCTTTGGCTTTGAATCTTGACGATGGCCTACTCTCACCTGGGGAAACCCCACACTAACATCGGCGGTGCATCGTTTCACTACTGAGTTc                                                     <  2:4316598/151‑1 (MQ=12)
                                                                                        gTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCggg                                                  <  2:1694347/151‑1 (MQ=18)
                                                                                         ttCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCggga                                                 <  1:1788530/151‑1 (MQ=18)
                                                                                         ttCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCggga                                                 <  2:4912311/151‑1 (MQ=18)
                                                                                          tCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGCTcacgat                                                >  1:600377/1‑146 (MQ=18)
                                                                                             aaCCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggatggg                                             >  1:2360072/1‑151 (MQ=18)
                                                                                              aCCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggatggga                                            <  2:2340929/151‑1 (MQ=18)
                                                                                              aCCGAGATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggatggga                                            <  2:5374502/151‑1 (MQ=18)
                                                                                                   tatTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAgg                                       >  2:3865966/1‑151 (MQ=17)
                                                                                                    atTGAACCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAggt                                      <  1:4784086/151‑1 (MQ=17)
                                                                                                    atTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAggt                                      >  1:4465164/1‑151 (MQ=17)
                                                                                                     ttGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAggtg                                     <  2:4110702/151‑1 (MQ=17)
                                                                                                         aaCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTc                                 <  1:1639904/151‑1 (MQ=17)
                                                                                                         aaCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGGTTCACTGCTGGGTTCGTGATTGGATCAGGTGGTTc                                 >  2:3251059/1‑151 (MQ=17)
                                                                                                              gtaaaaGACAAACCCTCTCTCTGCATCGCAGACAGGGGCTTGGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATg                            <  2:2039831/148‑1 (MQ=11)
                                                                                                                gAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGct                          <  1:4624538/151‑1 (MQ=14)
                                                                                                                 aaaGACAAAGCGGCTGTCTGCATCGCAGACAGGGGGTTTGGGTTTGAATCTTGGCGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCCTCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGctc                         <  2:1314043/151‑1 (MQ=14)
                                                                                                                 aaaGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGctc                         <  1:2357614/151‑1 (MQ=14)
                                                                                                                    gACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGATTTCGGTATGTGATCAGGTGGTTCCAATGCTCTAt                      >  2:2735921/1‑151 (MQ=14)
                                                                                                                     accaaGCCCCTGTCTGCGTCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATg                     <  2:5280955/148‑1 (MQ=14)
                                                                                                                       aaaGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGgt                   >  1:5252633/1‑151 (MQ=14)
                                                                                                                       aaaGCCACTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGgt                   <  2:866400/151‑1 (MQ=14)
                                                                                                                         aGCCCCTGTGTGTATCGCAGACAGGGTCTTTCGATTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGgtcg                 <  2:3275413/151‑1 (MQ=14)
                                                                                                                           ccccTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGctc                         <  2:1563405/141‑1 (MQ=14)
                                                                                                                           ccccTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGctc                         >  1:1563405/1‑141 (MQ=14)
                                                                                                                            cccTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCa              >  2:3691203/1‑151 (MQ=14)
                                                                                                                             ccTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCaa             >  1:2819534/1‑151 (MQ=11)
                                                                                                                             ccTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCaa             >  2:4710296/1‑151 (MQ=11)
                                                                                                                              cTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAg            <  1:2225694/151‑1 (MQ=11)
                                                                                                                               tGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGa           >  2:5183072/1‑151 (MQ=11)
                                                                                                                                gTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGaa          <  2:2526173/151‑1 (MQ=11)
                                                                                                                                 tCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGaaa         >  1:3601071/1‑151 (MQ=11)
                                                                                                                                  cTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAAt        >  2:1316820/1‑151 (MQ=11)
                                                                                                                                   cgCATCGCGGGCAGGGGCTTTCGATTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAAtt       <  2:1222810/150‑1 (MQ=11)
                                                                                                                                      aTCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATTGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTg    >  1:3227769/1‑151 (MQ=11)
                                                                                                                                      aTCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCACGAAATTccg    >  2:5272107/1‑149 (MQ=11)
                                                                                                                                      aTCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTg    <  2:1041130/151‑1 (MQ=11)
                                                                                                                                      aTCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTg    <  2:2713593/151‑1 (MQ=11)
                                                                                                                                        cGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTGTg  >  1:713011‑M2/1‑150 (MQ=255)
                                                                                                                                        cGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAggtgg                                    <  2:3156773/117‑1 (MQ=11)
                                                                                                                                                    |                                                                                                                                           
CCAATGAAAGCTACATTACCTCATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTGTT  >  CP009974/1915017‑1915304

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