Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,444,978 A→G 5.3% intergenic (‑151/‑115) PP_4781 ← / → thiD sensor histidine kinase/phosphomethylpyrimidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,444,9780AG 5.3% 224.7 / 2.2 75intergenic (‑151/‑115)PP_4781/thiDsensor histidine kinase/phosphomethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base A (37/34);  new base G (2/2);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

GCCTGTCCGGTCGCGTCGCGGTAGCCTCCCATCCGTTTGAGACGACGTGATCCTGCAGAACCCGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAGGTTTGGTGGTAAGCTCGCCGACCATGAATACCTACAGCTCCCGCCCCGTT  >  NC_002947/5444842‑5445119
                                                                                                                                        |                                                                                                                                              
gCCTGTCCGGTCGCGTCGCGGTAGCCTCCCATCCGTTTGAGACGACGTGATCCTGCAGAACCCGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCaa                                                                                                                                   >  1:441185/1‑150 (MQ=255)
gCCTGTCCGGTCGCGTCGCGGTAGCCTCCCATCCGTTTGAGACGACGTGATCCTGCAGAACCCGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCaa                                                                                                                                   <  1:354709/150‑1 (MQ=255)
     tCCGGTCGCGTCGCGGTAGCCTCCCATCCGTTTGAGACGACGTGATCCTGCAGAACCCGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGaaaa                                                                                                                                                <  1:511893/132‑1 (MQ=255)
     tCCGGTCGCGTCGCGGTAGCCTCCCATCCGTTTGAGACGACGTGATCCTGCAGAACCCGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGaaaa                                                                                                                                                >  2:511893/1‑132 (MQ=255)
             cgtcgcgGTAGCCTCCCATCCGTTTGAGACGACGTGATCCTGCAGAACCCGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAGGGCTCTAAATCAAt                                                                                                                                  <  1:690166/138‑1 (MQ=255)
             cgtcgcgGTAGCCTCCCATCCGTTTGAGACGACGTGATCCTGCAGAACCCGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAGGGCTCTAAATCAAt                                                                                                                                  >  2:690166/1‑138 (MQ=255)
                 gcgGTAGCCTCCCATCCGTTTGAGACGACGTGATCCTGCAGAACCCTTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAAc                                                                                                                  <  1:258697/150‑1 (MQ=255)
                        ccTCCCATCCGTTTGAGACGACGTGATCCTGCAGAACCCGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCttt                                                                                                                             <  1:466797/132‑1 (MQ=255)
                        ccTCCCATCCGTTTGAGACGACGTGATCCTGCAGAACCCGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCttt                                                                                                                             >  2:466797/1‑132 (MQ=255)
                            ccATCCGTTTGAGACGACGTGATCCTGCAGAACCCGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAAc                                                                                                       >  2:1079195/1‑150 (MQ=255)
                             cATCCGTTTGAGACGACGTGATCCTGCANAACCCGTTNCTCGAGTTANCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAAcc                                                                                                      >  1:329701/1‑150 (MQ=255)
                             cATCCGTTTGAGACGACGTGATCCTGCAGAACCCGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAAcc                                                                                                      >  2:546812/1‑150 (MQ=255)
                              aTCCGTTTGAGACGACGTGATCCTGCAGAACCCGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATa                                                                                                                        >  2:361581/1‑131 (MQ=255)
                              aTCCGTTTGAGACGACGTGATCCTGCAGAACCCGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATa                                                                                                                        <  1:361581/131‑1 (MQ=255)
                                 cGTTTGAGACGACGTGATCCTGCAGAACCCGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGAc                                                                                                  >  1:333452/1‑150 (MQ=255)
                                 cGTTTGAGACGACGTGATCCTGCAGAACCCGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGAc                                                                                                  <  1:125033/150‑1 (MQ=255)
                                  gTTTGAGACGACGTGATCCTGCAGAACCCGTTTATCGAGTTATGCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGAcc                                                                                                 >  2:873966/1‑150 (MQ=255)
                                   tttGAGACGACGTGATCCTGCAGAACCCGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACcg                                                                                                >  2:795477/1‑150 (MQ=255)
                                     tGAGACGACGTGATCCTGCAGAACCCGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACcgat                                                                                             <  1:795477/150‑1 (MQ=255)
                                              gTGATCCTGCAGAACCCGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCaa                                                                                                                          <  1:935369/113‑1 (MQ=255)
                                              gTGATCCTGCAGAACCCGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCaa                                                                                                                          >  2:935369/1‑113 (MQ=255)
                                                 aTCCTGCAGAACCCGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTaaaa                                                                                 <  1:127689/150‑1 (MQ=255)
                                                  tCCTGCAGAACCCGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTaaaa                                                                                 >  2:1040330/1‑149 (MQ=255)
                                                  tCCTGCAGAACCCGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTaaaa                                                                                 <  1:1040330/149‑1 (MQ=255)
                                                    cTGCAGAACCCGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTc                                                                              >  1:1004336/1‑150 (MQ=255)
                                                           aaccGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCACTCAATAACCCGACCGATCGATCAGGAAAGCTCGTTTGTc                                                                       <  2:707723/148‑1 (MQ=255)
                                                             ccGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCaa                                                                                                                                   >  1:99097/1‑89 (MQ=255)
                                                             ccGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCaa                                                                                                                                   <  2:99097/89‑1 (MQ=255)
                                                             ccGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCgg                                                                     >  1:324791/1‑150 (MQ=255)
                                                               gTTTATGGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGAc                                                                   >  2:203563/1‑150 (MQ=255)
                                                               gTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCttt                                                                                                                             >  2:1193968/1‑93 (MQ=255)
                                                               gTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCttt                                                                                                                             <  1:1193968/93‑1 (MQ=255)
                                                                  tATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGAt                                                                >  1:114094/1‑150 (MQ=255)
                                                                    tCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATcc                                                              >  2:212073/1‑150 (MQ=255)
                                                                       aGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGcc                                                           <  2:1004336/150‑1 (MQ=255)
                                                                              cGAGCACTCTAGCATAGCCTCACAGCCGGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTAACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAg                                                    <  2:90228/150‑1 (MQ=255)
                                                                               gAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATTCTTTCAATAGaaaa                                                                                                                   <  2:1199925/87‑1 (MQ=255)
                                                                               gAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGaaaa                                                                                                                   >  1:1199925/1‑87 (MQ=255)
                                                                               gAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGAc                                                                                                  >  1:335202/1‑104 (MQ=255)
                                                                               gAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGAc                                                                                                  <  2:335202/104‑1 (MQ=255)
                                                                               gAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAgg                                                   <  2:333452/150‑1 (MQ=255)
                                                                               gAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAgg                                                   >  2:1013289/1‑150 (MQ=255)
                                                                                 gCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAGGtt                                                 <  1:1079195/150‑1 (MQ=255)
                                                                                     tctAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAACCTCGTTTGTCGGACGATCCGCCCAAACAGGTTtggt                                             >  2:400409/1‑150 (MQ=255)
                                                                                     tctAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAGGTTtggt                                             <  2:441185/150‑1 (MQ=255)
                                                                                         gCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAGGTTTGGTGGTa                                         >  1:607411/1‑150 (MQ=255)
                                                                                            tAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAg                                                    <  2:538048/136‑1 (MQ=255)
                                                                                            tAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAg                                                    >  1:538048/1‑136 (MQ=255)
                                                                                             aGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATc                                                               <  1:1013475/124‑1 (MQ=255)
                                                                                             aGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATc                                                               >  2:1013475/1‑124 (MQ=255)
                                                                                                 tcaaagccTGCCAATCGGCAACATAAAAACACAAAGAAAAGGCTCTAAATCAATTCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAACCTCGTTTGTCGGACGATCCGCCCAAACAGGTTTGGTGGTAAGCTCGcc                                 <  2:329701/144‑1 (MQ=255)
                                                                                                 tCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAGGGCTCTAAATCAATGCTTTCAATAGaa                                                                                                                     >  1:844567/1‑67 (MQ=255)
                                                                                                 tCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAGGGCTCTAAATCAATGCTTTCAATAGaa                                                                                                                     <  2:844567/67‑1 (MQ=255)
                                                                                                 tCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGccc                                                          >  2:1143343/1‑125 (MQ=255)
                                                                                                 tCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGccc                                                          <  1:1143343/125‑1 (MQ=255)
                                                                                                 tCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAGGTTtggt                                             <  2:1148846/138‑1 (MQ=255)
                                                                                                 tCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAGGTTtggt                                             >  1:1148846/1‑138 (MQ=255)
                                                                                                 tCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAGGTTTGGTGGTAAGCTCGcc                                 >  2:643271/1‑150 (MQ=255)
                                                                                                  cacaACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGt                                                                            <  2:438122/106‑1 (MQ=255)
                                                                                                  cacaACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGt                                                                            >  1:438122/1‑106 (MQ=255)
                                                                                                   acaACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGTTCGAAATCAATGCGTTCAATAGAACACCAATCACTCAAGTGACTGCATCGGTAGTATGACTGGTTTGTCGGACGATCCGCCCAAACAGGTTTGGTGGTATGCTCGCCGa                               >  1:1177878/1‑150 (MQ=255)
                                                                                                     aaCCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAGGTTTGGTGGTAAGCTCGCCGAcc                             <  1:873966/150‑1 (MQ=255)
                                                                                                        cTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAg                                                    >  1:957935/1‑124 (MQ=255)
                                                                                                        cTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAg                                                    <  2:957935/124‑1 (MQ=255)
                                                                                                        cTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTAAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAGGTTTGGTGGTAAGCTCGCCGACCATg                          >  2:234184/1‑150 (MQ=255)
                                                                                                             aaTCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAGGTTTGGTGGTAAGCTCGCCGACCATGAATAc                     <  1:1170737/150‑1 (MQ=255)
                                                                                                               tCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAGGTTTGGTGGTAAGCTCGCCGACCATGAATACCt                   >  2:204479/1‑150 (MQ=255)
                                                                                                                cGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAGGTTTGGTGGTAAGCTCGCCGACCATGAATACCTa                  <  1:204479/150‑1 (MQ=255)
                                                                                                                 ggCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCATTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAGGTTTGGTGGTAAGCTCGCCGACCATGAATACCTAc                 >  1:1147928/1‑150 (MQ=255)
                                                                                                                  gCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAGGTTTGGTGGTAAGCTCGCCGACCATGAATACCTACa                >  2:677200/1‑150 (MQ=255)
                                                                                                                  gCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAGGTTTGGTGGTAAGCTCGCCGACCATGAATACCTACa                <  2:324791/150‑1 (MQ=255)
                                                                                                                      caTAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAGGTTTGGTGGTAAGCTCGCCGACCATGAATACCTACAGCTc            <  1:212073/150‑1 (MQ=255)
                                                                                                                       aTAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAGGTTTGGTGGTAAGCTCGCCGACCATGAATACCTACAGCTcc           >  1:44097/1‑150 (MQ=255)
                                                                                                                             acacAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAGGTTTGGTGGTAAGCTCGCCGACCATGAATACCTACAGCTCCCGcccc     <  1:546812/150‑1 (MQ=255)
                                                                                                                                caAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAGGTTTGGTGGTAAGCTCGCCGACCATGAATACCTACAGCTCCCGCCCCGtt  <  2:607411/150‑1 (MQ=255)
                                                                                                                                 aaaGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAGGTTTGGTGGTAAGCTCGCCGACCAt                           <  2:515579/124‑1 (MQ=255)
                                                                                                                                 aaaGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAGGTTTGGTGGTAAGCTCGCCGACCAt                           >  1:515579/1‑124 (MQ=255)
                                                                                                                                     aaaaGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAGGTTTGGTGGTAAGCTCGCCGACCATGAATACCTACAGCTCCCGCCCCGtt  <  2:584230/145‑1 (MQ=255)
                                                                                                                                     aaaaGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAGGTTTGGTGGTAAGCTCGCCGACCATGAATACCTACAGCTCCCGCCCCGtt  >  1:584230/1‑145 (MQ=255)
                                                                                                                                        |                                                                                                                                              
GCCTGTCCGGTCGCGTCGCGGTAGCCTCCCATCCGTTTGAGACGACGTGATCCTGCAGAACCCGTTTATCGAGTTATCCGAGCACTCTAGCATAGCCTCACAACCTGGCAATCGGCCACATAAAAACACAAAGAAAAGGCTCTAAATCAATGCTTTCAATAGAAAACCAATCAATAACCTGACCGATCGATCAGTAAAACTCGTTTGTCGGACGATCCGCCCAAACAGGTTTGGTGGTAAGCTCGCCGACCATGAATACCTACAGCTCCCGCCCCGTT  >  NC_002947/5444842‑5445119

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: