Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,968,770	G→A	R199C (CGT→TGT)	cheR ←	chemotaxis regulator, protein‑glutamate methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,968,770	0	G	A	98.9%	297.4 / ‑4.4	87	R199C (CGT→TGT)	cheR	chemotaxis regulator, protein‑glutamate methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (46/40);  total (46/41)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

TTACGACAGAAGATCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGAT  >  NC_000913/1968673‑1968866                                                                                                  |                                                                                                 ttACGACAGAAGATCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGacaca                                                                                                <  2:2015703/100‑1 (MQ=255)  tACGACAGAAGATCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGacacac                                                                                               >  2:1273841/1‑100 (MQ=255)    cGACAGAAGATCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGacacacgc                                                                                             >  2:2040650/1‑100 (MQ=255)    cGACAGAAGATCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGacacacgc                                                                                             >  1:804535/1‑100 (MQ=255)     gACAGAAGATCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACAcgcg                                                                                            >  2:333260/1‑100 (MQ=255)       cAGAAGATCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTa                                                                                          <  2:19736/100‑1 (MQ=255)        agaagaTCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTAc                                                                                         >  1:664927/1‑100 (MQ=255)         gaagaTCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTAcc                                                                                        >  1:792916/1‑100 (MQ=255)         gaagaTCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTAcc                                                                                        <  1:439585/100‑1 (MQ=255)           agaTCGCATCAAACGGCCCTGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAg                                                                                      >  2:2462789/1‑100 (MQ=255)              tCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGccc                                                                                   >  1:424343/1‑100 (MQ=255)               cGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCt                                                                                  >  2:2533475/1‑100 (MQ=255)                gCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCtt                                                                                 <  1:2321169/100‑1 (MQ=255)                  aTCAAACGGCCCCGGCACGGTGTCCTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCa                                                                               <  2:980890/100‑1 (MQ=255)                  aTCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCa                                                                               >  2:2171277/1‑100 (MQ=255)                     aaaCGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGc                                                                            >  1:618646/1‑100 (MQ=255)                      aaCGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCg                                                                           <  1:2549303/100‑1 (MQ=255)                      aaCGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCg                                                                           <  1:1487782/100‑1 (MQ=255)                       aCGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCgg                                                                          <  1:1863896/100‑1 (MQ=255)                        cGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGc                                                                         >  2:1886081/1‑100 (MQ=255)                        cGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGc                                                                         <  1:2375623/100‑1 (MQ=255)                         ggCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGcc                                                                        <  2:20676/100‑1 (MQ=255)                         ggCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGcc                                                                        <  1:592808/100‑1 (MQ=255)                          gCCCCGGTACGGTGTACAGTTTCGCCAGTAGATTCAGCGGGGCATAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGccc                                                                       <  1:2428656/100‑1 (MQ=255)                           ccccGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGcccc                                                                      >  2:2226534/1‑100 (MQ=255)                           ccccGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGcccc                                                                      >  2:16698/1‑100 (MQ=255)                           ccccGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGcccc                                                                      >  1:2059714/1‑100 (MQ=255)                            cccGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCg                                                                     <  2:1473783/100‑1 (MQ=255)                            cccGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCg                                                                     >  1:1867011/1‑100 (MQ=255)                            cccGGCACGGTGTACTGTTACGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCg                                                                     <  1:873664/100‑1 (MQ=255)                                  aCGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTc                                                               >  1:1016097/1‑100 (MQ=255)                                    ggTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGc                                                             <  1:971096/100‑1 (MQ=255)                                    ggTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGc                                                             >  1:1784953/1‑100 (MQ=255)                                    ggTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGc                                                             <  1:1614112/100‑1 (MQ=255)                                      tgtACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCAt                                                           <  1:642704/100‑1 (MQ=255)                                       gtACTGTTTCGCCCGTAGATTCAGCGGACCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATg                                                          >  1:2324639/1‑100 (MQ=255)                                       gtACTGTTTCGCCAGTAGCTTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATCCGGCCCCGTCCCCCCCATg                                                          >  2:2314372/1‑100 (MQ=255)                                       gtACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATg                                                          >  2:209390/1‑100 (MQ=255)                                        tACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGa                                                         <  2:1928093/100‑1 (MQ=255)                                          cTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGaaa                                                       >  2:174311/1‑100 (MQ=255)                                           tGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAAt                                                      >  1:1830380/1‑100 (MQ=255)                                                 gCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGtt                                                >  1:1797684/1‑100 (MQ=255)                                                     gTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAg                                            <  1:2623388/100‑1 (MQ=255)                                                     gTAGATTCAGCGGGACAAAATGAACATAGTTCGCCAGCTCCTGACACACGCGTAACAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAg                                            >  2:697960/1‑100 (MQ=255)                                                        gATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTtg                                         >  2:2273751/1‑100 (MQ=255)                                                        gATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTtg                                         >  2:1774514/1‑100 (MQ=255)                                                             aGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTgcg                                    <  1:531499/100‑1 (MQ=255)                                                             aGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTgcg                                    >  1:891961/1‑100 (MQ=255)                                                             aGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTgcg                                    <  2:790173/100‑1 (MQ=255)                                                              gCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTgcgg                                   >  1:1513898/1‑100 (MQ=255)                                                               cGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTgcggc                                  <  2:616126/100‑1 (MQ=255)                                                                ggggCTAAATCAACATAGTTCGCCAACCCCTGACACACGCGTACCAGCACTTCATGCGTCCCCCTCACTCGCTTGAAATACCGTTGCAGTTGCTgcggcg                                 >  2:106809/1‑100 (MQ=255)                                                                ggggCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTgcggcg                                 <  1:2183668/100‑1 (MQ=255)                                                                ggggCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTgcggcg                                 <  1:1401063/100‑1 (MQ=255)                                                                  ggCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTc                               >  2:1945678/1‑100 (MQ=255)                                                                  ggCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTc                               >  2:694010/1‑100 (MQ=255)                                                                  ggCAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTc                               >  1:1288215/1‑100 (MQ=255)                                                                    cAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAg                             >  1:1497794/1‑100 (MQ=255)                                                                    cAAAATCAACATAGTTCGCCAGCTCCTGACACACGCGCACCAGCCCCTCATGCGGCCCCGTCCCTTGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAg                             <  1:2456006/100‑1 (MQ=255)                                                                     aaaaTCAACATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAgg                            >  1:600495/1‑100 (MQ=255)                                                                            aCATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCa                     <  1:839204/100‑1 (MQ=255)                                                                            aCATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCa                     <  2:686971/100‑1 (MQ=255)                                                                             cATAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCaa                    >  2:1240792/1‑100 (MQ=255)                                                                              aTAGTTCGCCAGCTCCTGACACACTCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAAc                   <  1:1202680/100‑1 (MQ=255)                                                                              aTAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCaan                   <  1:2097850/100‑2 (MQ=255)                                                                              aTAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAAc                   >  1:1971165/1‑100 (MQ=255)                                                                              aTAGTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAAc                   <  2:1170163/100‑1 (MQ=255)                                                                               tAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAAct                  <  1:390915/100‑1 (MQ=255)                                                                               tAGTTCGCCAGCTCCTGACACACGCGAACCAGCCCTTCATGCGGCCCCGCCCCTCGCATGAACTACCGTTGCAGTTGCTCCGGCGTCAGGTTTTTCAAct                  >  2:309923/1‑100 (MQ=255)                                                                                 gTTCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAActct                >  2:668086/1‑100 (MQ=255)                                                                                  ttCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCtt               <  1:1925059/100‑1 (MQ=255)                                                                                   tCGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTc              >  1:1582488/1‑100 (MQ=255)                                                                                    cGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCa             >  2:2595292/1‑100 (MQ=255)                                                                                    cGCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCa             >  1:1318716/1‑100 (MQ=255)                                                                                     gCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCAt            >  1:372290/1‑100 (MQ=255)                                                                                     gCCAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCAt            <  1:1661017/100‑1 (MQ=255)                                                                                      ccAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATg           >  2:306520/1‑100 (MQ=255)                                                                                      ccAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATg           >  1:664549/1‑100 (MQ=255)                                                                                      ccAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATg           >  2:1992479/1‑100 (MQ=255)                                                                                      ccAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTACTGCGGCGTCAGGTTTTTCAACTCTTCATg           >  2:1132070/1‑100 (MQ=255)                                                                                       cAGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTCCTGCGGCGTCAGGTTTTTCAACTCTTCATgg          <  2:1069847/100‑1 (MQ=255)                                                                                        aGCTCCTGACACACGCGTTCCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGc         <  1:2478428/100‑1 (MQ=255)                                                                                        aGCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGc         <  1:2569038/100‑1 (MQ=255)                                                                                         gCTCCTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATCCCGTTGCAGTTGCTGCGGCGTCAGGGCTTTCAACTCTTCATGGCg        >  2:1602474/1‑100 (MQ=255)                                                                                            ccTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCgata     <  1:2013320/100‑1 (MQ=255)                                                                                            ccTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCgata     <  1:1040510/100‑1 (MQ=255)                                                                                             cTGACACGCGGGTACCAGACCTTCAGGCGGCCCTGTCCCTCGCATGAAATCCCGTTGCAGTTGCTGAGGCGTCAGGTTTTTCCACTCTTCATGGCgatag    <  2:247802/100‑1 (MQ=255)                                                                                             cTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCgatag    <  1:1878268/100‑1 (MQ=255)                                                                                             cTGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCgatag    <  1:1765963/100‑1 (MQ=255)                                                                                              tGACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCgataga   <  1:1701574/100‑1 (MQ=255)                                                                                               gACACACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCgatagat  >  1:1264500/1‑100 (MQ=255)                                                                                                  |                                                                                                 TTACGACAGAAGATCGCATCAAACGGCCCCGGCACGGTGTACTGTTTCGCCAGTAGATTCAGCGGGGCAAAATCAACATAGTTCGCCAGCTCCTGACGCACGCGTACCAGCCCTTCATGCGGCCCCGTCCCTCGCATGAAATACCGTTGCAGTTGCTGCGGCGTCAGGTTTTTCAACTCTTCATGGCGATAGAT  >  NC_000913/1968673‑1968866

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: ‑