Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,971,151 T→C 13.3% intergenic (+18/‑120) ygdQ → / → ygdR UPF0053 family inner membrane protein/DUF903 family verified lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,971,1510TC13.3% 237.9 / 22.6 82intergenic (+18/‑120)ygdQ/ygdRUPF0053 family inner membrane protein/DUF903 family verified lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (35/36);  new base C (5/6);  total (40/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

GCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATA  >  NC_000913/2971009‑2971265
                                                                                                                                              |                                                                                                                  
gCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCggg                                                                                                               <  4:622066/148‑1 (MQ=255)
    tACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAgag                                                                                                           <  4:1053371/148‑3 (MQ=255)
         tGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAGCAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTac                                                                                                      >  2:895956/1‑148 (MQ=255)
                 gAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTactc                                                                                                    >  2:372401/1‑142 (MQ=255)
                 gAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTactc                                                                                                    <  1:372401/142‑1 (MQ=255)
                 gAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTa                                                                                                       <  2:807520/139‑1 (MQ=255)
                 gAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTa                                                                                                       >  1:807520/1‑139 (MQ=255)
                 gAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCt                                                                                              >  4:182310/1‑148 (MQ=255)
                  aaaGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCtt                                                                                             >  3:322307/1‑148 (MQ=255)
                            cATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTcac                                                                                   <  4:322307/148‑1 (MQ=255)
                            cATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTcac                                                                                   <  3:182310/148‑1 (MQ=255)
                                    tACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAg                                                                           <  1:895956/148‑1 (MQ=255)
                                               tACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCt                                                                >  4:206537/1‑148 (MQ=255)
                                               tACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAgag                                                                                                           <  1:738279/105‑3 (MQ=255)
                                               tACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAgag                                                                                                           >  2:738279/1‑103 (MQ=255)
                                                   tCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCa                                                            >  1:706647/1‑148 (MQ=255)
                                                   tCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCa                                                            <  2:706647/148‑1 (MQ=255)
                                                     tACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAgg                                                                                                            >  1:7742/1‑98 (MQ=255)
                                                     tACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAgg                                                                                                            <  2:7742/98‑1 (MQ=255)
                                                      aCTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGa                                                         <  3:206537/148‑1 (MQ=255)
                                                         tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAggg                                                                                                           <  1:854101/95‑1 (MQ=255)
                                                         tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAggg                                                                                                           >  2:854101/1‑95 (MQ=255)
                                                         tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGc                                                                                                          <  1:524320/96‑1 (MQ=255)
                                                         tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGc                                                                                                          >  2:524320/1‑96 (MQ=255)
                                                         tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTactca                                                                                                   <  1:61342/103‑1 (MQ=255)
                                                         tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTactca                                                                                                   >  2:61342/1‑103 (MQ=255)
                                                         tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTact                                                                                                     >  1:277271/1‑101 (MQ=255)
                                                         tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTact                                                                                                     <  2:277271/101‑1 (MQ=255)
                                                         tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTa                                                                                                       <  3:124091/99‑1 (MQ=255)
                                                         tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTa                                                                                                       >  4:124091/1‑99 (MQ=255)
                                                              aTGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAg                                                          <  1:1125022/139‑1 (MQ=255)
                                                              aTGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAg                                                          >  2:1125022/1‑139 (MQ=255)
                                                               tGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAggg                                                                                                           <  3:392036/89‑1 (MQ=255)
                                                               tGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAggg                                                                                                           >  4:392036/1‑89 (MQ=255)
                                                                     tctcTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTa                                                                             >  2:1053637/1‑113 (MQ=255)
                                                                     tctcTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTa                                                                             <  1:1053637/113‑1 (MQ=255)
                                                                     tctcTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTa                                                                                                       <  2:207265/87‑1 (MQ=38)
                                                                     tctcTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTa                                                                                                       >  1:207265/1‑87 (MQ=38)
                                                                     tctcTATTGCGGTT‑‑AAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTta                                        <  1:128813/148‑1 (MQ=255)
                                                                     tctcTATTGCGGTT‑‑AAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTta                                        >  2:128813/1‑148 (MQ=255)
                                                                          aTTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTactca                                                                                                   >  4:641643/1‑86 (MQ=255)
                                                                          aTTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTactca                                                                                                   <  3:641643/86‑1 (MQ=255)
                                                                            tGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGa                                                                          >  2:960018/1‑109 (MQ=255)
                                                                            tGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGa                                                                          <  1:960018/109‑1 (MQ=255)
                                                                                     aaGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCa                                                            >  1:16240/1‑114 (MQ=255)
                                                                                     aaGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCa                                                            <  2:16240/114‑1 (MQ=255)
                                                                                        ccTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCg                                                                                                         >  3:737580/1‑66 (MQ=255)
                                                                                        ccTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCa                                                            >  1:648777/1‑111 (MQ=255)
                                                                                        ccTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCa                                                            <  2:648777/111‑1 (MQ=255)
                                                                                        acTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCg                                                                                                         <  4:737580/65‑1 (MQ=255)
                                                                                         cTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCaa                                                           <  1:581167/111‑1 (MQ=255)
                                                                                         cTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCaa                                                           >  2:581167/1‑111 (MQ=255)
                                                                                          tCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCt                                                                                              >  1:531646/1‑75 (MQ=255)
                                                                                          tCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCt                                                                                              <  2:531646/75‑1 (MQ=255)
                                                                                          tCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAc                                                                                         <  2:311657/80‑1 (MQ=255)
                                                                                          tCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAc                                                                                         >  1:311657/1‑80 (MQ=255)
                                                                                               ttGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGata                                                       >  3:366361/1‑109 (MQ=255)
                                                                                               ttGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGata                                                       <  4:366361/109‑1 (MQ=255)
                                                                                                tGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAgag                                                                                                           >  3:131695/1‑54 (MQ=25)
                                                                                                tGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAgag                                                                                                           <  4:131695/56‑3 (MQ=25)
                                                                                                    tCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtaggta                 >  4:973696/1‑142 (MQ=255)
                                                                                                    tCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtaggta                 <  3:973696/142‑1 (MQ=255)
                                                                                                       caacaaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGttt                                                                              <  2:44353/78‑1 (MQ=255)
                                                                                                       caacaaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGttt                                                                              >  1:44353/1‑78 (MQ=255)
                                                                                                       caacaaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCgg                                            <  3:233012/112‑1 (MQ=255)
                                                                                                       caacaaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCgg                                            >  4:233012/1‑112 (MQ=255)
                                                                                                           aaaaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGtt                                                                               <  2:2574/73‑1 (MQ=255)
                                                                                                           aaaaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGtt                                                                               >  1:2574/1‑73 (MQ=255)
                                                                                                            aaaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTa                                                                             >  3:819295/1‑74 (MQ=255)
                                                                                                            aaaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTa                                                                             <  4:819295/74‑1 (MQ=255)
                                                                                                             aaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGa       >  3:335031/1‑143 (MQ=255)
                                                                                                             aaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGa       <  4:335031/143‑1 (MQ=255)
                                                                                                                 aaTCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTta                                        >  3:677171/1‑106 (MQ=255)
                                                                                                                 aaTCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTta                                        <  4:677171/106‑1 (MQ=255)
                                                                                                                                tCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATa  <  4:828878/129‑1 (MQ=255)
                                                                                                                                tCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATa  >  3:828878/1‑129 (MQ=255)
                                                                                                                                   gTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAg                                 >  2:981382/1‑95 (MQ=38)
                                                                                                                                   gTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAg                                 <  1:981382/95‑1 (MQ=38)
                                                                                                                                     aCGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtagg                   >  3:749781/1‑107 (MQ=255)
                                                                                                                                     aCGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtagg                   <  4:749781/107‑1 (MQ=255)
                                                                                                                                      cGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATg        <  4:583881/117‑1 (MQ=255)
                                                                                                                                      cGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATg        >  3:583881/1‑117 (MQ=255)
                                                                                                                                              tGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAg                                                          >  2:254111/1‑59 (MQ=255)
                                                                                                                                              tGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAg                                                          <  1:254111/59‑1 (MQ=255)
                                                                                                                                              |                                                                                                                  
GCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATA  >  NC_000913/2971009‑2971265

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