Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,224,845 T→C 5.2% C340R (TGC→CGC)  metH → homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,224,8450TC 5.2% 241.3 / 2.8 77C340R (TGC→CGC) metHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base T (36/37);  new base C (2/2);  total (38/39)
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.65e-01

ATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGT  >  NC_000913/4224698‑4224984
                                                                                                                                                   |                                                                                                                                           
aTGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCt                                                                                                                                             <  3:378096/148‑1 (MQ=255)
 tGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTg                                                                                                                                            <  7:214994/148‑1 (MQ=255)
       aaCATATACGTGAATGGTCGCAAGCGGGTTTTCTAAATATCGTCTGCGGATGCTGTGGCACCACGACACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCCGAAATTCCCGTAGCCTGCCGttt                                                                                                                                      >  6:169001/1‑148 (MQ=255)
       aaCAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGttt                                                                                                                                      >  3:314757/1‑148 (MQ=255)
               cGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGc                                                                                                                               <  3:355945/147‑1 (MQ=255)
                gTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGcc                                                                                                                              <  8:266528/147‑1 (MQ=255)
                gTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCt                                                                                                                             >  1:83444/1‑148 (MQ=255)
                gTGAATGGGCGCAAGCGGGGTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCt                                                                                                                             <  4:363543/148‑1 (MQ=255)
                 tGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTg                                                                                                                            >  6:161795/1‑148 (MQ=255)
                   aaTGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGc                                                                                                                               <  4:306868/143‑1 (MQ=255)
                   aaTGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGc                                                                                                                               >  3:306868/1‑143 (MQ=255)
                           cAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGa                                                                                                                  <  4:314757/148‑1 (MQ=255)
                             aGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAAc                                                                                                                <  5:161795/148‑1 (MQ=255)
                             aGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAAc                                                                                                                <  4:73488/148‑1 (MQ=255)
                                gggTTTTCTCAATATCCTCGGCGGCTGCTGTGTCACCACGCCACAACATATTTCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACAtt                                                                                                             <  1:388436/148‑1 (MQ=255)
                                 ggTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGGAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTg                                                                                                            <  5:377231/148‑1 (MQ=255)
                                 ggTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTg                                                                                                            >  1:423913/1‑148 (MQ=255)
                                         cAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCTGCCTGGAGCCGCTGAACATTGGCGAAGAt                                                                                                    >  6:163813/1‑148 (MQ=255)
                                          aaTATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGAt                                                                                                    <  6:175379/147‑1 (MQ=255)
                                           atatCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAg                                                                                                  <  2:423913/148‑1 (MQ=255)
                                           atatCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTGTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAg                                                                                                  <  5:177312/148‑1 (MQ=255)
                                           atatCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCATTTGTCCGGCCTGGAGCCg                                                                                                                      >  2:367395/1‑128 (MQ=255)
                                           atatCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCATTTGTCCGGCCTGGAGCCg                                                                                                                      <  1:367395/128‑1 (MQ=255)
                                               cgtcgGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTg                                                                                              >  7:442679/1‑148 (MQ=255)
                                                  cggcggCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGTGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAg                                                                                                      >  4:34814/1‑137 (MQ=255)
                                                  cggcggCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGTGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAg                                                                                                      <  3:34814/137‑1 (MQ=255)
                                                  cggcggCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCCGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGa                                                                                                     >  2:244411/1‑138 (MQ=255)
                                                  cggcggCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCCGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGa                                                                                                     >  2:251583/1‑138 (MQ=255)
                                                  cggcggCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCCGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGa                                                                                                     <  1:251583/138‑1 (MQ=255)
                                                  cggcggCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCCGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGa                                                                                                     <  1:244411/138‑1 (MQ=255)
                                                    gcggcTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTtgt                                                                                         >  4:382309/1‑148 (MQ=255)
                                                       gctgctGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGaa                                                                                      <  3:382309/148‑1 (MQ=255)
                                                       gctgctGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTACCGTAGACTGCCGGTTGTCCGTCCTGGAGCCGCTTAACATTGGAGAAGATCGCCTGTTTGTGaa                                                                                      >  8:335765/1‑148 (MQ=255)
                                                         tgctgTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACg                                                                                    >  2:224701/1‑148 (MQ=255)
                                                          gctgTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGt                                                                                   >  1:146201/1‑148 (MQ=255)
                                                           ctgTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACg                                                                                    <  8:279559/146‑1 (MQ=255)
                                                              tgGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTgg                                                                                 <  2:398103/146‑1 (MQ=255)
                                                               ggCACCACGCCACAACATATTGCAGCGATGAGTCGTGAAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATAGCCGAAGATAGCCTGTTTGTGAACGTGGGt                                                                               >  8:458034/1‑147 (MQ=255)
                                                                   ccacGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTg                                                                                              >  6:395505/1‑128 (MQ=255)
                                                                   ccacGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTg                                                                                              <  5:395505/128‑1 (MQ=255)
                                                                          acaacaTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGtt                                                                                                                                       >  1:192301/1‑80 (MQ=255)
                                                                          acaacaTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGtt                                                                                                                                       <  2:192301/80‑1 (MQ=255)
                                                                          acaacaTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGCACGTGGGTGAACGCACCAAc                                                                   >  6:103988/1‑148 (MQ=255)
                                                                          acaacaTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAAc                                                                   >  7:183793/1‑148 (MQ=255)
                                                                            aacaTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGt                                                                 >  5:333995/1‑148 (MQ=255)
                                                                              caTATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCa                                                               <  7:82341/148‑1 (MQ=255)
                                                                                tatTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGc                                                                                                 >  5:169674/1‑112 (MQ=255)
                                                                                tatTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGc                                                                                                 <  6:169674/112‑1 (MQ=255)
                                                                                 atTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCg                                                            >  8:133594/1‑148 (MQ=255)
                                                                                       gCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCg                                                      >  1:449454/1‑148 (MQ=255)
                                                                                         gATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGc                                                     >  6:42516/1‑147 (MQ=255)
                                                                                           tGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCt                                                    >  7:512/1‑146 (MQ=255)
                                                                                                cGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGtt                                                         <  2:39236/136‑1 (MQ=255)
                                                                                                cGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGtt                                                         >  1:39236/1‑136 (MQ=255)
                                                                                                   gCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAgc                                          >  2:400307/1‑148 (MQ=255)
                                                                                                            ggATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCa                                 <  3:368783/148‑1 (MQ=255)
                                                                                                             gATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCa                                 >  1:179984/1‑147 (MQ=255)
                                                                                                              aTTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCa                                 >  7:32583/1‑146 (MQ=255)
                                                                                                               ttAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCa                                 >  2:351073/1‑145 (MQ=255)
                                                                                                               ttAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAaag                              <  2:449454/148‑1 (MQ=255)
                                                                                                                tAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAaag                              >  7:100296/1‑147 (MQ=255)
                                                                                                                    ggcgcgCATTCTGCCGGTAATTCCCTTAGCCTGCCGTTTGTCCGGCGTTGAGCCTGTGAATATTGTCTTAGATATCCTGTTTGTGAATGTGGGTTTACGCACCAATGTCAc                                                              <  5:297221/109‑1 (MQ=255)
                                                                                                                    gccgcGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAaagaa                            >  3:393823/1‑145 (MQ=255)
                                                                                                                     ccgcGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAagag                         >  6:398210/1‑147 (MQ=255)
                                                                                                                      cgcgCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGaa                       <  5:247687/148‑1 (MQ=255)
                                                                                                                      cgcgCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGaa                       <  5:146760/148‑1 (MQ=255)
                                                                                                                        cgcAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAAt                     <  2:146201/148‑1 (MQ=255)
                                                                                                                            aaCTGACGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAg                 <  5:42516/148‑1 (MQ=255)
                                                                                                                               tGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGc                >  5:38993/1‑146 (MQ=255)
                                                                                                                                 ccGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGa              <  6:396639/146‑1 (MQ=255)
                                                                                                                                 ccGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGa              <  8:442679/146‑1 (MQ=255)
                                                                                                                                  cGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGgc           >  1:283759/1‑148 (MQ=255)
                                                                                                                                   ggAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGgcg          >  4:112334/1‑148 (MQ=255)
                                                                                                                                    gAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGgcgc         <  6:333995/148‑1 (MQ=255)
                                                                                                                                     aaaTTCCCGTAGCCTGCCGTTGGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCt        <  8:141577/148‑1 (MQ=255)
                                                                                                                                       aTTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTgg      <  2:83444/148‑1 (MQ=255)
                                                                                                                                       aTTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTgg      <  8:32583/148‑1 (MQ=255)
                                                                                                                                           ccGGAGCCTGCCGTTTGGCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTCTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGt  <  5:163813/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
ATGGCAAAACAGATACGTGAATGGGCGCAAGCGGGTTTTCTCAATATCGTCGGCGGCTGCTGTGGCACCACGCCACAACATATTGCAGCGATGAGTCGTGCAGTAGAAGGATTAGCGCCGCGCAAACTGCCGGAAATTCCCGTAGCCTGCCGTTTGTCCGGCCTGGAGCCGCTGAACATTGGCGAAGATAGCCTGTTTGTGAACGTGGGTGAACGCACCAACGTCACCGGTTCCGCTAAGTTCAAGCGCCTGATCAAAGAAGAGAAATACAGCGAGGCGCTGGATGT  >  NC_000913/4224698‑4224984

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