Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,021,557 T→C 6.5% P444P (CCT→CCC ubiB → regulator of octaprenylphenol hydroxylation, ubiquinone synthesis; regulator of 2'‑N‑acetyltransferase; putative ABC1 family protein kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,021,5570TC 6.5% 186.2 / 3.7 61P444P (CCT→CCCubiBregulator of octaprenylphenol hydroxylation, ubiquinone synthesis; regulator of 2'‑N‑acetyltransferase; putative ABC1 family protein kinase
Reads supporting (aligned to +/- strand):  ref base T (33/24);  new base C (2/2);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

TGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCGCCAGGGCAAGTATTTACAGCACAGTGTTGATAAGATTGCCCGCGAGCTTCAGTCAAAT  >  NC_000913/4021412‑4021701
                                                                                                                                                 |                                                                                                                                                
tGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTgc                                                                                                                                                <  7:108795/148‑1 (MQ=255)
       gCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGt                                                                                                                                          >  1:133216/1‑147 (MQ=255)
       gCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGt                                                                                                                                          >  4:152883/1‑147 (MQ=255)
       gCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCCGCGCTGGt                                                                                                                                          >  3:57127/1‑147 (MQ=255)
        cAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTg                                                                                                                                         <  2:133216/147‑1 (MQ=255)
                    gtgtTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTa                                                                                                                               >  7:188672/1‑145 (MQ=255)
                        tACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGaaa                                                                                                                         >  8:15921/1‑147 (MQ=255)
                            ccAGAAAACCCTGCTCTACGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTga                                                                                                                                        <  4:180776/128‑1 (MQ=255)
                            ccAGAAAACCCTGCTCTACGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTga                                                                                                                                        >  3:180776/1‑128 (MQ=255)
                            ccAGAAAACCCTGCTCTACGTCGAAGGGGAAGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAgcg                                                                                                                    >  6:301666/1‑148 (MQ=255)
                               gAAAACCCTGCTCTACGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCg                                                                                                                 >  1:365190/1‑148 (MQ=255)
                                     ccTGCTCTACGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTgg                                                                                                           >  7:188065/1‑148 (MQ=255)
                                       tGCTCTACGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGt                                                                                                         >  8:82212/1‑148 (MQ=255)
                                               gTCGAAGGGGGAGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGaaaaaa                                                                                                 <  1:57233/148‑1 (MQ=255)
                                                tCGAAGGGGTAGGACGCCAGCTTTATCCGCAACTCGATTTCTGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGaaaa                                                                                                   <  5:100399/145‑1 (MQ=255)
                                                 cGAAGGGGTAGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGCGCATTTAAAGAAAAAGCGCCGTTCGGGGTCGcaaaa                                                                                                  >  6:100399/1‑145 (MQ=255)
                                                      gggTAGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGa                                                                                          <  8:188672/148‑1 (MQ=255)
                                                      gggTAGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGa                                                                                          <  8:188065/148‑1 (MQ=255)
                                                       ggTAGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGaaaa                                                                                                   <  8:128726/138‑1 (MQ=255)
                                                       ggTAGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGaaaa                                                                                                   >  7:128726/1‑138 (MQ=255)
                                                           ggACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGcc                                                                                                                  <  6:184380/119‑1 (MQ=255)
                                                           ggACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGcc                                                                                                                  >  5:184380/1‑119 (MQ=255)
                                                             aCGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTACTGCGCTGGTGAGAGCATTTAAAGAAAAAGAGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCt                                                                                   >  7:317937/1‑148 (MQ=255)
                                                                    cTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTgg                                                                            >  1:323581/1‑148 (MQ=255)
                                                                          ccGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTAc                                                                       >  3:134012/1‑147 (MQ=255)
                                                                             cAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGaaaaa                                                                                                  <  5:136397/117‑1 (MQ=255)
                                                                             cAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGaaaaa                                                                                                  >  6:136397/1‑117 (MQ=255)
                                                                               aCTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGt                                                                 <  3:152883/148‑1 (MQ=255)
                                                                                    aTTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGaaa                                                                                                                         >  1:4866/1‑87 (MQ=255)
                                                                                    aTTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGaaa                                                                                                                         <  2:4866/87‑1 (MQ=255)
                                                                                    aTTTAGGGAAAAAGGCGAAGCCTTTCCTGGAGTCGTTGATTAAAGATCAGGTCGGTATTCCTGCGCTGGGGAGAGCCTTTAAAGAAAAAGCGCCGTTCTGGTTCGAAAAAATTCCATAACTGCCTGAATTGGTTTACTACAGTTTgcg                                                            <  8:317937/148‑1 (MQ=255)
                                                                                       tATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCGcc                                                          <  2:365190/147‑1 (MQ=255)
                                                                                           gAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGACCTGCCTGAATTGGTTTACGACAGTTTGCGCCAggg                                                      >  5:6584/1‑147 (MQ=255)
                                                                                                      aGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCGCCAGGGCAAGTATTTa                                            >  5:125769/1‑146 (MQ=255)
                                                                                                         cTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCAGAACTGCCTGAATTGGTTTACGAc                                                                    >  3:26877/1‑118 (MQ=255)
                                                                                                         cTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCAGAACTGCCTGAATTGGTTTACGAc                                                                    <  4:26877/118‑1 (MQ=255)
                                                                                                           ttCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCGCCAGGGCAAGTATTTACAGca                                       >  7:165476/1‑146 (MQ=255)
                                                                                                           ttCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCGCCAGGGCAAGTATTTACAGca                                       >  6:134621/1‑146 (MQ=255)
                                                                                                           atacTGGACTAGTGGATTAAAGAAAACGTAGGTATTACCGCGATGGTGAGAGCATTTAAAGAAAAACCGCCGATCTCGCTAGAAAAAATGCCAGAACTGCCTGAATTGGTTTACTACAGTTTGCGCCAGGCCAAGTATTTACAGcaca                                     >  8:287355/4‑148 (MQ=255)
                                                                                                            tCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCGCCAGGGCAAGTATTTACAGcac                                      <  7:82212/146‑1 (MQ=255)
                                                                                                             ccTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCCGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAATAATTGCAAGAACTGCGTGAATTGGTTTACGACAGTTTGCGCCAGGGCAAGTATTTACAGCACAg                                    <  7:287355/147‑1 (MQ=255)
                                                                                                                ggAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCGCCAGGGCAAGTATTTACAGCACAGTGtt                                >  1:47116/1‑148 (MQ=255)
                                                                                                                 gAGTCGTGGATTAAAGATCAGGTCGGTATTCCCGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCGCCAGGGCAAGTATTTACAGCACAGTGTTg                               <  4:57127/148‑1 (MQ=255)
                                                                                                                    tCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCGCCAGGGCAAGTATTTACAGCACAGTGTTGATa                            <  7:15921/148‑1 (MQ=255)
                                                                                                                      gTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCGCCAGGGCAAGTATTTACAGCACAGTGTTGATAAg                          <  8:165476/148‑1 (MQ=255)
                                                                                                                       tGGATTAAAGCTCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCGCCAGGGCAAGTATTTACAGCACAGTGTTGATAAGa                         <  2:323581/148‑1 (MQ=255)
                                                                                                                       tGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCGCCAGGGCAAGTATTTACAGCACAGTGTTGAt                             >  7:115166/1‑144 (MQ=255)
                                                                                                                       tGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCGCCAGGGCAAGTATTTACAGCACAGTGTTGAt                             >  7:371999/1‑144 (MQ=255)
                                                                                                                            tAAAGATCAGGTCGGTATTCCTGCTCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCTCCAGGCAAGTATTTACAGCACAGTGTTGATAAGATTGccc                   >  4:14479/1‑148 (MQ=255)
                                                                                                                             aaaGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCGCCAGGGCAAGTATTTACAGCACAGTGTTGATAAGATTGcc                    <  8:371999/147‑1 (MQ=255)
                                                                                                                                  tCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCGCCAGGGCAAGTATTTACAGCACAGTGTTGATAAGATTGCCCGCGAg              >  3:127531/1‑148 (MQ=255)
                                                                                                                                       tCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCGcc                                                          >  5:376265/1‑99 (MQ=255)
                                                                                                                                       tCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCGcc                                                          <  6:376265/99‑1 (MQ=255)
                                                                                                                                       tCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCGCCAGGGCAAGTATTTACAGCACAGTGTTGATAAGATTg                      >  8:103641/1‑135 (MQ=255)
                                                                                                                                       tCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCGCCAGGGCAAGTATTTACAGCACAGTGTTGATAAGATTg                      <  7:103641/135‑1 (MQ=255)
                                                                                                                                       tCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCGCCAGGGCAAGTATTTACAGCACAGTGTTGATAAGATTGCCCGCGAGCtt           >  8:271100/1‑146 (MQ=255)
                                                                                                                                        cGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCGCCAGGGCAAGTATTTACAGCACAGTGTTGATAAGATTGCCCGCGAGCTTCAg        <  2:47116/148‑1 (MQ=255)
                                                                                                                                           tATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCGCCAGGGCAAGTATTTACAGCACAGTGTTGATAAGATTGCCCGCGAGCTTCAGTc      >  6:222058/1‑147 (MQ=255)
                                                                                                                                             ttCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCGCCAGGGCAAGTATTTACAGCACAGTGTTGATAAGATTGCCCGCGAGCTTCAGTCaa    <  5:301666/147‑1 (MQ=255)
                                                                                                                                              tCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCGCCAGGGCAAGTATTTACAGCACAGTGTTGATAAGATTGCCCGCGAGCTTCAGTCAAAt  >  1:81681/1‑148 (MQ=255)
                                                                                                                                              tCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCGCCAGGGCAAGTATTTACAGCACAGTGTTGATAAGATTGCCCGCGAGCTTCAGTCAAAt  >  6:226940/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TGGAAGTGCAGCCGCAACTGGTGTTACTCCAGAAAACCCTGCTCTACGTCGAAGGGGTAGGACGCCAGCTTTATCCGCAACTCGATTTATGGAAAACGGCGAAGCCTTTCCTGGAGTCGTGGATTAAAGATCAGGTCGGTATTCCTGCGCTGGTGAGAGCATTTAAAGAAAAAGCGCCGTTCTGGGTCGAAAAAATGCCAGAACTGCCTGAATTGGTTTACGACAGTTTGCGCCAGGGCAAGTATTTACAGCACAGTGTTGATAAGATTGCCCGCGAGCTTCAGTCAAAT  >  NC_000913/4021412‑4021701

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