Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	243,573	A→C	N11H (AAC→CAC)	gmhA →	D‑sedoheptulose 7‑phosphate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	243,573	0	A	C	100.0%	234.7 / NA	64	N11H (AAC→CAC)	gmhA	D‑sedoheptulose 7‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (38/26);  total (38/26)

TTTATCCGCACTTCAGGTCAAAAAGTCCTGGTCATAGCACCTGCCCGTACTTCTCGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATT  >  NC_000913/243426‑243711                                                                                                                                                    |                                                                                                                                           tttATCCGCACTTCAGGTCAAAAAGTCCTGGTCATAGCACCTGCCCGTACTTCTCGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCAc                                                                                                                                          <  2:614272/150‑1 (MQ=255)   tATCCGCACTTCAGGTCAAAAAGTCCTGGTCATAGCACCTGCCCGTACTTCTCGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGa                                                                                                                                        <  2:778901/150‑1 (MQ=255)    aTCCGCACTTCAGGTCAAAAAGTCCTGGTCATAGCACCTGCCCGTCCTTCTCGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATCTATACTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGa                                                                                                                                        >  1:158194/1‑149 (MQ=255)           cTTCAGGTCAAAAAGTCCTGGTCATAGCACCTGCCCGTACTTCTCGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGcggcgg                                                                                                                                >  2:357013/1‑150 (MQ=255)           cTTCAGGTCAAAAAGTCCTGGTCATAGCACCTGCCCGTACTTCTCGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGcggcgg                                                                                                                                >  2:379635/1‑150 (MQ=255)            ttCAGGTCAAAAAGTCCTGGTCATAGCACCTGCCCGTACTTCTCGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGa                                                                                                                               >  2:301590/1‑150 (MQ=255)             tCAGGTCAAAAAGTCCTGGTCATAGCACCTGCCCGTACTTCTCGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGaa                                                                                                                              <  2:53698/150‑1 (MQ=255)             tCAGGTCAAAAAGTCCTGGTCATAGCACCTGCCCGTACTTCTCGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGaa                                                                                                                              <  1:587072/150‑1 (MQ=255)             tCAGGTCAAAAAGTCCTGGTCATAGCACCTGCCCGTACTTCTCGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGaa                                                                                                                              <  1:263896/150‑1 (MQ=255)                   cAAAAAGTCCTGGTCATAGCACCTGCCCGTACTTCTCGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTgg                                                                                                                       >  2:757006/1‑151 (MQ=255)                    aaaaaGTCCTGGTCATAGCACCTGCCCGTACTTCTCGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTgg                                                                                                                       <  2:769865/150‑1 (MQ=255)                        aGTCCTGGTCATAGCACCTGCCCGTACTTCTCGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAAc                                                                                                                  <  1:337488/151‑1 (MQ=255)                         gTCCTGGTCATAGCACCTGCCCGTACTTCTCGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACt                                                                                                                 >  1:231391/1‑151 (MQ=255)                           ccTGGTCATAGCACCTGCCCGTACTTCTCGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACttt                                                                                                               >  2:636557/1‑151 (MQ=255)                            cTGGTCATAGCACCTGCCCGTACTTCTCGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACtttt                                                                                                              >  2:397977/1‑151 (MQ=255)                                    aGCACCTGCCCGTACTTCTCGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTaaa                                                                                                          >  2:426035/1‑147 (MQ=255)                                       aCCTGCCCGTACTTCTCGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGAc                                                                                                   <  2:876204/151‑1 (MQ=255)                                       aCCTGCCCGTACTTCTCGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGAc                                                                                                   <  1:682626/151‑1 (MQ=255)                                           gCCCGTACTTCTCGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGcc                                                                                                >  2:355688/1‑150 (MQ=255)                                             ccGTACTTCTCGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAat                                                                                             >  2:26550/1‑151 (MQ=255)                                             ccGTACTTCTCGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAat                                                                                             <  2:300755/151‑1 (MQ=255)                                              cGTACTTCTCGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAata                                                                                            <  1:618116/151‑1 (MQ=255)                                               gTACTTCTCGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAatat                                                                                           >  1:317367/1‑151 (MQ=255)                                                 aCTTCTCGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTc                                                                                         >  2:158169/1‑151 (MQ=255)                                                       cGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGcc                                                                                    >  1:928672/1‑150 (MQ=255)                                                       cGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGc                                                                                     <  1:869643/149‑1 (MQ=255)                                                        gCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCAt                                                                                  <  1:155660/151‑1 (MQ=255)                                                             tGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAgc                                                                             >  1:455374/1‑151 (MQ=255)                                                                cGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAgcgcg                                                                          >  1:628704/1‑151 (MQ=255)                                                                  gTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGcgcgcg                                                                        >  1:993/1‑151 (MQ=255)                                                                     tCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCgcggcg                                                                     >  1:71660/1‑151 (MQ=255)                                                                     tCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCgcggcg                                                                     >  1:48943/1‑151 (MQ=255)                                                                      ccGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGcggcgg                                                                    <  1:821667/151‑1 (MQ=255)                                                                      ccGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGcggcgg                                                                    <  2:442922/151‑1 (MQ=255)                                                                          tACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCt                                                                >  2:169644/1‑151 (MQ=255)                                                                            cacTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGt                                                              <  1:108063/151‑1 (MQ=255)                                                                             acTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGtt                                                             >  1:658409/1‑151 (MQ=255)                                                                               tGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGtt                                                             >  1:54691/1‑149 (MQ=255)                                                                                 cATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGc                                                          >  2:18303/1‑150 (MQ=255)                                                                                  aTTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGc                                                          >  1:723367/1‑149 (MQ=255)                                                                                  aTTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAg                                                        <  2:295311/151‑1 (MQ=255)                                                                                      tGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACa                                                     <  1:226157/150‑1 (MQ=255)                                                                                             tACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTaa                                              <  2:595846/150‑1 (MQ=255)                                                                                               cATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGc                                           >  2:204374/1‑151 (MQ=255)                                                                                                 tttATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGcc                                          <  1:348676/150‑1 (MQ=255)                                                                                                 tttATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGcc                                          >  1:320533/1‑150 (MQ=255)                                                                                                    aTGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGt                                       <  2:754039/150‑1 (MQ=255)                                                                                                    aTGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGt                                       <  1:922679/150‑1 (MQ=255)                                                                                                        tGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCa                                   >  1:96272/1‑150 (MQ=255)                                                                                                        tGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCa                                   >  2:778215/1‑150 (MQ=255)                                                                                                         gAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCaaa                                 <  1:130531/151‑1 (MQ=255)                                                                                                            ggATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTg                              >  2:96984/1‑151 (MQ=255)                                                                                                               tatCCTCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCtt                           >  1:259081/1‑151 (MQ=255)                                                                                                                    tCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTg                      >  2:59743/1‑151 (MQ=255)                                                                                                                    tCATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTg                      >  1:398497/1‑151 (MQ=255)                                                                                                                     cATGTACCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCt                       <  2:761047/149‑1 (MQ=255)                                                                                                                          aCCAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCa                 >  1:222618/1‑150 (MQ=255)                                                                                                                           ccAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAAc               <  1:426020/151‑1 (MQ=255)                                                                                                                            cAGGATCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAAcg              >  2:275167/1‑151 (MQ=255)                                                                                                                                aTCTTATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAAcggcgg          >  2:704773/1‑151 (MQ=255)                                                                                                                                   ttATTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGtt        >  2:357704/1‑150 (MQ=255)                                                                                                                                     aTTCGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTccc     >  2:400443/1‑151 (MQ=255)                                                                                                                                        cGTAACGAACTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATAGTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGGGCTTTCCTGCGGCAACGGCGGGGCCCAtt  <  1:375601/151‑1 (MQ=255)                                                                                                                                               aaCTGCACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTccccct  >  1:493809/1‑141 (MQ=255)                                                                                                                                                    |                                                                                                                                           TTTATCCGCACTTCAGGTCAAAAAGTCCTGGTCATAGCACCTGCCCGTACTTCTCGCTTTTGGCGGTATCCGGTACACTGCATTTTGTCTATTACATTTATGCTGAAGGATATCCTCATGTACCAGGATCTTATTCGTAACGAACTGAACGAAGCGGCGGAAACGCTGGCTAACTTTTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATT  >  NC_000913/243426‑243711

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑