Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,081,044 T→A 12.6% intergenic (+35/+27) PP_5332 → / ← PP_5333 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,081,0440TA12.6% 193.0 / 19.1 72intergenic (+35/+27)PP_5332/PP_5333hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/35);  new base A (5/4);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.01e-01

TGGATGGCAGCGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCGGGGTGTAGAT  >  NC_002947/6080895‑6081193
                                                                                                                                                     |                                                                                                                                                     
tggatggCAGCGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCt                                                                                                                                                       <  2:286041/150‑1 (MQ=255)
 ggatggCAGCGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCtt                                                                                                                                                      <  1:1038871/150‑1 (MQ=255)
  gatggCAGCGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTg                                                                                                                                                     >  1:347410/1‑150 (MQ=255)
          cGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCAAGCCCGCGaa                                                                                                                                             <  1:1193102/150‑1 (MQ=255)
            aaTGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAaga                                                                                                                                           <  2:1936902/150‑1 (MQ=255)
            aaTGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAaga                                                                                                                                           <  1:2148412/150‑1 (MQ=255)
                 tgtgCATCGAGACAGCGAATGTTTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCg                                                                                                                                      <  2:2379211/150‑1 (MQ=255)
                 tgtgCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCg                                                                                                                                      >  2:1824782/1‑150 (MQ=255)
                   tgCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGCCTTTCCCGCGAAGAGGCCggg                                                                                                                                    <  2:2041031/150‑1 (MQ=255)
                    gCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCAAGCCCGCGa                                                                                                                                              <  1:621967/139‑1 (MQ=255)
                    gCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCAAGCCCGCGa                                                                                                                                              >  2:621967/1‑139 (MQ=255)
                              aGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCaa                                                                                                                         >  1:2158422/1‑150 (MQ=255)
                                 gAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAgg                                                                                                                      >  1:1478061/1‑150 (MQ=255)
                                  aaTGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGGGTAGAGGCCGGGGCAATTATCAAAGGt                                                                                                                     >  1:1074/1‑150 (MQ=255)
                                      tgtgGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGc                                                                                                                 <  2:1088350/150‑1 (MQ=255)
                                      tgtgGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGc                                                                                                                 <  1:2138091/150‑1 (MQ=255)
                                        tgGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGt                                                                                                               <  1:2109486/150‑1 (MQ=255)
                                                      ggAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTATCGGACGCTTCGCGGGCGTGCCCGCGATGAGGCAGGAGATATTATCAACGCACTGCGACTTCCTCCACCCTc                                                                                                 >  1:1223633/1‑150 (MQ=255)
                                                          ttGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCt                                                                                             >  1:453527/1‑150 (MQ=255)
                                                              aGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTc                                                                                         <  2:347410/150‑1 (MQ=255)
                                                                ccTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTg                                                                                       >  1:2346254/1‑150 (MQ=255)
                                                                 cTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGc                                                                                      <  1:1462800/150‑1 (MQ=255)
                                                                  tGGGGCGAGTCATTCGCGGGGCGTCCAGATCGGTTGTGGAAGTATCTGAGCCCTGCTCTTGTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCa                                                                                     <  2:2361899/150‑1 (MQ=255)
                                                                             aTTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCaccacc                                                                                                    >  2:566835/1‑124 (MQ=255)
                                                                             aTTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCaccacc                                                                                                    <  1:566835/124‑1 (MQ=255)
                                                                              ttCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAg                                                                                >  2:800475/1‑143 (MQ=255)
                                                                              ttCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAg                                                                                <  1:800475/143‑1 (MQ=255)
                                                                                 gCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGcc                                                                      >  1:2444328/1‑150 (MQ=255)
                                                                                         tCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGc                                                              <  2:2346254/150‑1 (MQ=255)
                                                                                         tCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGc                                                              >  2:701580/1‑150 (MQ=255)
                                                                                           gaGATCGGTTGTGAAATTATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCa                                                            <  2:1223633/149‑1 (MQ=255)
                                                                                           cAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCa                                                            <  2:1630731/150‑1 (MQ=255)
                                                                                           cAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCa                                                            >  1:1665514/1‑150 (MQ=255)
                                                                                            aGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAagag                                                                                                                                          <  1:253197/71‑1 (MQ=255)
                                                                                            aGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAagag                                                                                                                                          >  2:253197/1‑71 (MQ=255)
                                                                                            aGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAg                                                           <  1:1274840/150‑1 (MQ=255)
                                                                                            aGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAg                                                           <  2:317226/150‑1 (MQ=255)
                                                                                                      tgAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCCGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTg                                                 >  1:2200904/1‑150 (MQ=255)
                                                                                                        aaaaTATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGTGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATTCAGCGCATAGGCCGCATTGGCCAGGTCATTGCTgcc                                               >  2:1004499/1‑148 (MQ=255)
                                                                                                        aaaaTATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGAc                                               <  2:1478061/150‑1 (MQ=255)
                                                                                                         aaaTATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATGGGCCGCATCGGCCAGGTCATTGCTGAcc                                              <  2:2152032/150‑1 (MQ=255)
                                                                                                         aaaTATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGAcc                                              <  2:1665514/150‑1 (MQ=255)
                                                                                                                tGAGCCCTGTTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCa                                       <  2:1074/150‑1 (MQ=255)
                                                                                                                tGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCa                                       >  2:1476712/1‑150 (MQ=255)
                                                                                                                      cTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCa                                 <  1:1824782/150‑1 (MQ=255)
                                                                                                                      cTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCa                                 >  2:174862/1‑150 (MQ=255)
                                                                                                                           ctTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCATCTGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTg                            <  1:1208231/150‑1 (MQ=255)
                                                                                                                            tttttttACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGt                                                                                          >  1:264874/1‑87 (MQ=255)
                                                                                                                            tttttttACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGt                                                                                          <  2:264874/87‑1 (MQ=255)
                                                                                                                            tttttttACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGc                           <  1:1476712/150‑1 (MQ=255)
                                                                                                                            tttttttACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGc                           <  1:174862/150‑1 (MQ=255)
                                                                                                                                         ctctTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTtc                                           <  1:427445/121‑1 (MQ=255)
                                                                                                                                         ctctTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTtc                                           >  2:427445/1‑121 (MQ=255)
                                                                                                                                          tctTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCTTAAACATAGGCCGCATCGGCCAGGGCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGc             >  1:643834/1‑150 (MQ=255)
                                                                                                                                          tctTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGc             >  2:1441895/1‑150 (MQ=255)
                                                                                                                                          tctTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGc             >  1:1170359/1‑150 (MQ=255)
                                                                                                                                          tctTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGc             >  2:870636/1‑150 (MQ=255)
                                                                                                                                          tctTCGCGGGCAAGCCCGCGAAGAGGCCGGAGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGc             >  1:1217445/1‑150 (MQ=255)
                                                                                                                                           ctTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGAc                                               <  2:686128/115‑1 (MQ=255)
                                                                                                                                           ctTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGAc                                               >  1:686128/1‑115 (MQ=255)
                                                                                                                                            ttCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCgg           >  1:1889107/1‑150 (MQ=255)
                                                                                                                                            ttCGCGGGCAAGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCAcc                                                                                              <  2:915803/67‑1 (MQ=37)
                                                                                                                                            ttCGCGGGCAAGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCAcc                                                                                              >  1:915803/1‑67 (MQ=37)
                                                                                                                                            ttCGCGGGCAAGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCa                                                       >  1:613374/1‑106 (MQ=255)
                                                                                                                                            ttCGCGGGCAAGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCa                                                       <  2:613374/106‑1 (MQ=255)
                                                                                                                                             tCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCggg          <  1:870636/150‑1 (MQ=255)
                                                                                                                                             tCGCGGGCAAGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCggg          >  2:2188302/1‑150 (MQ=255)
                                                                                                                                              cgcgGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCgggg         <  1:1441895/150‑1 (MQ=255)
                                                                                                                                               gcgGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCGGGgt        <  1:701580/150‑1 (MQ=255)
                                                                                                                                               gcgGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCGGGgt        <  2:453527/150‑1 (MQ=255)
                                                                                                                                                cgGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCgg                                                               <  1:1101159/94‑1 (MQ=255)
                                                                                                                                                cgGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCgg                                                               <  1:1102499/94‑1 (MQ=255)
                                                                                                                                                cgGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCgg                                                               >  2:1101159/1‑94 (MQ=255)
                                                                                                                                                cgGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCgg                                                               >  2:1102499/1‑94 (MQ=255)
                                                                                                                                                    tttGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCGGGGTGTAGa   <  2:1856110/149‑1 (MQ=255)
                                                                                                                                                     ttGNCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCGGGGTGTAGat  >  2:1506930/1‑150 (MQ=255)
                                                                                                                                                     ttGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCGGGGTGTAGat  >  2:1352993/1‑150 (MQ=255)
                                                                                                                                                     ttGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCGGGGTGTAGat  <  2:1889107/150‑1 (MQ=255)
                                                                                                                                                     ttGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGATTGACGCTGACCCACAGCGGGGCGTAGat  >  1:1697577/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGATGGCAGCGAATGCTGTGCATCGAGACAGCGAATGTGTGGGATGACCTGGTGGAGTTGAAGCCTGGGGCGAGTCATTCGCTGGGCGTCCAGATCGGTTGTGAAAATATCTGAGCCCTGCTCTTTTTTTACCGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGGGCAATTATCAAAGGTCTGCGTCTTCCACCACCCTCACCTTGTCTGCATCCAGCGCATAGGCCGCATCGGCCAGGTCATTGCTGACCTTCTCCACCTTCAGCTTGCCGCTGACCCACAGCGGGGTGTAGAT  >  NC_002947/6080895‑6081193

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