Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 503,832 A→T 11.1% intergenic (+35/‑116) PP_0414 → / → rpe polyamine ABC transporter permease/ribulose‑5‑phosphate 3‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947503,8320AT11.1% 268.8 / 19.7 99intergenic (+35/‑116)PP_0414/rpepolyamine ABC transporter permease/ribulose‑5‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (40/48);  new base T (5/6);  total (45/54)
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.13e-01

GACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCATTC  >  NC_002947/503686‑503978
                                                                                                                                                   |                                                                                                                                                   
gACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAg                                                                                                                                                   >  7:238087/1‑149 (MQ=255)
    aTTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCg                                                                                                                                                >  1:189276/1‑148 (MQ=255)
       gCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑Gaa                                                                                                                                            <  1:12455/148‑1 (MQ=255)
       gCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑Gaa                                                                                                                                            <  1:335861/148‑1 (MQ=255)
        ccGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTGCCGGCCTCTTC‑C‑GGCTTTAGCCGCGAaga                                                                                                                                          <  7:285886/149‑1 (MQ=39)
         cgcgGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCGGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTGCCGGCCTCTTC‑C‑GGCTTTAGCCGCGAagag                                                                                                                                         <  5:242482/149‑1 (MQ=38)
           cggcggcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTTGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTGCCGGCCTCTTC‑C‑GGCTTTAGCCGCGAagag                                                                                                                                         <  5:139090/147‑1 (MQ=38)
           cggcggcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTGCCGGCCTCTTC‑C‑GGCTTTAGCCGCGAagag                                                                                                                                         >  6:139090/1‑147 (MQ=39)
           cggcggcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTGCCGGCCTCTTC‑C‑GGCTTTAGCCGCGAAGAGGc                                                                                                                                       <  7:54827/149‑1 (MQ=38)
           cggcggcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGc                                                                                                                                       <  8:238087/149‑1 (MQ=255)
            ggcggcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTGCCGGCCTCTTC‑C‑GGCTTTAGCCGCGAAGAGGcc                                                                                                                                      <  4:223794/149‑1 (MQ=38)
             gcggcgACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCg                                                                                                                                     <  7:287646/149‑1 (MQ=255)
                  gACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGcaca                                                                                                                                <  7:199750/149‑1 (MQ=255)
                     cTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGt                                                                                                                             >  5:129630/1‑149 (MQ=255)
                      ttGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGtt                                                                                                                            >  3:247203/1‑149 (MQ=255)
                          tGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGt                                                                                                                        >  1:242508/1‑149 (MQ=255)
                           gATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGtt                                                                                                                       >  6:264100/1‑149 (MQ=255)
                           gATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGT‑GGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTg                                                                                                                      <  5:64487/149‑1 (MQ=255)
                             ttGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTgg                                                                                                                     >  7:185280/1‑149 (MQ=255)
                                cGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTc                                                                                                                   <  8:8383/148‑1 (MQ=255)
                                cGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTc                                                                                                                   <  2:74971/148‑1 (MQ=255)
                                cGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTc                                                                                                                   <  2:189276/148‑1 (MQ=255)
                                 gTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGt                                                                                                                        >  4:85307/1‑142 (MQ=255)
                                 gTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGt                                                                                                                        <  3:85307/142‑1 (MQ=255)
                                   tctcGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACCGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGaca                                                                                                               >  6:81452/1‑149 (MQ=255)
                                    ctcGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGaca                                                                                                               <  3:131605/148‑1 (MQ=255)
                                     tcGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGACAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGacac                                                                                                              >  6:1196/1‑148 (MQ=255)
                                        gtgtTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑Gaa                                                                                                                                            <  4:264178/115‑1 (MQ=255)
                                        gtgtTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑Gaa                                                                                                                                            >  3:264178/1‑115 (MQ=255)
                                        gtgtTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAAt                                                                                                          >  3:83541/1‑149 (MQ=255)
                                        gtgtTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAAt                                                                                                          >  5:47993/1‑149 (MQ=255)
                                        gtgtTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAAt                                                                                                          >  6:10961/1‑149 (MQ=255)
                                          gttgttGCTTCTGACCGTGGAGTGGTTGGGTGGGCGTAGCGATAAGCTGAGGACTCAACAGCCGGCTTGATGTTGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑TAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAAtt                                                                                                         <  2:242508/148‑1 (MQ=255)
                                           ttgttgCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGa                                                                                                       <  4:223088/149‑1 (MQ=255)
                                             gttgCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGAt                                                                                                      <  4:83541/148‑1 (MQ=255)
                                              ttgCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATa                                                                                                     <  7:48544/148‑1 (MQ=255)
                                              ttgCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATAc                                                                                                    <  3:25648/149‑1 (MQ=255)
                                              ttgCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATAc                                                                                                    <  5:10961/149‑1 (MQ=255)
                                               tgctgcTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATAc                                                                                                    <  2:85916/148‑1 (MQ=255)
                                                gctgctGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATAcc                                                                                                   >  5:163396/1‑148 (MQ=255)
                                                   gctgACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGa                                                                                                       >  6:90665/1‑141 (MQ=255)
                                                   gctgACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGa                                                                                                       <  5:90665/141‑1 (MQ=255)
                                                            ggAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGAt                                                                                                      <  1:229977/133‑1 (MQ=255)
                                                            ggAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGAt                                                                                                      >  2:229977/1‑133 (MQ=255)
                                                             gAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATc                                                                                     <  6:47993/149‑1 (MQ=255)
                                                                  gTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCtt                                                                                 <  1:304560/148‑1 (MQ=255)
                                                                      cGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTccc                                                                             <  6:163396/148‑1 (MQ=255)
                                                                       gTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGAt                                                                                                      >  5:188887/1‑122 (MQ=255)
                                                                       gTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGAt                                                                                                      <  6:188887/122‑1 (MQ=255)
                                                                       gTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCt                                                                            >  6:87032/1‑148 (MQ=255)
                                                                        tGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGa                                                                          <  4:247203/149‑1 (MQ=255)
                                                                              tAGCGAGAAGCTGAGGACTCAACAGCCTGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACa                                                                    <  8:272257/149‑1 (MQ=255)
                                                                                gCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAAt                                                                  <  6:129630/149‑1 (MQ=255)
                                                                                         tGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCAc                                                         >  1:112188/1‑149 (MQ=255)
                                                                                          gAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCAcc                                                        >  3:223863/1‑149 (MQ=255)
                                                                                           aGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCAc                                                         <  5:264100/147‑1 (MQ=255)
                                                                                                tCAACAGCCGGTTTGATGATGGCTTGCCTGTACCGGCCTCTTCGCGGGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCCACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGaa                                                   >  1:250052/1‑149 (MQ=255)
                                                                                                tCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCg                                                                                                                         <  3:115628/78‑1 (MQ=255)
                                                                                                tCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCg                                                                                                                         >  4:115628/1‑78 (MQ=255)
                                                                                                tCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGGGTTACAATgc                                                                <  3:228099/135‑1 (MQ=255)
                                                                                                tCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGGGTTACAATgc                                                                >  4:228099/1‑135 (MQ=255)
                                                                                                tCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAAt                                                  >  3:255874/1‑149 (MQ=255)
                                                                                                  aaCAGCGGGCTTGATGATGGATTGCCGGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGATGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAAttt                                                <  4:255874/149‑1 (MQ=255)
                                                                                                  aaCAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAAttt                                                <  5:1196/149‑1 (MQ=255)
                                                                                                  aaCAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAAt                                                  <  3:98683/147‑1 (MQ=255)
                                                                                                  aaCAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCAGAGTTACAATtcgc                                                              >  4:338670/1‑135 (MQ=255)
                                                                                                    cAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGc                                              <  2:112188/149‑1 (MQ=255)
                                                                                                       ccGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCaa                                           >  7:182620/1‑149 (MQ=255)
                                                                                                         ggCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTcc                                                                              <  3:168215/112‑1 (MQ=255)
                                                                                                         ggCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTcc                                                                              >  4:168215/1‑112 (MQ=255)
                                                                                                          gtttGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAg                                        <  2:250052/147‑1 (MQ=255)
                                                                                                                tgatgGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTgcgc                                  >  7:239131/1‑149 (MQ=255)
                                                                                                                tgatgGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGCGTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTgcgc                                  >  4:100143/1‑149 (MQ=255)
                                                                                                                tgatgGCTTGCCTGTACCGGCCCCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTgcgc                                  >  8:113550/1‑149 (MQ=255)
                                                                                                                  atgGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCa                                >  7:234384/1‑149 (MQ=255)
                                                                                                                     gCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGc                             <  8:182620/149‑1 (MQ=255)
                                                                                                                      cTTGCCTGTACCGGCCTCTTCGC‑GCCTTAAGCCG‑GAAGAGGCCGGCCCAGGTTTGGTTGGTGGACACCAATTGATACGAATCAGCGATCATCGCATTCCCTGAGTTAAAAAGGCTGGCACCGTGAATTTGCAAACAGGTGCGCCAtt                              >  8:209596/1‑146 (MQ=255)
                                                                                                                       ttGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTgcg                                   <  3:302843/141‑1 (MQ=255)
                                                                                                                       ttGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTgcg                                   >  4:302843/1‑141 (MQ=255)
                                                                                                                       ttGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAg                           >  8:8348/1‑149 (MQ=255)
                                                                                                                          ccTGTACCGGCCTCTTCGC‑GGCTAAAGCAG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGccc                        >  1:19871/1‑149 (MQ=255)
                                                                                                                           cTGTGCCGGCCTCTTC‑C‑GGCTTTAGCCGCGAAGAGGCCGGTACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAg                                        <  3:132768/132‑1 (MQ=37)
                                                                                                                           cTGTGCCGGCCTCTTC‑C‑GGCTTTAGCCGCGAAGAGGCCGGTACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAg                                        >  4:132768/1‑132 (MQ=37)
                                                                                                                           cTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCAcc                                                        <  4:121356/116‑1 (MQ=255)
                                                                                                                           cTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCAcc                                                        >  3:121356/1‑116 (MQ=255)
                                                                                                                              tACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACg                    <  8:239131/149‑1 (MQ=255)
                                                                                                                                 cGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCAt                               <  2:221655/135‑1 (MQ=255)
                                                                                                                                 cGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTa                 >  1:345856/1‑149 (MQ=255)
                                                                                                                                   gCCTCTTC‑C‑GGCTTTAGCCGCGAAGAGGCCGGTACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTAtt               >  5:136100/1‑149 (MQ=38)
                                                                                                                                    ccTCTTC‑C‑GGCTTTAGCCGCGAAGAGGCCGGTACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTg              >  4:249799/1‑149 (MQ=39)
                                                                                                                                     ctattCGC‑GGATAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGATGGTCGACACCAATTGATACCAATCACCCTTCATCGCTTACCCTGAGTTACAATGCGTGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGc             <  7:209596/146‑1 (MQ=255)
                                                                                                                                         ttcc‑gGCTTTAGCCGCGAAGAGGCCGGTACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGcccc          >  5:238365/4‑149 (MQ=255)
                                                                                                                                          cgc‑gGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCTTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGccccc         <  8:185280/148‑1 (MQ=255)
                                                                                                                                          cgc‑gGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGcc                         <  1:27099/132‑1 (MQ=255)
                                                                                                                                          cgc‑gGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGcc                         >  2:27099/1‑132 (MQ=255)
                                                                                                                                          cgc‑gGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGccccc         <  7:113550/148‑1 (MQ=255)
                                                                                                                                           gc‑gGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTc       <  4:223863/149‑1 (MQ=255)
                                                                                                                                              ggCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCa     <  5:87032/149‑1 (MQ=255)
                                                                                                                                                 tAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCATtc  >  3:72395/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGACCCTGGAGTGGTTGCGTGGGCGTAGCGAGAAGCTGAGGACTCAACAGCCGGCTTGATGATGGCTTGCCTGTACCGGCCTCTTCGC‑GGCTAAAGCCG‑GAAGAGGCCGGCACAGGTTGCGTTGGTCGACACCAATTGATACCAATCAGCCTTCATCGCTTTCCCTGAGTTACAATGCGCGCCACCGTGAATTTGCCAACAGGTGCGCCATGCAGCCCTACGCTATTGCCCCCTCCATTC  >  NC_002947/503686‑503978

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