Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,271,827	0	A	T	40.3%	95.0 / 31.2	68	V8E (GTA→GAA)	chaA	calcium/sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base A (6/34);  new base T (27/0);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 2.70e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACACTAAAGGCGCTACTTAAAATACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAA‑CATTTTACGAATTTTTACCCCGCTGTCGATTTTTTACTCAT  >  NC_000913/1271679‑1271974                                                                                                                                                     |                                                                                                                                                       aacaCTAAAGGCGCTACTTAAAATACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTAcc                                                                                                                                                      <  2:247721/151‑1 (MQ=255)     cTAAAGGCGCTACTTAAAATACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCCGcc                                                                                                                                                   >  2:388995/1‑150 (MQ=255)       acagccgcTACTTAAAATACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGCAGCGTCTGGCCGCTTCACCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGACGTCTCCCTGCGCCGGCTCTTTACACCctct                                                                                                                                                >  1:854110/6‑151 (MQ=255)         aGGCGCTACTTAAAATACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTTTCGCCTCtt                                                                                                                                               >  1:123331/1‑150 (MQ=255)         aGGCGCTACTTAAAATACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTACGCCCCtt                                                                                                                                               >  1:689005/1‑150 (MQ=255)               tACTTAAAATACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTTTCGCCCCTTTAGCAtt                                                                                                                                        >  2:449265/1‑151 (MQ=255)                aCTTAAAATACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCGCCTCTTTAGGAttt                                                                                                                                       >  2:876440/1‑151 (MQ=255)                  ttAAAATACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCGCCCCTTTAGCAttt                                                                                                                                       >  2:998500/1‑149 (MQ=255)                   tAAAATACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTg                                                                                                                                      <  2:163869/149‑1 (MQ=255)                    aaaaTACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACa                                                                                                                                   <  2:966073/151‑1 (MQ=255)                     aaaTACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTTTCCCCTTTTTAGCATTTGTCa                                                                                                                                   >  1:334893/1‑150 (MQ=255)                      aaTACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACAtt                                                                                                                                 <  2:974489/151‑1 (MQ=255)                            aaTAAGCGCAAGAAGATTGATGGCAATGACCCCTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCACCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGtt                                                                                                                          <  1:78995/151‑1 (MQ=255)                            aaTAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGtt                                                                                                                          <  2:241093/151‑1 (MQ=255)                            aaTAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCCGCCTCTTTAGCATTTTTCATTA‑TTGt                                                                                                                           >  2:923553/1‑150 (MQ=255)                              tAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCCCCTCTTTAGCATTTTACAtttttggtat                                                                                                                         >  2:626138/1‑143 (MQ=255)                              tAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCCGACTCCTTAGCATTTGTAATTA‑ttggta                                                                                                                         >  1:365614/1‑145 (MQ=255)                                     aagaagATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTTCCGCCTCTTTAGCATTTTTCATTA‑TGGGTATTCCtttt                                                                                                                 >  1:407435/1‑150 (MQ=255)                                       gaagaTTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCa                                                                                                               <  1:273037/151‑1 (MQ=255)                                        aagaTTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAg                                                                                                              <  1:296001/151‑1 (MQ=255)                                        aagaTTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCGCCCCTTTGGCATTTTGCATTT‑TTGTTTTTCCttt                                                                                                                  >  2:276943/1‑147 (MQ=255)                                             tGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑a                                                                                                         <  2:820428/150‑1 (MQ=255)                                             tGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCGCCTCCTTAGGATTTTACATTT‑TGGTTATCCCttt                                                                                                                  >  2:968552/1‑142 (MQ=255)                                                 ggCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCTCTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTa                                                                                                    <  1:164710/151‑1 (MQ=255)                                                  gCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCCCCCCCTTTTGCATTTTTAATTT‑TTGTTATTCCTTTGCAGATG‑AAttt                                                                                                     >  1:1076955/1‑149 (MQ=255)                                                   cAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCCCCTCCTGAGCATTTTACATTA‑TTTTTATTCCTTTTCTGATT‑AAtttttt                                                                                                  >  1:360428/1‑148 (MQ=255)                                                    aaTGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTTTCCCCTCTTTTGCATTTTTCATTT‑TTGTTATTCCTTTTTCGATG‑AAttttt                                                                                                   >  1:632991/1‑147 (MQ=255)                                                       gACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGaaa                                                                                              <  1:473503/151‑1 (MQ=255)                                                       gACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGaa                                                                                               <  2:165108/150‑1 (MQ=255)                                                        aCCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGaaaa                                                                                             <  2:129050/151‑1 (MQ=255)                                                         ccACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGAAAAt                                                                                            <  2:579314/151‑1 (MQ=255)                                                           aCTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGAAAATGt                                                                                          <  2:354653/151‑1 (MQ=255)                                                              ggTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGAAAATGTaa                                                                                        <  2:881935/150‑1 (MQ=255)                                                              ggTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCCCCCTCTTTAGGATTTTACCATT‑TTGTTATTCCCTTTCAGATG‑AATTTTTTGaaaaaataa                                                                                        >  2:66430/1‑145 (MQ=255)                                                                 aGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGAAAATGTaaaaaa                                                                                    <  1:794023/151‑1 (MQ=255)                                                                  gtgtCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGAAAATGTAAAAAAt                                                                                   <  2:517631/151‑1 (MQ=255)                                                                    gtCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGAAAATGTAAAAAAta                                                                                  <  2:171409/150‑1 (MQ=255)                                                                       tGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTTTCGCCCCTTGTGCATTTTATATTT‑TGGTTTTCCCTTTTTAGATG‑AATTTTTTTAAAATGTAAAAAAta                                                                                  >  1:142350/1‑147 (MQ=255)                                                                       tGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCCCCCCCTTTAGCATATGATCTTT‑TTGTTTTTCCCTTGCGGCTG‑AATTTTTTGAAAATGTaaaaaaaa                                                                                  >  1:596111/1‑145 (MQ=255)                                                                        ggCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGAAAATGTAAAAAAtaggta                                                                              <  2:288195/150‑1 (MQ=255)                                                                        ggCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGAAAATGTAAAAAAtagg                                                                                <  2:800138/148‑1 (MQ=255)                                                                           tgctTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTATCCCCCCTTGAGCCTTTTTCATTT‑TTGTTTTTCCTTTTTAGATG‑AATTTTTTTAaaaaataaaaaataag                                                                                >  1:1047375/1‑143 (MQ=255)                                                                           tgctTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTAACCCCCCTTTAGCATTTTACCTTA‑TTGTTTTCCCCTTTCAGATG‑AAATTTTTGAAAAAGTaaaaaa                                                                                    >  1:532291/1‑141 (MQ=255)                                                                            gctTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGAAAATGTAAAAAATAGGTAGGa                                                                           <  2:102266/149‑1 (MQ=255)                                                                              ttCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTTTCGCCTCTTGTGGATTTTACATTT‑TTGTTTTCCCCTTTCCGAAG‑TATTTTTCACAAATGTTAAAAATAGGGAGGgaaa                                                                        >  2:359436/1‑150 (MQ=255)                                                                               tCCCCACAGGAACAGCATTACCGGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGAAAATGTAAAAAATAGGTAGGaaaaa                                                                       <  1:600041/150‑1 (MQ=255)                                                                                 cccACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATaa                                                                    <  2:721220/151‑1 (MQ=255)                                                                                     caGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAg                                                                 <  1:922094/150‑1 (MQ=255)                                                                                       ggAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGaaa                                                              <  1:923188/151‑1 (MQ=255)                                                                                           cAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTTCCGCCCCTTTAGCATTTTACATTA‑TTGTTATTCCTTTTCAGATG‑AATTTATCGAAAAAGTAAAAAAAAGGTAGGaaaaaaaaaaaaaa                                                              >  2:326433/1‑147 (MQ=255)                                                                                                cTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCCGCCCCTTTAGGATTTTACATTATTGGTTTTCCCTTT‑CCGATG‑AATTTTTTGAAAAAGTAAAAAAAAGGGAGGGAAAAAAAAAGaaaagggctgg                                                      >  2:874060/1‑142 (MQ=255)                                                                                                        gcCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTa                                             <  2:51012/151‑1 (MQ=255)                                                                                                         ccAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTaa                                            <  2:272216/151‑1 (MQ=255)                                                                                                             aaCCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTTTCCCCTCTTTTGCATTTTACATTA‑TTGTTTTCCCTTTTCAGAAG‑AATTTATTGAAAAAGTAAAAAAAAGGTAGGAAAAAAAAAAAAAAAAGGCTGGGTATCGGTAA‑AAtt                                       >  2:881086/1‑151 (MQ=255)                                                                                                                  ggAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTTCCCCCCTTTAACATTTGTCATTT‑TTGTTGTTCCCTTTTAGATG‑AAATTTTTTTAAATGTAAAAAATAAGGAGGGAAAATAAAAAAAAAAGTCTTGATATTGGTTA‑AAtttt                                     >  1:656310/1‑148 (MQ=255)                                                                                                                   gAAAATAAGCGAAGTCTCCTTGTGGCGGGGTTTTTCCGCCCCATATAAATTTTGAATTT‑TTGTTGTGTCTTTTCAGATG‑AAATTTTTTCAAATGTAAAAAAAAGGGAGGAAAAAAAAAAAAAAAAGTCTGGTTATTGGTAA‑AAttt                                      >  1:545765/1‑146 (MQ=255)                                                                                                                     aagtAAGCGGAGTCTCCTGGTGGCGCGTTTTTCCCGCCTCTTGAGCACTGGACAGTA‑GGGGTGAGTGTTTGCAGATG‑AATTTATCGAAAATGTAAAAAATAGTTAGGAAAAATAACAGAAAGTGTCTGTATATCGTTAA‑CATTTTACaaa                                <  1:649286/147‑1 (MQ=255)                                                                                                                         aaGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAA‑CATTTTACGAAttttt                           <  2:760275/151‑1 (MQ=255)                                                                                                                         aaGCGAAGTCTCCTTGTGGCGGATTTTTTTCGCCCCCTTAGTATATGACATAA‑TTGTTTTCCGTTTTTAGATGTAATTTTTCAAAAATTTAAAATA‑AGGTAGGTAAAAAAAAACAAAATGTCTTGGTAAACGTAATCATTTTA‑GAAttttt                           >  1:1000151/1‑151 (MQ=255)                                                                                                                              aaGTCTCCTTGTGGCGGGGTTTTTCCGCCTCTTGAGCATTTGACATTA‑TTGGTTTTCCTTTTTAGATG‑AATTTTTTTAAAAAATAAAAAAAAGGTAGGGAAAATAAAAAAAAATGTCTGGGTATTGGTAA‑AATTTTTAGAAatttttatcc                      >  1:363270/1‑146 (MQ=255)                                                                                                                                    ccTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAA‑CATTTTACGAATTTTTACCCCGCt                   <  1:716671/148‑1 (MQ=255)                                                                                                                                         tggtgggggtttACCGCCTCTTGAGCATTTGACATTA‑TGGCTATCGCTTTGTAGATG‑GATTTATGGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAA‑CATTTTAAGAAATTTTACCCCGCTGTCGAttt           >  1:487893/10‑151 (MQ=255)                                                                                                                                          ggCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAA‑CATTTTACGAATTTTTACCCCGCTGTCGAtttt          <  2:684933/151‑1 (MQ=255)                                                                                                                                           gCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGTTTATCCCTTTGCAGATG‑AATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAA‑CATTTTACGAATTTTTACCCCGCTGTCGAttttt         >  1:558545/1‑151 (MQ=255)                                                                                                                                            cGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAA‑CATTTTACGAATTTTTACCCCGCTGTCGAtttttt        <  2:764368/151‑1 (MQ=255)                                                                                                                                                tttttACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAA‑CATTTTACGAATTTTTACCCCGCTGTCGATTTTTTACt     <  1:388927/150‑1 (MQ=255)                                                                                                                                                tttttACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAA‑CATTTTACGAATTTTTACCCCGCTGTCGATTTTTTACTc    >  1:4309/1‑151 (MQ=255)                                                                                                                                                   ttACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAA‑CATTTTACGAATTTTTACCCCGCTGTCGATTTTTTACTCAt  <  1:758616/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                       AACACTAAAGGCGCTACTTAAAATACCAATAAGCGCAAGAAGATTGATGGCAATGACCACTGGTAGTGTCTGGCTGCTTCCCCACAGGAACAGCACTACCAGCGCCAGAACCGGGAAAATAAGCGAAGTCTCCTTGTGGCGGGTTTTTACCGCCTCTTGAGCATTTGACATTA‑TGGTTATCCCTTTGCAGATG‑AATTTATCGAAAATGTAAAAAATAGGTAGGAAAAATAACAGAAAGTGTCTGGATATCGGTAA‑CATTTTACGAATTTTTACCCCGCTGTCGATTTTTTACTCAT  >  NC_000913/1271679‑1271974

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