Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,071,143 T→A 5.9% intergenic (+95/+132) PP_1846 → / ← PP_1847 group II intron‑encoding maturase/TonB‑dependent ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,071,1430TA 5.9% 195.9 / 2.6 68intergenic (+95/+132)PP_1846/PP_1847group II intron‑encoding maturase/TonB‑dependent ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (32/31);  new base A (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

GGGCTGGTCTCGATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCG  >  NC_002947/2070998‑2071288
                                                                                                                                                 |                                                                                                                                                 
gggCTGGTCTCGATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCa                                                                                                                                                <  2:328795/149‑1 (MQ=255)
      gtttcGATTCTGGATACGATTAATCGGCTTAGCCGCGTAACCTTAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGATTCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCa                                                                                                                                           <  2:72453/145‑1 (MQ=17)
          cGATACTGGATACGATTAACCGGCTTAGCCGTGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCACAGCAGccccaatgg                                                                                                                                      >  6:227280/1‑143 (MQ=14)
           gATACTGGATACGATTAACCGGCTTAGCCGTGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCACAGCAGccccaatgg                                                                                                                                      <  5:227280/148‑7 (MQ=14)
                tGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAg                                                                                                                                               <  6:50948/134‑1 (MQ=255)
                tGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAg                                                                                                                                               >  5:50948/1‑134 (MQ=255)
                   aTACGATTAACCGGCTTAGCCGCGAAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCG                                                                                                                             >  6:315010‑M1/1‑149 (MQ=255)
                     aCGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGG                                                                                                                           >  3:95290‑M1/1‑149 (MQ=255)
                       gATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCa                                                                                                                                                <  4:306935/126‑1 (MQ=255)
                       gATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCa                                                                                                                                                >  3:306935/1‑126 (MQ=255)
                       gATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCT                                                                                                                         >  1:103908‑M1/1‑149 (MQ=255)
                          tAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGC                                                                                                                      >  5:105764‑M1/1‑149 (MQ=255)
                           aaCCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCAGGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGAATTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGC                                                                                                                      <  1:257494‑M1/148‑1 (MQ=255)
                            aCCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCCCGCTTGtgt                                                                                                                    <  3:8111‑M1/149‑4 (MQ=255)
                            aCCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCACAGCAGccc                                                                                                                                            >  5:390375/1‑125 (MQ=255)
                            aCCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCACAGCAGccc                                                                                                                                            <  6:390375/125‑1 (MQ=255)
                              cGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGC                                                                                                                  >  6:208929‑M1/1‑149 (MQ=255)
                              cGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGC                                                                                                                  >  2:227090‑M1/1‑149 (MQ=255)
                                gCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCG                                                                                                                 >  1:113981‑M1/1‑148 (MQ=255)
                                 cTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGTGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAA                                                                                                               <  3:185729‑M1/149‑1 (MQ=255)
                                                 gAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCgtggtg                                                                                                >  7:353304/1‑148 (MQ=255)
                                                  aaCCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGa                                                                                              <  5:315010/149‑1 (MQ=255)
                                                     cgccgTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATgg                                                                                           >  2:88560/1‑149 (MQ=255)
                                                     cgccgTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATgg                                                                                           >  4:380717/1‑149 (MQ=255)
                                                        cgTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCAc                                                                                        <  2:113981/149‑1 (MQ=255)
                                                          tgtACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACgg                                                                                      <  1:88560/149‑1 (MQ=255)
                                                            tACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGc                                                                                    <  8:367147/149‑1 (MQ=255)
                                                              cGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCt                                                                                   <  3:155596/148‑1 (MQ=255)
                                                                   ccGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGccc                                                                             <  2:103908/149‑1 (MQ=255)
                                                                     gTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCgt                                                                           <  8:48364/149‑1 (MQ=255)
                                                                         gCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTACAGGGCCATTGGGGCTGCTGGGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCGAAGCCCgtgcca                                                                       <  3:387473/149‑4 (MQ=38)
                                                                             ggtggtGTGAGAGGTCGGCGGGTGTGAATCCGGCTCCTAGTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACAGAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGttcgcgg                                                                   >  3:29795‑M2/1‑142 (MQ=255)
                                                                             ggtggtGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAcacgggcgaagcccgtgccatccaccacggtgttggcttcgcgg                                                                   <  6:4595‑M1/149‑45 (MQ=255)
                                                                                  tgtgAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGttcgcgggcaag                                                              >  6:105764‑M2/1‑137 (MQ=255)
                                                                                    tgAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGttcgcgggcaagcc                                                            <  3:306495‑M2/149‑15 (MQ=255)
                                                                                     gagaGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGG                                                                                                                           <  3:73727‑M1/85‑1 (MQ=255)
                                                                                     gagaGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGG                                                                                                                           >  4:73727‑M1/1‑85 (MQ=255)
                                                                                        gggaCGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGCCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTTGTGGATGGCACGGGCTTCCCCCCTGttcgcggtcaagcccgct                                                        <  3:244383‑M2/148‑19 (MQ=255)
                                                                                        agGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGttcgcgggcaagcccgct                                                        <  7:84897‑M2/149‑19 (MQ=255)
                                                                                        agGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAA                                                                                                               <  5:78900‑M1/94‑1 (MQ=255)
                                                                                        agGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAA                                                                                                               >  6:78900‑M1/1‑94 (MQ=255)
                                                                                         ggACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGttcgcgggcaa                                                               >  7:20540‑M2/1‑130 (MQ=255)
                                                                                         ggACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGttcgcgggcaa                                                               <  8:20540‑M2/141‑12 (MQ=255)
                                                                                         ggACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTc                                                       <  5:208929/149‑1 (MQ=255)
                                                                                         ggACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTc                                                       <  3:380717/149‑1 (MQ=255)
                                                                                           aCGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTccc                                                     <  4:95290/149‑1 (MQ=255)
                                                                                            cggcggGTGTGAAGCCGCCTCCTACTAAATTTTCAGGGCCATTGGGGCGGCTGTCCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTGCGCCCGTGTTCGCGGGCAAGCCCGCTCCca                                                    <  6:140355/149‑1 (MQ=25)
                                                                                                  gtgtGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTc                                              >  5:225062/1‑149 (MQ=255)
                                                                                                   tgtgAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTcg                                             >  7:371493/1‑149 (MQ=255)
                                                                                                    gtgAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAcacgggcgaagcccgtgccatccaccacggtgttggct                                                                         <  1:19881‑M1/120‑39 (MQ=255)
                                                                                                    gtgAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAcacgggcgaagcccgtgccatccaccacggtgttggct                                                                         >  2:19881‑M1/1‑82 (MQ=255)
                                                                                                      gAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGAGGGCAAGCCCGCTCCCACAGGTcgcg                                           >  4:149249/1‑148 (MQ=255)
                                                                                                            gCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTCCCCGCGAAGCCAACACCGTGGTGGATGCCACGGGCTACGCCCGTTTTCGCGGGCAACCCCGAGCCCCAAGGTCGCGAAAGGGc                                    >  3:210498/1‑149 (MQ=25)
                                                                                                              ctcctACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCTAAGCCAACACCGTGGTGGAGGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCg                                  <  3:381702/149‑1 (MQ=255)
                                                                                                                 ctACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCg                                  >  4:381702/1‑146 (MQ=255)
                                                                                                                  tACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCa                                >  7:308894/1‑147 (MQ=255)
                                                                                                                    cTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCagag                             <  8:353304/148‑1 (MQ=255)
                                                                                                                     tCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCg                           >  2:167256/1‑149 (MQ=255)
                                                                                                                     tCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCg                           >  2:232314/1‑149 (MQ=255)
                                                                                                                     tCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCg                           >  4:262805/1‑149 (MQ=255)
                                                                                                                     tCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCg                           >  5:395658/1‑149 (MQ=255)
                                                                                                                              aGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCTAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCa                  <  1:227090/149‑1 (MQ=255)
                                                                                                                                ggCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACAGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAATCCCGCTCCCACCGGTCGCGAAAGGGCCGCAGAGAGGCCCCGGAAAt                >  3:322196/1‑149 (MQ=25)
                                                                                                                                 gCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATc               >  1:120138/1‑149 (MQ=255)
                                                                                                                                        gggCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGttcgcgggcaag                                                              <  4:108882‑M2/95‑13 (MQ=255)
                                                                                                                                        gggCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGttcgcgggcaag                                                              >  3:108882‑M2/1‑83 (MQ=255)
                                                                                                                                          gctgctGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCGGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCt       >  1:322943/1‑148 (MQ=255)
                                                                                                                                           ctgctgTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGa     <  1:232314/149‑1 (MQ=255)
                                                                                                                                               tgtgCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCg  >  8:365153/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GGGCTGGTCTCGATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCG  >  NC_002947/2070998‑2071288

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: