Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,180,945	A→G	5.8%	E270E (GAA→GAG)	PP_3678 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,180,945	0	A	G	5.8%	194.8 / 2.3	69	E270E (GAA→GAG)	PP_3678	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/31);  new base G (2/2);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.58e-01

CAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAAAGTCAATGCCCCCTCAGCACGCCAGTTCGATGCCTATTTAGACGCTAGAAAATCCTATATAGATTTTTTGCACTCGCGGGGTAT  >  NC_002947/4180801‑4181091                                                                                                                                                  |                                                                                                                                                   cAAGCACAACCCAGATCAATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑aaaaaa                                                                                                                                                <  3:70081/149‑1 (MQ=255)   aGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑aaaaaaaa                                                                                                                                              <  4:221933/149‑1 (MQ=255)    gCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑aaaaaa                                                                                                                                                >  6:208219/1‑146 (MQ=255)         aCCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAAtt                                                                                                                                        <  5:208219/149‑1 (MQ=255)          cccAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACATCCATGTATTAG‑AAAAAAAAGAAATTc                                                                                                                                       >  5:373765/1‑149 (MQ=255)            cAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATTAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCtt                                                                                                                                     <  3:274905/149‑1 (MQ=255)                 cTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTAcc                                                                                                                                >  2:76649/1‑149 (MQ=255)                 cTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTAcc                                                                                                                                >  3:214962/1‑149 (MQ=255)                  tataTGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGAAGATGAAACTAATGGGCATGCGACAAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTAcc                                                                                                                                >  5:289829/1‑148 (MQ=255)                   atatGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCt                                                                                                                               >  4:261482/1‑148 (MQ=255)                    tatGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTa                                                                                                                              >  5:172392/1‑148 (MQ=255)                      tGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACaa                                                                                                                           <  8:282476/149‑1 (MQ=255)                      tGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACaa                                                                                                                           <  2:219342/149‑1 (MQ=255)                        tCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAGAAAAAAAAAGAAATTCTTGTACCTACaa                                                                                                                           >  6:57838/1‑148 (MQ=255)                          gCCACAAATCTTTATATTTATGAAAGACGCATTTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAGAAAAAAAAAGAAATTCTTGTACCTACaaaa                                                                                                                         <  5:57838/148‑1 (MQ=255)                           ccACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCtt                                                                                                                                     <  6:282875/134‑1 (MQ=255)                           ccACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCtt                                                                                                                                     >  5:282875/1‑134 (MQ=255)                            cacaAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AGAAGAAAGAAATTCTTGt                                                                                                                                   <  6:11728/135‑1 (MQ=255)                            cacaAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AGAAGAAAGAAATTCTTGt                                                                                                                                   <  6:5883/135‑1 (MQ=255)                            cacaAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AGAAGAAAGAAATTCTTGt                                                                                                                                   >  5:5883/1‑135 (MQ=255)                            cacaAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AGAAGAAAGAAATTCTTGt                                                                                                                                   >  5:11728/1‑135 (MQ=255)                            cacaAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAATAAATTCTTGTACCTACAAAAAAAc                                                                                                                     >  7:273295/1‑149 (MQ=255)                             acaAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACg                                                                                                                    >  6:285987/1‑149 (MQ=255)                              caAATCTTTATATTTATGAAATACGCATGTTTATCATCATTATGATGAAACTCATGGGCATGCGACCATATGAACTATACGACCTTGACTTAAATGCCAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACgg                                                                                                                   <  1:313414/149‑1 (MQ=255)                              caAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACgg                                                                                                                   <  4:214962/149‑1 (MQ=255)                               aaaTCTTTATATTTATGCAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACggg                                                                                                                  <  6:410828/149‑1 (MQ=255)                               aaaTATTTATATTTATGTAAGACGCATGTTTATCATCAGGTTGATGAAACTCATGGGCATGCGACCAGAGTAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAATCggg                                                                                                                  <  6:373765/149‑1 (MQ=255)                                aaTCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACgggg                                                                                                                 <  8:252583/149‑1 (MQ=255)                                   cTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGTGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGaa                                                                                                               >  7:133017/1‑148 (MQ=255)                                        tatTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCt                                                                                                          <  1:76649/148‑1 (MQ=255)                                            tATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCAc                                                                                                      <  3:24514/148‑1 (MQ=255)                                               gAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCtt                                                                                                                                     >  6:115841/1‑114 (MQ=255)                                               gAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCtt                                                                                                                                     <  5:115841/114‑1 (MQ=255)                                                     cGCATGTTTATCATGAGTATTATGAAACTCATGGGCATGCGACCATAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCa                                                                                            <  8:273295/149‑1 (MQ=255)                                                      tcATGTTTATCATCAGGATGATGAAACTCATGGTCAAGCGACCAGAGGAACTATACGCCGATGACTAAAAGGAGAACCTCCATGTATTAG‑AAAAAATAGATATTCTTGTACGTACTATAAAACTGGGAATACCTGCACCCATACGCCaa                                                                                           >  6:426118/2‑148 (MQ=255)                                                      gCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCAc                                                                                           >  5:274367/1‑149 (MQ=255)                                                      gCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTAACTAAAAAAAAACGGGGAATACCTTCACCCATACGCCa                                                                                            >  6:404848/1‑148 (MQ=255)                                                        aTGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTAcc                                                                                                                                >  1:231783/1‑110 (MQ=255)                                                        aTGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTAcc                                                                                                                                <  2:231783/110‑1 (MQ=255)                                                        aTGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACtt                                                                                         <  3:261482/149‑1 (MQ=255)                                                          gTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACggg                                                                                                                  <  5:394999/122‑1 (MQ=255)                                                          gTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACggg                                                                                                                  >  6:394999/1‑122 (MQ=255)                                                          gTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTc                                                                                        >  5:102359/1‑148 (MQ=255)                                                             tATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCaaa                                                                                     <  6:274367/148‑1 (MQ=255)                                                                           tgaAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAAAg                                                                                    >  5:21992/1‑135 (MQ=255)                                                                           tgaAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAAAg                                                                                    <  6:21992/135‑1 (MQ=255)                                                                              aaCTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAAAGTCAATGCCCCCTCAGCa                                                                   >  6:208236/1‑149 (MQ=255)                                                                                 tCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAAAGTCAATGCCCCCTCAGCACGc                                                                >  4:102898/1‑149 (MQ=255)                                                                                       gCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCTACTTCAAAGTCAATGCCCCCTCAGCACGCCAGTTc                                                          <  5:285987/149‑1 (MQ=255)                                                                                         aTGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAAAGTCAATGCCCCCTCAGCACGCCAGTTCGa                                                        >  2:230038/1‑149 (MQ=255)                                                                                             gACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAAAGTCAATGCCCCCTCAGCACGCCAGTTCGATGcc                                                    >  2:195760/1‑149 (MQ=255)                                                                                               ccAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAAAGTCAATCCCCCCTCAGCACGCCAGTTCGATGCCTa                                                  <  1:230038/149‑1 (MQ=255)                                                                                                cAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAAAGTCAATGcccc                                                                          <  2:346608/124‑1 (MQ=255)                                                                                                  gagGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAAAGTCAATGCCCCCTCAGCACGCCAGTTCGATGCCTAtt                                                >  6:193680/1‑148 (MQ=255)                                                                                                      aaCTATACGACCTTGACTTAAAGTACAACCTCCATTTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAAAGTCAATGCCCCCTCAGCACGCCAGTTCGATGCCTATTTAGa                                            <  5:404848/148‑1 (MQ=255)                                                                                                       aCTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAAAGTCAATGCCCCCTCAGCACGCCAGTTCGATGCCTATTTAGACg                                          >  4:211350/1‑149 (MQ=255)                                                                                                        cTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAATGTCAATGCCCCCTCAGCATGCGAGTTCGATGGCTATTTAGACGc                                         >  6:142620/1‑149 (MQ=255)                                                                                                         tataCGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAAAGTCAATGCCCCCTCAGCACGCCAGTTCGATGCCTATTTAGACGCt                                        <  5:208236/149‑1 (MQ=255)                                                                                                           taCGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAAAGTCAATGCCCCCTCAGCACGCCAGTTCGATGCCTATTTAGACGc                                         >  7:128692/1‑146 (MQ=255)                                                                                                             cGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAAAGTCAATGCCCCCTCa                                                                      >  8:273075/1‑115 (MQ=255)                                                                                                             cGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAAAGTCAATGCCCCCTCa                                                                      <  7:273075/115‑1 (MQ=255)                                                                                                                 cTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAAAGTCAATGCCCCCTCAGCACGCCAGTTCGATGCCTATTTAGACGCTAGAAAATc                                >  5:290132/1‑149 (MQ=255)                                                                                                                    gACTTAAAGGACAACCTCCATTTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAAAGTCAATGCCCCCTCAGCACGCCAGTTCGATGCCTATTTAGACGCTAGAAAATcc                               <  6:290132/147‑1 (MQ=255)                                                                                                                    gACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAAAGTCAATGCCCCCTCAGCACGCCAGTTCGATGCCTATTTAGACGCTAGAAAATCCt                              <  3:102898/148‑1 (MQ=255)                                                                                                                      cTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAGATTCTTGTACCt                                                                                                                               <  3:38462/49‑1 (MQ=255)                                                                                                                      cTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCt                                                                                                                               >  4:38462/1‑49 (MQ=255)                                                                                                                                aaCCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAAAGTCCATGCCCCCTCAGCACGCCAGTTCGATGCCTATTTAGACGCTAGAAAATCCTATATAGAtttttt                 <  6:102359/149‑1 (MQ=255)                                                                                                                                  ccTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAAAGTCAATGCCCCCTCAGCACGCCAGTTCGATGCCTATTTAGACGCTAGAAAATCCTATATAGATTTTTTg                >  4:70560/1‑148 (MQ=255)                                                                                                                                   cTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAAAGTCAATGCCCCCTCAGCACGCCAGTTCGAt                                                       >  3:67541/1‑108 (MQ=255)                                                                                                                                   cTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAAAGTCAATGCCCCCTCAGCACGCCAGTTCGAt                                                       <  4:67541/108‑1 (MQ=255)                                                                                                                                       aTGTATTAG‑AAAAAAAAGAAATTCTAGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAAAGTAAATGCCCCCTCAGCACGCAAATTCGATGCCTATATAGACGCTAGAAAATCATATATAGATTATATGCaatag          >  6:212009/1‑145 (MQ=255)                                                                                                                                               g‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAAAGTCAATGCCCCCTCAGCACGCCAGTTCGATGCCTATTTAGACGCTAGAAAATCCTATATAGATTTTTTGCACTCGCGGGGTAt  >  4:92792/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                   CAAGCACAACCCAGATCTATATGTCGCCACAAATCTTTATATTTATGAAAGACGCATGTTTATCATCAGGATGATGAAACTCATGGGCATGCGACCAGAGGAACTATACGACCTTGACTTAAAGGACAACCTCCATGTATTAG‑AAAAAAAAGAAATTCTTGTACCTACAAAAAAACGGGGAATACCTGCACCCATACGCCACTTCAAAGTCAATGCCCCCTCAGCACGCCAGTTCGATGCCTATTTAGACGCTAGAAAATCCTATATAGATTTTTTGCACTCGCGGGGTAT  >  NC_002947/4180801‑4181091

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