Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,803,324	(T)5→4	pseudogene (231/426 nt)	arpB →	pseudogene, ankyrin repeats

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,803,320	0	T	.	100.0%	315.0 / NA	65	pseudogene (227/426 nt)	arpB	pseudogene, ankyrin repeats
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (24/41);  total (24/41)

ACCTTCATCAGAGCCTATTAATCAATTTGCGCGCCAGCTCATCACCCTACTTGATACCTCTGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCCTTTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGAAATCATTTCCTTATTACTGGAAACCACGCCTATTAAAGCTAACCATTTGCGAAGCATACTGGATTATATTGA  >  NC_000913/1803171‑1803464                                                                                                                                                      |                                                                                                                                                 acCTTCATCAGAGCCTATTAATCAATTTGCGCGCCAGCTCATCACCCTACTTGATACCTCTGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCCtt                                                                                                                                                 <  1:593606/151‑1 (MQ=255)     tcatcaGAGCCTATTAATCAATTTGCGCGCCAGCTCATCACCCTACTTGATACCTCTGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCCttttac                                                                                                                                             >  1:471501/1‑149 (MQ=255)            agCCTATTAATCAATTTGCGCGCCAGCTCATCACCCTACTTGATACCTCTGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAAc                                                                                                                                     <  2:3910/151‑1 (MQ=255)             gCCTATTAATCAATTTGCGCGCCAGCTCATCACCCTACTTGATACCTCTGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACt                                                                                                                                    <  1:883846/151‑1 (MQ=255)             gCCTATTAATCAATTTGCGCGCCAGCTCATCACCCTACTTGATACCTCTGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACt                                                                                                                                    >  1:408189/1‑151 (MQ=255)              ccTATTAATCAATTTGCGCGCCAGCTCATCACCCTACTTGATACCTCTGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACtt                                                                                                                                   >  2:187576/1‑151 (MQ=255)              ccTATTAATCAATTTGCGCGCCAGCTCATCACCCTACTTGATACCTCTGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACtt                                                                                                                                   >  2:81882/1‑151 (MQ=255)                 aTTAATCAATTTGCGCGCCAGCTCATCACCCTACTTGATACCTCTGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACttt                                                                                                                                  <  1:136753/149‑1 (MQ=255)                 aTTAATCAATTTGCGCGCCAGCTCATCACCCTACTTGATACCTCTGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGc                                                                                                                                <  1:105975/151‑1 (MQ=255)                  ttAATCAATTTGCGCGCCAGCTCATCACCCTACTTGATACCTCTGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGc                                                                                                                                >  2:498574/1‑150 (MQ=255)                  ttAATCAATTTGCGCGCCAGCTCATCACCCTACTTGATACCTCTGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGc                                                                                                                                <  1:498623/150‑1 (MQ=255)                    aaTCAATTTGCGCGCCAGCTCATCACCCTACTTGATACCTCTGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCAc                                                                                                                             <  1:626022/151‑1 (MQ=255)                        aaTTTGCGCGCCAGCTCATCACCCTACTTGATACCTCTGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCa                                                                                                                         <  1:720551/151‑1 (MQ=255)                          tttGCGCGTCAGCTCATCACCCTACTTGATACCTCTGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCAcc                                                                                                                       >  2:677504/1‑151 (MQ=255)                                 gcCAGCTCATCACCCTACTTGATACCTCTGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACccatccat                                                                                                                <  2:159743/151‑1 (MQ=255)                                  ccAGCTCATCACCCTACTTGATACCTCTGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATa                                                                                                               <  1:432573/151‑1 (MQ=255)                                   cAGCTCATCACCCTACTTGATACCTCTGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATAc                                                                                                              <  2:448602/151‑1 (MQ=255)                                      cTCATCACCCTACTTGATACCTCTGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATAc                                                                                                              >  1:764175/1‑148 (MQ=255)                                        catcaCCCTACTTGATACCTCTGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGaaa                                                                                                         <  1:500970/151‑1 (MQ=255)                                                 aCTTGATACCTCTGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTca                                                                                                <  1:108707/151‑1 (MQ=255)                                                     gATACCTCTGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTcatta                                                                                             <  2:582542/150‑1 (MQ=255)                                                        aCCTCTGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTcattacatt                                                                                         <  1:93743/151‑1 (MQ=255)                                                         ccTCTGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTcattacatt                                                                                         >  2:774704/1‑150 (MQ=255)                                                          ctctGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTcattacatta                                                                                        >  1:587906/1‑150 (MQ=255)                                                            ctGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTcattacatta                                                                                        <  2:587857/148‑1 (MQ=255)                                                            ctGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCattacattac                                                                                       <  1:165684/149‑1 (MQ=255)                                                            ctGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCattacattac                                                                                       <  1:509384/149‑1 (MQ=255)                                                            ctGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTacac                                                                                     >  1:883424/1‑151 (MQ=255)                                                              gACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACt                                                                                    >  2:313422/1‑150 (MQ=255)                                                               aCTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACta                                                                                   <  1:668775/150‑1 (MQ=255)                                                                  tAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACtata                                                                                 <  1:323660/149‑1 (MQ=255)                                                                  tAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACtata                                                                                 <  2:184303/149‑1 (MQ=255)                                                                  tAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATaa                                                                                <  1:758004/150‑1 (MQ=255)                                                                   aaGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATaaaa                                                                              <  2:612105/151‑1 (MQ=255)                                                                     gTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATaaaaa                                                                             >  1:698152/1‑150 (MQ=255)                                                                     gTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATaaaaa                                                                             >  2:663701/1‑150 (MQ=255)                                                                        cGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGaa                                                                         <  1:47487/151‑1 (MQ=255)                                                                        cGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGaa                                                                         <  2:102927/151‑1 (MQ=255)                                                                        cGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGaa                                                                         <  1:210392/151‑1 (MQ=255)                                                                        cGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGaa                                                                         <  2:886680/151‑1 (MQ=255)                                                                        cGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAAGATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGa                                                                          <  2:494457/150‑1 (MQ=255)                                                                          aTGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTcattacat                                                                                          >  1:583914/1‑132 (MQ=255)                                                                          aTGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTcattacat                                                                                          <  2:583864/132‑1 (MQ=255)                                                                            gCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGAAATCa                                                                     <  2:764125/151‑1 (MQ=255)                                                                                 cATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGAAATCATTTc                                                                 <  2:346751/150‑1 (MQ=255)                                                                                  aTACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGAAATCATTTCCt                                                               >  1:337544/1‑151 (MQ=255)                                                                                          ttACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGAAATCATTTCCTTATTACTg                                                       >  1:859731/1‑151 (MQ=255)                                                                                               cAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGAAATCATTTCCTTATTACTGGaa                                                    <  2:303369/149‑1 (MQ=255)                                                                                               cAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGAAATCATTTCCTTATTACTGGAAAc                                                  >  1:316227/1‑151 (MQ=255)                                                                                               cAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGAAATCATTTCCTTATTACTGGAAAc                                                  >  1:59842/1‑151 (MQ=255)                                                                                                 ggAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGAAATCATTTCCTTATTACTGGAAACCa                                                <  1:464598/151‑1 (MQ=255)                                                                                                         ccGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGAAATCATTTCCTTATTACTGGAAACCACGCCTAtt                                        >  2:890288/1‑151 (MQ=255)                                                                                                         ccGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGAAATCATTTCCTTATTACTGGAAACCACGCCTAtt                                        <  2:880250/151‑1 (MQ=255)                                                                                                         ccGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGAAATCATTTCCTTATTACTGGAAACCACGCCTAt                                         <  2:882029/150‑1 (MQ=255)                                                                                                          cGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGAAATCATTTCCTTATTACTGGAAACCACGCCTATTa                                       <  2:723033/151‑1 (MQ=255)                                                                                                           gCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGAAATCATTTCCTTATTACTGGAAACCACGCCTATTa                                       >  1:804866/1‑150 (MQ=255)                                                                                                            cAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGAAATCATTTCCTTATTACTGGAAACCACGCCTATTaaa                                     <  2:60391/151‑1 (MQ=255)                                                                                                              aaCTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGAAATCATTTCCTTATTACTGGAAACCACGCCTATTAAAGc                                   >  1:762460/1‑151 (MQ=255)                                                                                                                    cGAAAAATATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGAAATCATTTCCTTATTACTGGAAACCACGCCTATTAAAGCTAACCa                             <  1:616741/151‑1 (MQ=255)                                                                                                                      aaaaaTATCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGAAATCATTTCCTTATTACTGGAAACCACGCCTATTAAAGCTAACCAt                            <  2:641245/150‑1 (MQ=255)                                                                                                                            atCACAAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGAAATCATTTCCTTATTACTGGAAACCACGCCTATTAAAGCTAACCATTTGCGaa                     <  1:733864/151‑1 (MQ=255)                                                                                                                                caAAATTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGAAATCATTTCCTTATTACTGGAAACCACGCCTATTAAAGCTAACCATTTGCGAAGCa                  <  2:314021/150‑1 (MQ=255)                                                                                                                                    aTTGTTATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGAAATCATTTCCTTATTACTGGAAACCACGCCTATTAAAGCTAACCATTTGCGAAGCATACTg             >  1:669778/1‑151 (MQ=255)                                                                                                                                        ttATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGAAATCATTTCCTTATTACTGGAAACCACGCCTATTAAAGCTAACCATTTGCGAAGCATACTGGAtt         >  2:8266/1‑151 (MQ=255)                                                                                                                                         tATTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGAAATCATTTCCTTATTACTGGAAACCACGCCTATTAAAGCTAACCATTTGCGAAGCATACTGGAtt         <  1:193307/150‑1 (MQ=255)                                                                                                                                          aTTCCACTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGAAATCATTTCCTTATTACTGGAAACCACGCCTATTAAAGCTAACCATTTGCGAAGCATACTGGATtat       >  1:543870/1‑151 (MQ=255)                                                                                                                                                cTGCCC‑TTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGAAATCATTTCCTTATTACTGGAAACCACGCCTATTAAAGCTAACCATTTGCGAAGCATACTGGATTATATTga  >  2:431994/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                 ACCTTCATCAGAGCCTATTAATCAATTTGCGCGCCAGCTCATCACCCTACTTGATACCTCTGACTTAAGTTCGATGCTGTCATACTGTGTTACTCAGGAATTTACCGCAAACTGTCGAAAAATATCACAAAATTGTTATTCCACTGCCCTTTTTACCATTAACTTTGCCACTTCACCCATCCATACAGAAAATATACTCATTACATTACACTATAAAAAAGAAATCATTTCCTTATTACTGGAAACCACGCCTATTAAAGCTAACCATTTGCGAAGCATACTGGATTATATTGA  >  NC_000913/1803171‑1803464

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