Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,958,833 T→C 15.3% intergenic (+56/+86) PP_2589 → / ← PP_2590 aldehyde dehydrogenase family protein/ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,958,8330TC15.3% 167.5 / 20.3 72intergenic (+56/+86)PP_2589/PP_2590aldehyde dehydrogenase family protein/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (35/26);  new base C (5/6);  total (40/32)
Fisher's exact test for biased strand distribution p-value = 5.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

AAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGG  >  NC_002947/2958688‑2958980
                                                                                                                                                 |                                                                                                                                                   
aaGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGc                                                                                                                                                 >  2:1101922/1‑150 (MQ=255)
      tCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGt                                                                                                                                           >  1:1017151/1‑150 (MQ=255)
           gTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAg                                                                                                                                      <  1:687897/150‑1 (MQ=255)
                     tgacccagaccctggggttgggggcgccgcggacctatcagctttcggctagcgtggacttctaaccgggtcccgccgcggcgccagatacggcctcttcgcgggcatGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTg                                                                                                                            >  1:1446419‑M1/109‑150 (MQ=255)
                     caTCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTg                                                                                                                            <  1:388704/149‑1 (MQ=255)
                      aTCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAACCCTCAGACTTGTGCTga                                                                                                                           <  1:1074687/150‑2 (MQ=38)
                        cTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGAc                                                                                                                         <  1:1101922/150‑1 (MQ=255)
                          cagaccctggggttgggggcgccgcggacctatcagctttcggctagcgtggacttctaaccgggtcccgccgcggcgccagatacggcctcttcgcgggcatGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCt                                                                                                                       <  1:1065168‑M1/47‑1 (MQ=255)
                              cTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTgg                                                                                                                   <  1:103672/150‑1 (MQ=255)
                                           aCAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGGATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGcc                                                                                                      >  1:1559543/1‑150 (MQ=255)
                                             aaGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCcg                                                                                                    <  2:1816489/150‑1 (MQ=255)
                                                       aaCTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgt                                                                                                                     <  1:672997/123‑1 (MQ=255)
                                                       aaCTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgt                                                                                                                     >  2:672997/1‑123 (MQ=255)
                                                         cTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGt                                                                                        >  1:1261371/1‑150 (MQ=255)
                                                          tCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgt                                                                                                                     >  1:886462/1‑120 (MQ=255)
                                                          tCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgt                                                                                                                     <  2:886462/120‑1 (MQ=255)
                                                          tCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGgcgcgc                                                                                                       <  1:1309384/134‑1 (MQ=255)
                                                          tCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGgcgcgc                                                                                                       >  2:1309384/1‑134 (MQ=255)
                                                          tCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCAGCTCCAGATGGGTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTa                                                                                       >  1:1090148/1‑150 (MQ=255)
                                                           cAAGACCACCTGGTACCAGTTGCGCGTCTGATGTATTCATCATTGACCGCGCCGTCTGCTCCGCTGGCCCTCTCTCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTAt                                                                                      <  1:1790376/150‑1 (MQ=38)
                                                              gACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTg                                                                                   >  1:1761793/1‑150 (MQ=255)
                                                                ccaccTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGgc                                                                                 >  1:507026/1‑150 (MQ=255)
                                                                ccaccTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGgc                                                                                 <  1:1665978/150‑1 (MQ=255)
                                                                 caccTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGgcg                                                                                >  2:1766048/1‑150 (MQ=255)
                                                                 caccTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGgcg                                                                                >  1:624938/1‑150 (MQ=255)
                                                                    cTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgt                                                                                                                     >  1:1311313/1‑110 (MQ=255)
                                                                    cTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgt                                                                                                                     <  2:1311313/110‑1 (MQ=255)
                                                                         aCCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCgcggcg                                                                        >  2:38897/1‑150 (MQ=255)
                                                                                ttctaaccgggtcccgccgcggcgccagatacggcctcttcgcgggcgcGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGAc                                                                                                                         <  2:1852265‑M1/45‑1 (MQ=255)
                                                                                ttctaaccgggtcccgccgcggcgccagatacggcctcttcgcgggcgcGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGAc                                                                                                                         >  1:1852265‑M1/50‑94 (MQ=255)
                                                                                  gcgGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGt                                                               >  1:1316459/1‑150 (MQ=255)
                                                                                       tGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTgg                                                                                                                   <  2:1592363/93‑1 (MQ=14)
                                                                                       tGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCGCAGACTTGTGCTGAACCTGTgg                                                                                                                   >  1:1592363/1‑93 (MQ=17)
                                                                                           tcccgccgcggcgccagatacggcctcttcgcgggcgcGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCt                                                                                                                       >  1:417054‑M1/39‑85 (MQ=255)
                                                                                             aTTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCAc                                                    >  1:1178429/1‑150 (MQ=255)
                                                                                                 aTCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTa                                                <  2:507026/150‑1 (MQ=255)
                                                                                                 aTCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTa                                                >  1:898187/1‑150 (MQ=255)
                                                                                                    aTTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGcc                                             <  2:898187/150‑1 (MQ=255)
                                                                                                     ttGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCa                                                     >  2:1618786/1‑141 (MQ=255)
                                                                                                     ttGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCa                                                     <  1:1618786/141‑1 (MQ=255)
                                                                                                     ttGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCg                                            >  1:454967/1‑150 (MQ=255)
                                                                                                       gACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCGGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGaa                                          >  1:443939/1‑150 (MQ=255)
                                                                                                          cgcgTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTtgtg                                                                                                                               <  2:1645086/62‑1 (MQ=255)
                                                                                                          cgcgTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTtgtg                                                                                                                               >  1:1645086/1‑62 (MQ=255)
                                                                                                                 tcTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGCGGGTTGTGCTTGACCTGTGGGAGCGTGCGCGCCCGCGAAGAGGCCGTATCTGCCGCCGCGGCGGGACCCGGTTGGAAGTCCACGCTAGCCGAAAGCTGATAgg                                >  2:84759/1‑150 (MQ=39)
                                                                                                                 tcTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAgg                                <  1:38897/150‑1 (MQ=255)
                                                                                                                    gCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTcc                             >  1:984286/1‑150 (MQ=255)
                                                                                                                         gcgGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCgcggc                        <  2:454967/150‑1 (MQ=255)
                                                                                                                             gCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGccc                    >  1:408667/1‑150 (MQ=255)
                                                                                                                               aTGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTGGCCGAAAGCTGATAGGTCCGCGGCGccccc                  >  2:644358/1‑150 (MQ=255)
                                                                                                                                     GCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAAcccc            <  2:1065168‑M1/150‑1 (MQ=255)
                                                                                                                                       TCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGAt                                   <  2:1113636‑M1/125‑1 (MQ=255)
                                                                                                                                       TCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGAt                                   >  1:1113636‑M1/1‑125 (MQ=255)
                                                                                                                                        CCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAgg         >  1:1766048‑M1/1‑150 (MQ=255)
                                                                                                                                         CCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGa                                    >  1:1867496‑M1/1‑122 (MQ=255)
                                                                                                                                         CCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGa                                    <  2:1867496‑M1/122‑1 (MQ=255)
                                                                                                                                           ACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCgg                                                                                                             >  2:1914998‑M1/1‑47 (MQ=255)
                                                                                                                                           ACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCgg                                                                                                             <  1:1914998‑M1/47‑1 (MQ=255)
                                                                                                                                           ACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGTGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGaa                                                          <  1:747068‑M1/98‑1 (MQ=255)
                                                                                                                                           ACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGTGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGaa                                                          >  2:747068‑M1/1‑98 (MQ=255)
                                                                                                                                           ACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTa                                                >  2:1175102‑M1/1‑108 (MQ=255)
                                                                                                                                           ACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTa                                                <  1:1175102‑M1/108‑1 (MQ=255)
                                                                                                                                           ACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATa                                  <  1:22480‑M1/122‑1 (MQ=255)
                                                                                                                                           ACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATa                                  >  2:22480‑M1/1‑122 (MQ=255)
                                                                                                                                           ACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGACCGCGAAGAGGCCGTATATGGCGACGCGGCGGGACCAGGTTAGAAGTTCACTGGAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTc      <  2:707798‑M1/150‑1 (MQ=255)
                                                                                                                                           ACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGAt                                   <  1:349053‑M1/121‑1 (MQ=255)
                                                                                                                                           ACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGAt                                   >  2:349053‑M1/1‑121 (MQ=255)
                                                                                                                                           ACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTc      <  2:1850225‑M1/150‑1 (MQ=255)
                                                                                                                                             AGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGa                                                                     >  2:417054‑M1/1‑85 (MQ=255)
                                                                                                                                              GGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTgg   >  2:1178429‑M1/1‑150 (MQ=255)
                                                                                                                                              GGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTgg   >  2:984286‑M1/1‑150 (MQ=255)
                                                                                                                                               gTTCAGCACTCGTCTGAGGTTTGTGCTGGACCTGTGGGAGGGGGCGCGCGCGCGAAGAGGCCGTATCTGCCGCCGCGGCGTGTCCCGGTTAGACATCCACGCTACCCGAAAGCGGTTAGTGCCGCGGCGCCCCCCACCCCAGGGTCTggg  <  2:615677/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
AAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGG  >  NC_002947/2958688‑2958980

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