Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,060,651 T→C 7.1% intergenic (+5/‑199) argP → / → scpA DNA‑binding transcriptional dual regulator ArgP/methylmalonyl‑CoA mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,060,6510TC 7.1% 172.0 / 3.0 66intergenic (+5/‑199)argP/scpADNA‑binding transcriptional dual regulator ArgP/methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base T (37/24);  new base C (2/3);  total (39/27)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

GATTGACTTAACGCCTGGGCTATTTCAACGACGGATGCTCTACTGGCACCGCTTTGCTCCTGAAAGCCGCATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTGCTGGCGACAAACCCACGTA  >  NC_000913/3060514‑3060792
                                                                                                                                         |                                                                                                                                             
gATTGACTTAACGCCTGGGCTATTTCAACGACGGATGCTCTACTGGCACCGCTTTGCTCCTGAAAGCCGCATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAaat                                                                                                                                          >  1:821871/1‑143 (MQ=255)
     aCTTAACGCCTGGGCTATTTCAACGACGGATGCTCTACTGGCACCGCTTTGCTCCTGAAAGCCGCATGATGCGTAAGGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCCAATAATGc                                                                                                                                     >  2:159866/1‑143 (MQ=255)
       ttAACGCCTGGGCTATTTCAACGACGGATGCTCTACTGGCACCGCTTTGCTCCTGAAAGCCGCATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCt                                                                                                                                   >  2:843106/1‑143 (MQ=255)
           cGCCTGGGCTATTTCAACGACGGATGCTCTACTGGCACCGCTTTGCTCCTGAAAGCCGCATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATa                                                                                                                               <  1:356529/143‑1 (MQ=255)
            gCCTGGGCTATTTCAACGACGGATGCTCTACTGGCACCGCTTTGCTCCTGAAAGCCGCATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAg                                                                                                                              >  2:532695/1‑143 (MQ=255)
            gCCTGGGCTATTTCAACGACGGATGCTCTACTGGCACCGCTTTGCTCCTGAAAGCCGCATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGACTGATAg                                                                                                                              >  1:551020/1‑143 (MQ=255)
                 ggCTATTTCAACGACGGATGCTCTACTGGCACCGCTTTGCTCCTGAAAGCCGCATGATGCGTAAAGTCACTGCTGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCCAATCATGCCCTATAGcacct                                                                                                                         >  1:294884/1‑141 (MQ=255)
                    tATTTCAACGACGTATGCTCTACTGGCACCGCTTTGCTCCTGAAAGCCGCATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATc                                                                                                                      <  1:917179/143‑1 (MQ=255)
                    tATTTCAACGACGGATGCTCTACTGGCACCGCTTTGCTCCTGAAAGCCGCATGGTGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCGGGATTAATCCATCAAATAATGCCTGATAGCACAtctc                                                                                                                      >  2:945773/1‑143 (MQ=255)
                      tttCAACGACGGATGCGCTACTGGCACCGCTTTGCTCCTGAAAGCCGCATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAg                                                                                                                    <  2:1012395/143‑1 (MQ=255)
                          aaCGACGGATGCTCTACTGGCACCGCTTTGCTCCTGAAAGCCGCATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGt                                                                                                                >  2:53358/1‑143 (MQ=255)
                             gaccGATGCTCTACTGCCACCGCTTTGCTCCTGAAAGCCGCATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGt                                                                                                             <  2:416916/143‑1 (MQ=255)
                              acgGATGCTCTACTGGCACCGCTTTGCTCCTGAAAGCCGCATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATACCATGCGTTgcc                                                                                                            >  2:344097/1‑141 (MQ=255)
                                  acggtttACTGGCACCGCTTTGCTCCTGGAAGCCGCATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCa                                                                                                        <  2:6689/137‑1 (MQ=255)
                                  aTGCTCTACTGGCACCGCTTTGCTCCTGAAAGCCGCATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCa                                                                                                        <  1:1064888/143‑1 (MQ=255)
                                      tctACTGGCACCGCTTTGCTCCTGAAAGCCGCATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCACCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCActtc                                                                                                    >  2:870876/1‑143 (MQ=255)
                                        tACTGGCACCGCTTTGCTCCTGAAAGCCGCATGATGCGTAAAGTCACTGATGCGTTCCTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCCTCACATAATGCCTGATAGCACCTACCAGGCGCGGTCCTCAcctctc                                                                                                  >  1:501711/1‑142 (MQ=255)
                                        tACTGGCACCGCTTTGCTCCTGAAAGCCGCATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCActtctt                                                                                                  >  1:999638/1‑143 (MQ=255)
                                          cTGGCACCGCTTTGCTCCTGAAAGCCGCATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCtttt                                                                                                >  2:377023/1‑143 (MQ=255)
                                          cTGGCACCGCTTTGCTCCTGAAAGCCGCATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCtttt                                                                                                <  1:675377/143‑1 (MQ=255)
                                          cTGGCACCGCTTTGCGCCTGAAAGCCGCATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGGCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCtttt                                                                                                <  1:675393/143‑1 (MQ=255)
                                                 cGCTTTGCTCCTGAAAGCCGCATGATGCGTAAAGTCACTGATACGTTACTAGATTATGGTCACAAAGTCCTTCGCCAGGATTTATCTATCAAATAATGGCTGAGAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTc                                                                                         >  1:145561/1‑143 (MQ=255)
                                                   cTTTGCTCCTGAAAGCCGCATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGGATTCCt                                                                                       <  1:696301/143‑1 (MQ=255)
                                                        cTCCTGAAAGCCGCATGATGCGTCAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTACTCCATCCAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGCat                                                                                  >  1:1045906/1‑143 (MQ=255)
                                                         tCCTGAAAGCCGCATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTATTTCTATTACTTGAatc                                                                                 >  1:857435/1‑143 (MQ=255)
                                                           cTGAAAGCCGCATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAatcac                                                                               >  1:824843/1‑143 (MQ=255)
                                                                  ccGCATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACaa                                                                        >  1:583332/1‑143 (MQ=255)
                                                                  ccGCATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACaa                                                                        >  2:597277/1‑143 (MQ=255)
                                                                    gCATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACaaaa                                                                      >  1:623712/1‑143 (MQ=255)
                                                                     cATGCTGCGTAACGCCACTGACGCGTTCCTCGATTATTGTCACAAAGTCCTTCGTCAGGATTAATCCATCCAACAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTACTCCTCGAATCACATCACAAAAt                                                                     <  2:632558/143‑1 (MQ=255)
                                                                     cATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCCGGATTAATCCATCACATAATGCCTGATAGCCCATATCGGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGGATCACACCACAAAAt                                                                     >  2:999295/1‑143 (MQ=255)
                                                                     cATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATACGGTGTTTTGCGCGGTTTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACAACAt                                                                     >  1:1054414/1‑143 (MQ=255)
                                                                      atgatgCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCACCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATa                                                                    <  1:870876/143‑1 (MQ=255)
                                                                       tgatgCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAg                                                                   >  2:634324/1‑143 (MQ=255)
                                                                           gCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACaa                                                               >  1:174207/1‑143 (MQ=255)
                                                                                 aGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACAAATCTCa                                                         <  1:755460/143‑1 (MQ=255)
                                                                                   tCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACAAATCTCaag                                                       >  1:146620/1‑141 (MQ=255)
                                                                                      cggaggCGTGACTCGATTATGGTCACAAACATCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACAAATCTCAggcgg                                                    <  2:564604/138‑1 (MQ=255)
                                                                                             gTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACAAATCTCAGGCGGCaaaaaa                                             <  1:843106/143‑1 (MQ=255)
                                                                                                  tCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACACATCTCAGGCGGCAAAAAccgacg                                        >  2:903915/1‑143 (MQ=255)
                                                                                                  tCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACAAATCTCAGGCGGCAAAAAacgacg                                        >  1:867999/1‑143 (MQ=255)
                                                                                                   cGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACAAATCTCAGGCGGCAAAAAACGACGt                                       <  2:342535/143‑1 (MQ=255)
                                                                                                         tGGTCACAAAGCCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAAt                                 <  1:368813/143‑1 (MQ=255)
                                                                                                            tCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACAAATCTCGGGCGGCAAAAAACGACGTCTGAATGCa                              >  1:53029/1‑143 (MQ=255)
                                                                                                            tCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCa                              <  1:579720/143‑1 (MQ=255)
                                                                                                            tCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCATGAACCACATCACAACATAGACAAATCTCAGGCGGCAAAAAACAAAGTCTGaaacca                              >  2:12154/1‑139 (MQ=255)
                                                                                                            tCACAAAGTCCTGCGCCAGGATTACTCCATCAACTAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTGTTTGTATTCCTTGAATCACATCACAAAATAGACAAATCTCAGGCGCCAAAAAACGACGTCTGAATGCa                              <  1:963432/143‑1 (MQ=255)
                                                                                                             cacaGAGTCCTGCGTCGAGATTAATCCACCAACTACTGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTCTTGTATTCCTTGAATCACATCACAAAATATACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCAt                             <  1:172032/143‑1 (MQ=255)
                                                                                                               caAAGTCCTTCGTCAGGATTAATCCACCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACAAAt                                                             >  1:719536/1‑109 (MQ=255)
                                                                                                               caAAGTCCTTCGTCAGGATTAATCCACCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACAAAt                                                             <  2:719536/109‑1 (MQ=255)
                                                                                                                aaaGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTGACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACAAATCTCAGGCGGCAAAAAAAGAAGTTTGTATGCAtttg                          >  1:388769/1‑142 (MQ=255)
                                                                                                                   gTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCAGTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCAttttttt                       <  1:532695/143‑1 (MQ=255)
                                                                                                                      cTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACAAATCTCAGGCTGCAAACAACGACGGCTGAATGCATTTTTTTTgg                    >  2:768629/1‑142 (MQ=255)
                                                                                                                      cTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAAACATATCTCAGGCGGCAGAAAACGACGTCTGAATGCATTTTTTTTGc                    >  1:1051866/1‑143 (MQ=255)
                                                                                                                       gtcataaGGATTGATCCATCGAATCATGAGTGATAACACAGATCAGGCAATGTACGAACTTCTATTTAGACACCATGAATCACATCACAAGATAGACAAATAAAAAAAAAAAAAAAACGACGTCTGAATGCATTTTTTTTGCt                   <  2:123810/137‑1 (MQ=255)
                                                                                                                          gacAGGAGTAATCCATCAAATCATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACATATCTCAGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTGCTGGc                <  1:43854/141‑1 (MQ=255)
                                                                                                                          gTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTGCTGGc                <  2:536563/143‑1 (MQ=255)
                                                                                                                          gTCAGAATTAATCCATCGAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATTACAAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTGCTGGc                <  2:29813/143‑1 (MQ=255)
                                                                                                                            aaGGATTGATCCCTCAAATAAGGCCGGATAGCACAGATCAGGCGTTGTCCTCACTTCTTATTGTAGTCCTTGAATCACATCACAAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTGCTGGCGa              <  2:298118/142‑1 (MQ=255)
                                                                                                                              ggATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTGCTGGCGACa            >  1:309682/1‑143 (MQ=255)
                                                                                                                                      ccATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTGCTGGCGACAAACCCACg    <  2:623712/143‑1 (MQ=255)
                                                                                                                                       cATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCGTCACTTCTTGTTGTATTCCTTGAATCACATCACAAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCAGTTTTTTTGCTGGCGACAAACCCACGt   >  1:581866/1‑143 (MQ=255)
                                                                                                                                       cATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTGCTGGCGACAAACCCACGt   >  1:565584/1‑143 (MQ=255)
                                                                                                                                       cATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTGCTGGCGACAAACCCACGt   >  1:584970/1‑143 (MQ=255)
                                                                                                                                       cATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTGTGTCTTCCTTGAATCACATCACAAAATAGACAAGTTTGAGGCGGCAATAAACGACGTCTGAACGAATTTTTCGTGCTGGAGAAAAACCCACGt   >  2:862917/1‑143 (MQ=255)
                                                                                                                                        aTCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATACCATCACAAAATAGACAAATCTCAGGCGGCAAACAACGACGTCTGAATGCATTTTTTTTGCTGTCGACAAAccctcata  >  2:1040479/1‑138 (MQ=255)
                                                                                                                                         |                                                                                                                                             
GATTGACTTAACGCCTGGGCTATTTCAACGACGGATGCTCTACTGGCACCGCTTTGCTCCTGAAAGCCGCATGATGCGTAAAGTCACTGATGCGTTACTCGATTATGGTCACAAAGTCCTTCGTCAGGATTAATCCATCAAATAATGCCTGATAGCACATATCAGGCGTTGTCCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTGCTGGCGACAAACCCACGTA  >  NC_000913/3060514‑3060792

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: