Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 1,915,172 3 bp→GTT intergenic (+61/+5) 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,1720AG83.0% 399.5 / 122.2 176intergenic (+61/+7)RPPX_08535/RPPX_08540ferrioxamine B receptor/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (76/70);  minor base . (10/20);  total (86/90)
Fisher's exact test for biased strand distribution p-value = 7.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01
*CP0099741,915,1730AT99.4% 486.1 / NA 176intergenic (+62/+6)RPPX_08535/RPPX_08540ferrioxamine B receptor/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (85/90);  minor base C (1/0);  total (86/90)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*CP0099741,915,1740AT99.4% 488.5 / NA 176intergenic (+63/+5)RPPX_08535/RPPX_08540ferrioxamine B receptor/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (85/90);  total (86/90)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 176 total aligned reads displayed.

ACCTCATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTGT  >  CP009974/1915034‑1915303
                                                                                                                                          |||                                                                                                                                 
aCCTCATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGAc                                                                                                                        >  2:1232937/1‑151 (MQ=255)
    cATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATg                                                                                                                    <  1:618365/151‑1 (MQ=255)
    cATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATg                                                                                                                    <  1:513304/151‑1 (MQ=255)
        tgCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCt                                                                                                                >  1:1283766/1‑151 (MQ=255)
        tgCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCt                                                                                                                <  1:3870131/151‑1 (MQ=255)
         gCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTa                                                                                                               <  1:4461359/151‑1 (MQ=255)
          ccAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAGGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTAc                                                                                                              <  2:4434909/151‑1 (MQ=255)
            aGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTActc                                                                                                            >  1:848319/1‑151 (MQ=255)
             gCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTActct                                                                                                           >  2:1953036/1‑151 (MQ=255)
             gCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTActct                                                                                                           <  2:1004034/151‑1 (MQ=255)
               cTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTca                                                                                                         >  1:1295140/1‑151 (MQ=255)
                tGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTcac                                                                                                        >  1:4731854/1‑151 (MQ=255)
                tGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTcac                                                                                                        >  1:912075/1‑151 (MQ=255)
                       ctGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGa                                                                                                 <  2:1496916/151‑1 (MQ=38)
                       ctGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATCGGGa                                                                                                 >  2:4042488/1‑151 (MQ=255)
                         ggtACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGaaa                                                                                               <  2:3177736/149‑1 (MQ=38)
                          cTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAAc                                                                                              <  1:3199114/151‑1 (MQ=38)
                           tACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAAcc                                                                                             >  2:2994147/1‑151 (MQ=38)
                             cATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAAcccc                                                                                           <  2:2952975/151‑1 (MQ=38)
                             cATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAAcccc                                                                                           <  2:1783400/151‑1 (MQ=38)
                               tGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCcac                                                                                         >  1:3413774/1‑151 (MQ=38)
                               tGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATCGGGAAGCCCcac                                                                                         >  2:4455340/1‑151 (MQ=38)
                               gggggGAAGGGCGCAATGGAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCGCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCcac                                                                                         <  2:3277292/150‑1 (MQ=37)
                                    gAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCGTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTAc                                                                                    >  2:5198101/1‑151 (MQ=37)
                                       gagCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCAt                                                                                 <  1:5036586/151‑1 (MQ=37)
                                         gcgcAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCg                                                                               <  1:3115000/151‑1 (MQ=37)
                                         gcgcAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCg                                                                               <  2:1723338/151‑1 (MQ=37)
                                           gcAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGc                                                                             <  2:2653078/151‑1 (MQ=37)
                                           gcAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGc                                                                             >  1:779145/1‑151 (MQ=37)
                                           gcAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCACCGGc                                                                             >  2:1677487/1‑151 (MQ=37)
                                            cacggTAAGCGCTCCGGTCAGCTATCAGTTCTAACCGATATTGAAACACGTGAAAGACAAAGCCTCTTTCTGCATTGAAGACAGGGTCTTAGGGTTTGAATCTTGAGGATGCTGTACTCTCTCCTGCGAAAATTCCACACTGCCAGCcgcg                                                                            >  2:2774495/5‑151 (MQ=255)
                                             aaTGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGa                                                                           <  1:1592049/151‑1 (MQ=37)
                                              aTGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGAt                                                                          >  1:2204045/1‑151 (MQ=37)
                                              aTGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGAt                                                                          <  1:1662856/151‑1 (MQ=37)
                                                     cgcTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGt                                                                   <  2:4526758/151‑1 (MQ=25)
                                                     cgcTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGt                                                                   >  2:4984520/1‑151 (MQ=25)
                                                      gcTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGtt                                                                  <  1:473573/151‑1 (MQ=25)
                                                       cTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGttt                                                                 <  2:678060/151‑1 (MQ=25)
                                                       cTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGttt                                                                 >  1:1885070/1‑151 (MQ=25)
                                                       cTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGttt                                                                 <  1:152739/151‑1 (MQ=25)
                                                        tACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTc                                                                >  2:2020101/1‑151 (MQ=25)
                                                         aCGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCa                                                               >  1:4632998/1‑151 (MQ=21)
                                                         aCGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCa                                                               >  1:377110/1‑151 (MQ=25)
                                                           ggTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCact                                                             <  2:2011764/151‑1 (MQ=21)
                                                           ggTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCact                                                             <  1:1166881/151‑1 (MQ=21)
                                                           ggTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCact                                                             >  1:3570733/1‑151 (MQ=25)
                                                            gTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCacta                                                            >  1:2593894/1‑151 (MQ=21)
                                                             tCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGTCCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCactgc                                                           <  2:5046621/151‑3 (MQ=21)
                                                             tCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCactac                                                           <  2:1218770/151‑1 (MQ=21)
                                                               aGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCActgctg                                                         <  2:3047724/151‑1 (MQ=21)
                                                                gCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGa                                                        >  2:4719881/1‑151 (MQ=21)
                                                                  tATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGt                                                      <  2:2565155/151‑1 (MQ=21)
                                                                   aTCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGtt                                                     >  1:4697500/1‑151 (MQ=21)
                                                                   aTCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGtt                                                     <  2:906068/151‑1 (MQ=21)
                                                                    tCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTc                                                    <  2:5107954/151‑1 (MQ=21)
                                                                     cAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCg                                                   <  2:1871116/151‑1 (MQ=21)
                                                                      aGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgg                                                  >  2:2561087/1‑151 (MQ=21)
                                                                      aGTTCGAACCGCTATTGAAACCCGTGAAATACAAAGGCGCTGTGTGCATCGCAGACAGGGGCTGTGGGTATGAATCTTGAAGATGACCTACTCACACAGGGGGAAACCCCACACTACCATCGGCGATACATCGTTTCACTATTGCGCTCgg                                                  >  2:4121604/1‑151 (MQ=17)
                                                                        ttCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCggga                                                >  1:4980868/1‑151 (MQ=18)
                                                                           tAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggatgg                                             >  1:3209776/1‑151 (MQ=18)
                                                                            aaCCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCCGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCgggatggg                                            <  2:4569838/151‑1 (MQ=18)
                                                                               ctatatTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGACATCGGGGTGGTATc                                         >  1:1346947/3‑151 (MQ=18)
                                                                               ctatatTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTCAGTTCTGGATGGTATc                                         >  1:1347251/3‑151 (MQ=18)
                                                                               cGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATc                                         <  2:4682660/151‑1 (MQ=18)
                                                                                gATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGAGCa                                        >  1:2127244/1‑151 (MQ=18)
                                                                                 atatTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAg                                       <  1:2761973/151‑1 (MQ=18)
                                                                                  tatTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAgg                                      >  1:4581941/1‑151 (MQ=18)
                                                                                   cctgcaaCCCGTGAAAGACAAAGCCCCTCTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACGCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAggt                                     <  1:3794996/146‑1 (MQ=14)
                                                                                   atTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAggt                                     >  2:417320/1‑151 (MQ=17)
                                                                                        aaCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTc                                <  1:436911/151‑1 (MQ=17)
                                                                                        aaCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCact                                                             <  1:2547496/122‑1 (MQ=17)
                                                                                        aaCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGCTCAGGTGGTTc                                >  2:1519210/1‑151 (MQ=17)
                                                                                         aCCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGTGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTcc                               >  1:1625215/1‑151 (MQ=17)
                                                                                            cGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAAt                            >  1:1451994/1‑151 (MQ=17)
                                                                                             gTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACTATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATg                           >  1:2419672/1‑151 (MQ=14)
                                                                                               gAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGct                         >  1:4539503/1‑151 (MQ=14)
                                                                                                 aaGCCATAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCGCATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGctct                       <  2:1228599/151‑1 (MQ=14)
                                                                                                 aaGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGctct                       >  2:3768484/1‑151 (MQ=14)
                                                                                                  aGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCcaca                                                                                        <  2:3342569/85‑1 (MQ=17)
                                                                                                  aGACAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCcaca                                                                                        >  1:3342569/1‑85 (MQ=18)
                                                                                                    aCAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATg                    <  2:922014/151‑1 (MQ=14)
                                                                                                     cAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATgg                   <  2:611540/151‑1 (MQ=14)
                                                                                                     cAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATgg                   >  1:544301/1‑151 (MQ=14)
                                                                                                       aaGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGgtc                 <  1:230514/151‑1 (MQ=14)
                                                                                                        aGCCCCTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGgtcg                <  2:3612983/151‑1 (MQ=14)
                                                                                                           cccTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCa             <  1:2752445/151‑1 (MQ=14)
                                                                                                           cccTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCa             <  2:779987/151‑1 (MQ=14)
                                                                                                           cacTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCa             <  2:5112398/149‑1 (MQ=14)
                                                                                                            ccTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCaa            <  1:2930093/151‑1 (MQ=14)
                                                                                                            ccTGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCaa            <  2:3414013/151‑1 (MQ=14)
                                                                                                              tGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGa          >  1:1035301/1‑151 (MQ=11)
                                                                                                              tGTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGgtc                 >  1:183351/1‑144 (MQ=11)
                                                                                                               gTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCa             >  2:4251452/1‑147 (MQ=11)
                                                                                                               gTCTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCa             <  1:4251452/147‑1 (MQ=11)
                                                                                                                 cTGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAAt       >  1:3255914/1‑151 (MQ=11)
                                                                                                                  tGCATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAAtt      >  2:2286624/1‑151 (MQ=11)
                                                                                                                    cATCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCt    <  2:750341/151‑1 (MQ=11)
                                                                                                                     aTCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTg   >  1:3711624/1‑151 (MQ=11)
                                                                                                                      tCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTGt  <  1:3535669/151‑1 (MQ=11)
                                                                                                                      tCGCAGACAGGGGCTTTGGGTTTGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTGCTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTGt  >  1:2012034/1‑151 (MQ=11)
                                                                                                                                          |||                                                                                                                                 
ACCTCATGTGCCAGCCTGAACTACTGCTACATGGGTGAAGAGCGCAATGTAAGCGCTACGGTCAGCTATCAGTTCTAACCGATATTGAAACCCGTGAAAGACAAAGCCCCTGTCTGCATCGCAGACAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGATGCATCGTTTCACTACTGAGTTCGGGATGGGATCAGGTGGTTCCAATGCTCTATGGTCGTCAAGAAATTCTGT  >  CP009974/1915034‑1915303

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