Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,061,041 T→G 5.1% R64R (CGT→CGG scpA → methylmalonyl‑CoA mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,061,0410TG 5.1% 287.2 / 6.9 79R64R (CGT→CGGscpAmethylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base T (36/39);  new base G (2/2);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

GAGCCGTCGGGAGAAAACTGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTG  >  NC_000913/3060894‑3061186
                                                                                                                                                   |                                                                                                                                                 
gAGCCGTCGGGAGAAAACTGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGt                                                                                                                                                   <  2:90319/148‑1 (MQ=255)
 aGCCGTCGGGAGAAAACTGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTg                                                                                                                                                  <  4:650602/148‑1 (MQ=255)
       cGGGAGAAAACTGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTa                                                                                                                                            <  2:629941/148‑1 (MQ=255)
          gagaAAACTGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGt                                                                                                                                         >  3:205153/1‑148 (MQ=255)
                    tCGACTCGCTGGTTCATCAAACCGCGGAAGGTATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCa                                                                                                                               >  4:727240/1‑148 (MQ=255)
                        cTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGt                                                                                                                                         <  2:664262/134‑1 (MQ=255)
                        cTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGt                                                                                                                                         >  1:664262/1‑134 (MQ=255)
                          cGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATg                                                                                                                                          >  2:466396/1‑131 (MQ=255)
                          cGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATg                                                                                                                                          <  1:466396/131‑1 (MQ=255)
                           gCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCCTGCCACTATGTATACCGCCCAACCGTgg                                                                                                                        <  3:727240/148‑1 (MQ=255)
                                    tcaAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTa                                                                                                                                            <  4:1026962/119‑1 (MQ=255)
                                    tcaAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTa                                                                                                                                            >  3:1026962/1‑119 (MQ=255)
                                    tcaAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGt                                                                                                               >  2:145323/1‑148 (MQ=255)
                                    tcaAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACAATCCGt                                                                                                               >  3:746230/1‑148 (MQ=255)
                                         ccGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTa                                                                                                          >  4:631889/1‑148 (MQ=255)
                                              gAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGt                                                                                                                                         <  4:859143/112‑1 (MQ=255)
                                              gAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGt                                                                                                                                         >  3:859143/1‑112 (MQ=255)
                                               aaGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATAcc                                                                                                                                    >  1:1045432/1‑116 (MQ=255)
                                               aaGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATAcc                                                                                                                                    <  2:1045432/116‑1 (MQ=255)
                                               aaGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTc                                                                                                              >  1:428208/1‑138 (MQ=255)
                                               aaGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTc                                                                                                              <  2:428208/138‑1 (MQ=255)
                                                aGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGt                                                                                                   >  4:426312/1‑148 (MQ=255)
                                                   gATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGtat                                                                                                                                       <  3:389408/109‑1 (MQ=255)
                                                   gATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGtat                                                                                                                                       >  4:389408/1‑109 (MQ=255)
                                                   gATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGt                                                                                                                          >  3:74633/1‑122 (MQ=255)
                                                   gATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGt                                                                                                                          <  4:74633/122‑1 (MQ=255)
                                                   gATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGtttt                                                                                                <  3:631889/148‑1 (MQ=255)
                                                    aTCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGttttt                                                                                               <  4:205153/148‑1 (MQ=255)
                                                     tCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGt                                                                                                   >  4:1023134/1‑143 (MQ=255)
                                                     tCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGt                                                                                                   <  3:1023134/143‑1 (MQ=255)
                                                         cATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACa                                                                                          >  1:1126708/1‑148 (MQ=255)
                                                          aTCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAg                                                                                         >  2:860503/1‑148 (MQ=255)
                                                           tCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGt                                                                                                                                         >  1:12510/1‑99 (MQ=255)
                                                           tCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGt                                                                                                                                         <  2:12510/99‑1 (MQ=255)
                                                           tCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGc                                                                                        >  2:1000610/1‑148 (MQ=255)
                                                             aaGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCaa                                                                                      <  1:145323/148‑1 (MQ=255)
                                                              aGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCg                                                                                                                                   <  3:865566/102‑1 (MQ=255)
                                                              aGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCg                                                                                                                                   >  4:865566/1‑102 (MQ=255)
                                                              aGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCaaa                                                                                     <  2:1126708/148‑1 (MQ=255)
                                                               gCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCaaaa                                                                                    <  1:1067238/148‑1 (MQ=255)
                                                                  gCTGTATACCGAAGCCGATCTCGATAATCTTGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAagag                                                                                 <  4:746230/148‑1 (MQ=255)
                                                                  gCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAagag                                                                                 <  1:1000610/148‑1 (MQ=255)
                                                                   cTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTAt                                                                                                                                           >  1:1072094/1‑89 (MQ=255)
                                                                   cTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTAt                                                                                                                                           <  2:1072094/89‑1 (MQ=255)
                                                                   cTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGt                                                                                <  1:860503/148‑1 (MQ=255)
                                                                          ccGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTc                                                                                                              <  2:444652/111‑1 (MQ=255)
                                                                          ccGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTc                                                                                                              >  1:444652/1‑111 (MQ=255)
                                                                                    tctcGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGtttt                                                                                                <  3:348078/115‑1 (MQ=255)
                                                                                    tctcGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGtttt                                                                                                >  4:348078/1‑115 (MQ=255)
                                                                                    tctcGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGc                                                               >  4:858960/1‑148 (MQ=255)
                                                                                       cGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGt                                                                                <  2:301494/128‑1 (MQ=255)
                                                                                       cGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGt                                                                                >  1:301494/1‑128 (MQ=255)
                                                                                         aTAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGt                                                                                                                                                   <  3:588180/59‑1 (MQ=255)
                                                                                         aTAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGt                                                                                                                                                   >  4:588180/1‑59 (MQ=255)
                                                                                           aaTCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAAcc                                                        <  3:858960/148‑1 (MQ=255)
                                                                                                        cAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAAc                                                                                           >  4:1018341/1‑100 (MQ=255)
                                                                                                        cAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAAc                                                                                           <  3:1018341/100‑1 (MQ=255)
                                                                                                            tACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGtttt                                                                                                <  3:372089/91‑1 (MQ=255)
                                                                                                            tACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGtttt                                                                                                >  4:372089/1‑91 (MQ=255)
                                                                                                            tACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTc                                                                               >  2:1010791/1‑108 (MQ=255)
                                                                                                            tACCCTTCCTGGTATGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCAAACCGTGAACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTc                                                                               <  1:1010791/108‑1 (MQ=255)
                                                                                                               ccTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGt                                                                                <  1:692802/104‑1 (MQ=255)
                                                                                                               ccTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGt                                                                                >  2:692802/1‑104 (MQ=255)
                                                                                                                  tCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGTACCATCCGTCAGTATGCTGGTATTTCAACAGCAAAAGAATCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCtt                                 >  4:325613/1‑148 (MQ=255)
                                                                                                                  tCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTc                                   <  4:1049671/146‑1 (MQ=255)
                                                                                                                  tCCTGGTCTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTc                                   >  3:1049671/1‑146 (MQ=255)
                                                                                                                      ggTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCaa                                                                                      <  3:842344/91‑1 (MQ=255)
                                                                                                                      ggTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCaa                                                                                      >  4:842344/1‑91 (MQ=255)
                                                                                                                      ggTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGt                                                                                >  4:183436/1‑97 (MQ=255)
                                                                                                                      ggTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGt                                                                                <  3:183436/97‑1 (MQ=255)
                                                                                                                      ggTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTAt                                                                  <  4:279172/111‑1 (MQ=255)
                                                                                                                      ggTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTAt                                                                  >  3:279172/1‑111 (MQ=255)
                                                                                                                        tttGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAAcc                                                        <  3:309974/119‑1 (MQ=255)
                                                                                                                        tttGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAAcc                                                        >  4:309974/1‑119 (MQ=255)
                                                                                                                         ttGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGc                                                               >  4:411904/1‑111 (MQ=255)
                                                                                                                         ttGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGc                                                               <  3:411904/111‑1 (MQ=255)
                                                                                                                           gccgccCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCg                        <  3:325613/148‑1 (MQ=255)
                                                                                                                                        agTGGCACGCGGGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGt                                                                                <  2:804595/78‑1 (MQ=255)
                                                                                                                                        agTGGCACGCGGGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGt                                                                                >  1:804595/2‑79 (MQ=255)
                                                                                                                                        agTGGCACGCGGGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGgccgcc                                                >  2:165215/2‑111 (MQ=255)
                                                                                                                                        agTGGCACGCGGGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGgccgcc                                                <  1:165215/110‑1 (MQ=255)
                                                                                                                                                 cgTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTg  <  3:426312/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GAGCCGTCGGGAGAAAACTGTCGACTCGCTGGTTCATCAAACCGCGGAAGGGATCGCCATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTG  >  NC_000913/3060894‑3061186

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