Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,061,818 G→C 5.7% W323C (TGG→TGC scpA → methylmalonyl‑CoA mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,061,8180GC 5.7% 185.6 / 2.1 74W323C (TGG→TGCscpAmethylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base G (36/32);  new base C (3/2);  total (40/34)
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.90e-01

CTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGC  >  NC_000913/3061677‑3061959
                                                                                                                                             |                                                                                                                                             
cttcGGCATCGGCATGGATCTGTGTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGa                                                                                                                                              >  2:704677/1‑143 (MQ=255)
 ttcGACATCGGCATCGATCTGTTTATGACCCTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAACCGGACAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAg                                                                                                                                             <  2:1107920/143‑1 (MQ=255)
  gcggcatcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGc                                                                                                                                            <  1:503072/142‑1 (MQ=255)
    ggcatcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGGGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCct                                                                                                                                          <  2:407372/143‑1 (MQ=255)
      cgtcggcatGGAGCTGTTTATGGACGTCGCCGTGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCctga                                                                                                                                        <  1:474097/141‑1 (MQ=255)
      catcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGACGCGGTCAGTGGATTTGGCGCCCAGGACCCGAAATCACTGGCGCTGCCTACCCCCTGCCAGACCTCAGGCTGGAGCctgc                                                                                                                                        >  2:244170/1‑142 (MQ=255)
        tcggcatGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCctgact                                                                                                                                      <  1:472713/143‑1 (MQ=255)
          ggcatGCATCTGTTTATGACCGTCGCCCTGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTCCAGCctgaggtc                                                                                                                                    <  2:770523/143‑5 (MQ=255)
           gcatGGCTCTGTTTATGAGCGTCGCCATGTTGCGTGCGGCACGTTATGTAGGGAGCGAAGCGGTCAGTGGATTCGGCGCACAGGACCCGAAAGCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGaa                                                                                                                                   <  2:35108/143‑1 (MQ=255)
             gtGGATCTGTTTACGAACGTCGCCATGTTGCGTGCGGCACGTTATCTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACa                                                                                                                                 <  2:648238/142‑1 (MQ=255)
             atGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTCGCGCAGCGAAACCACTGCCACACCTCAGGCTCGAGCCCGACCGcaca                                                                                                                                 >  2:363380/1‑143 (MQ=255)
               ggATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGACCAgg                                                                                                                               >  1:604822/1‑143 (MQ=255)
                    tGTTTATGAGCGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCg                                                                                                                          <  1:598413/143‑1 (MQ=255)
                    tGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCg                                                                                                                          >  1:994091/1‑143 (MQ=255)
                     gTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGt                                                                                                                         <  2:437688/143‑1 (MQ=255)
                     gTTGATGAACGTCGCGCTGTTGCGTGCGGCACGTTATTTAGGGAGCGAAGCGGTCAGTGGACTTGGCGCACAGGACCCGAAATCACTGGCGCCGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGt                                                                                                                         <  2:991980/143‑1 (MQ=255)
                         aTGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaa                                                                                                                     >  2:362640/1‑143 (MQ=255)
                         aTGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaa                                                                                                                     >  1:94121/1‑143 (MQ=255)
                           gAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCCCTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaaca                                                                                                                   >  2:893455/1‑143 (MQ=255)
                           gAACGTCGCCATGTAGCGTGCGACACGTTATGTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGCACCCGAAATCACTGGCGCTTCGTGCCCACTGCCAGACCGCACGCTGCAGCCTGCCTGAACAAGCTCCggcacaca                                                                                                                   >  2:1085737/1‑136 (MQ=255)
                             aCGTCTCCCTGTTGCGTGGCGCACTGTATTTATAGTGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACATACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaacaac                                                                                                                 <  2:394676/143‑1 (MQ=255)
                             aCGTCGCCGTGTGGCGTGCGGCACGTTATTGATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATaacaac                                                                                                                 <  2:813541/143‑1 (MQ=255)
                                 cGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTa                                                                                                             <  1:908607/143‑1 (MQ=255)
                                  gCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAGCAGGATCCGTATGACAACGTCAt                                                                                                            >  1:160764/1‑143 (MQ=255)
                                            gtgcgGCACGTTATTTATGGAGCGAAGCAGTCAGTGGACTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGACACGTATAACGACGTCATACGAcccacc                                                                                                  >  2:120361/1‑143 (MQ=255)
                                             tgcgGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGTCTGGAGCCTGCCTGCACAGGATCCGTATACCCACGTTAGCCGCAccacca                                                                                                 >  2:347522/1‑143 (MQ=255)
                                              gcgGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGGCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCAt                                                                                                >  1:135261/1‑143 (MQ=255)
                                               cgGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCAtt                                                                                               >  1:385188/1‑143 (MQ=255)
                                                  cACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACGCGAAATCACAAGCGCTGCATACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGaa                                                                                            <  1:509841/143‑1 (MQ=255)
                                                  cACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCCGACCGAACACGATCCGTATAACAACGTTATCCGCCCCACCATTCaa                                                                                            >  2:286916/1‑143 (MQ=255)
                                                  cACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTCGCGCCTGACTGAACAGGATCCGTATAACAACGTCATCCGCaccaccctcgca                                                                                            >  2:42987/1‑137 (MQ=255)
                                                     gTTCTTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGCACCCACTGCCAGACCTCAGCCTGGAGCCTGACTGAACGGCATCCGTCTAACCACGTTCTCCGCCCCCCCATTCAAgcg                                                                                         >  1:508260/1‑143 (MQ=255)
                                                     gTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAgcg                                                                                         >  1:214761/1‑143 (MQ=255)
                                                       tacttatGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCgct                                                                                       <  2:341863/140‑1 (MQ=255)
                                                               gggCGAAGCGGCCGGTGGATCTGCCGCACAGGACCCGAAAGCACTGGCGCTGCGTACCCACTGCCAGACCCCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGa                                                                               <  2:860508/141‑1 (MQ=255)
                                                               gAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGa                                                                               >  1:679012/1‑143 (MQ=255)
                                                                        ggTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACAGCCAGACCTCAGCCTGGAGCCTGACCGAACAGGATCACTAGAACACCGCTACCCGCACGACCATGGCAGCGCAGGCTGCGACGCTGGGcc                                                                      >  1:254868/1‑142 (MQ=255)
                                                                         gTCAGTGGATTTGGCGCACAGGACCCGAAATCACTTGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCAGTGCAGCGCTGGCTGCGACGCTGGGCgg                                                                     >  1:259096/1‑143 (MQ=255)
                                                                                 acttGGCGCACAGGACCCGAAATCACCGGCGCTGCGTACCCACTGCCAGACCGCAGGCTGGAGCCTGACTGAACAGGCTCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGt                                                             <  1:563372/141‑1 (MQ=255)
                                                                                     ggCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCAGCAGCAGTGATGCGCTGGCTGCGAGGCTGGGCGGTAGTCAGacacg                                                         >  1:203977/1‑142 (MQ=255)
                                                                                      gcgcACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGACGGACTGAACAGGATCCGTATAACAACGTTATCCGCACCCCCATTGAAGCGCTGGCTGCGACGCTGGGCGCTACCCAGCCACTg                                                        >  2:385123/1‑143 (MQ=255)
                                                                                        gcACTGTACCCTAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCa                                                      <  2:102437/143‑1 (MQ=255)
                                                                                        gcACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCCCTGCa                                                      >  1:122506/1‑143 (MQ=255)
                                                                                        gcACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTgta                                                      >  1:16309/1‑141 (MQ=255)
                                                                                          acaGGACCCGAAATCACTGACTTGATTTTCTCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGCTCCGTATAACAACGTTATCCGCACCACCACTGAAGCGCTGGCTGCGAAGCCGGACAGTACTCAGTCACTGCGTa                                                    >  2:137968/1‑143 (MQ=255)
                                                                                               aCCCGACATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAAc                                               <  1:882870/143‑1 (MQ=255)
                                                                                                 ccGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTGTAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCCGTCACTGCATACCACCGc                                             >  1:376938/1‑143 (MQ=255)
                                                                                                 ccGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTCTCCGCACCCCCATTGACGCGCTCGCTGCGCCGCTGGGCGGTACTCACTCACTGCGTACCCAccc                                             >  2:838966/1‑141 (MQ=255)
                                                                                                 ccGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGACTGAAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACAATTAAAGCGCGGTCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCGACGc                                             >  1:691722/1‑143 (MQ=255)
                                                                                                  cGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTCTCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGcc                                            <  1:42478/143‑1 (MQ=255)
                                                                                                   gAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCt                                           <  1:115597/143‑1 (MQ=255)
                                                                                                     aagcACTGGCGCGGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCttt                                         <  1:88493/140‑1 (MQ=255)
                                                                                                     aaTCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCttt                                         <  1:921999/143‑1 (MQ=255)
                                                                                                      aTCGCTGGCGCTGCGGACCCACTGCCAGACCGCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGGTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTg                                        <  2:355892/143‑1 (MQ=255)
                                                                                                      aTCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCGGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTg                                        <  1:959945/143‑1 (MQ=255)
                                                                                                                 tGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGACACTGCATACCAACGCCTATGACGCAGCGCtt                             >  2:208422/1‑143 (MQ=255)
                                                                                                                  gCGTACGCACTGCCAGACCTCAGGCGGCAGGCTCACCGAACAGGATGCGTATAACGACGTTATCCGCACGACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTg                            <  2:984821/143‑1 (MQ=255)
                                                                                                                  gCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTg                            <  1:605251/143‑1 (MQ=255)
                                                                                                                   cGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTgg                           <  2:503007/143‑1 (MQ=255)
                                                                                                                     tACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCAGTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGCCGAAGCGCTTgggt                         >  2:735626/1‑141 (MQ=255)
                                                                                                                      ccccACGGCTAGGCCTCAGGCTGGAGCCGGACTGCACAGGATCGGTGTGACAGCGGTATCCGCACCACCAGTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGGCGAAGCGCTTGGttt                        <  2:932525/142‑1 (MQ=255)
                                                                                                                        ccACTGCCAGACCTCAGGCTGGCGCCTGACTGAACAGGATCCGTATCCCACCGTTATCCGCACCACCCGTCAAGCGCTGGCCGCGCCGCGGGGCGGTACTCAGTCACTGCATACCCCCCCCTTTGCCGAAGCGCGTGGTTTGc                      >  1:543719/1‑143 (MQ=255)
                                                                                                                        ccACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCGGCACCACAATTGAAGAGGTGGCTTCGACAATCGTTGTTTCTTCGGCACTGCATACCAACGCCCTTGGCGAAGCGCATGGtctcc                      >  2:1075579/1‑139 (MQ=255)
                                                                                                                        ccACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATCCCCACGCCTTTGACGAAGCGCTTGGTTTGc                      >  1:267667/1‑143 (MQ=255)
                                                                                                                        ccACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGc                      <  1:893455/143‑1 (MQ=255)
                                                                                                                        ccACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGc                      <  2:12461/143‑1 (MQ=255)
                                                                                                                         cACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCCCCATTGAAGCGCTGGCCGCGACGCTGGGCGGTACACCGTCACTGCATACCAACGCCTTTGCCGACGCGCTTggcctccc                     >  2:230663/1‑137 (MQ=255)
                                                                                                                          aCTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTAATCAGTCACGTCATGACAACGCCGCTGAGGTAGCCCCTGGGTTGCCt                    >  2:941541/1‑143 (MQ=255)
                                                                                                                              tgttcagtCAGGCTCCAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCg                >  1:537666/8‑143 (MQ=255)
                                                                                                                               cAGACCTCAGGCCGGAGCCCGAGTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTTACGAAGCGCTTGGTTTGTCTACCGa               <  1:216186/143‑1 (MQ=255)
                                                                                                                                aGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGAt              >  1:1070924/1‑143 (MQ=255)
                                                                                                                                   agtCAGGCTCCAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTCCCTACCGACTtc           >  2:101371/3‑143 (MQ=255)
                                                                                                                                       aGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAg       >  1:941855/1‑143 (MQ=255)
                                                                                                                                       aGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAg       <  1:227877/143‑1 (MQ=255)
                                                                                                                                            ggAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGc  >  1:1021274/1‑143 (MQ=255)
                                                                                                                                            ggAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTACATACCAACGCCTTTGAAGAAGCGCTTGGATTGCCTACCGATTACTCAGCACGc  >  2:531223/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                             
CTTCGGCATCGGCATGGATCTGTTTATGAACGTCGCCATGTTGCGTGCGGCACGTTATTTATGGAGCGAAGCGGTCAGTGGATTTGGCGCACAGGACCCGAAATCACTGGCGCTGCGTACCCACTGCCAGACCTCAGGCTGGAGCCTGACTGAACAGGATCCGTATAACAACGTTATCCGCACCACCATTGAAGCGCTGGCTGCGACGCTGGGCGGTACTCAGTCACTGCATACCAACGCCTTTGACGAAGCGCTTGGTTTGCCTACCGATTTCTCAGCACGC  >  NC_000913/3061677‑3061959

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