Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,324,650 T→C 6.0% M181T (ATG→ACG)  atoA → acetyl‑CoA:acetoacetyl‑CoA transferase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,324,6500TC 6.0% 218.8 / 3.3 67M181T (ATG→ACG) atoAacetyl‑CoA:acetoacetyl‑CoA transferase, beta subunit
Reads supporting (aligned to +/- strand):  ref base T (28/35);  new base C (2/2);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.59e-01

CGCAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTATGATTGGTCGCATATCGCGTTTTATGACGCGTTTTGTCAGCC  >  NC_000913/2324505‑2324798
                                                                                                                                                 |                                                                                                                                                    
cgcAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAAtgtg                                                                                                                                                   >  2:131154/1‑149 (MQ=255)
 gcAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTgg                                                                                                                                                  <  1:581139/149‑1 (MQ=255)
   aaaGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCt                                                                                                                                                <  1:477684/149‑1 (MQ=255)
           catcGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGaaa                                                                                                                                        <  3:389254/149‑1 (MQ=255)
                ccATGGAAGATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGcc                                                                                                                                   <  2:288179/149‑1 (MQ=255)
                     gAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACgg                                                                                                                              <  1:21179/149‑1 (MQ=255)
                       aCATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGgt                                                                                                                            <  3:418419/149‑1 (MQ=255)
                       aCATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAACGTGGCTCACCGAAATTGCCGACGGgt                                                                                                                            <  1:157140/149‑1 (MQ=255)
                         aTTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCACTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGgtgt                                                                                                                          >  4:586506/1‑149 (MQ=255)
                              gcCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGAttt                                                                                                                     >  3:7331/1‑149 (MQ=255)
                                       ggTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACcgt                                                                                                            <  2:75274/149‑1 (MQ=255)
                                        gTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAACGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACcgtg                                                                                                           >  1:71690/1‑149 (MQ=255)
                                         ttCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGTAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACcgtgc                                                                                                          >  1:121048/1‑149 (MQ=255)
                                             gCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCgtgcgtgc                                                                                                      >  4:77792/1‑149 (MQ=255)
                                              cAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGcc                                                                                                     <  4:150433/149‑1 (MQ=255)
                                                  aaTTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCaaaa                                                                                                 <  1:131154/149‑1 (MQ=255)
                                                  aaTTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCaaaa                                                                                                 <  3:273819/149‑1 (MQ=255)
                                                     tttGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTACCAAAACAg                                                                                              <  4:311782/149‑1 (MQ=255)
                                                          gccgcTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCt                                                                                         >  3:14053/1‑149 (MQ=255)
                                                                   ccATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAACGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGaa                                                                                <  2:71690/149‑1 (MQ=255)
                                                                    cATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTa                                                                                                                    >  1:437847/1‑112 (MQ=255)
                                                                    cATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTa                                                                                                                    <  2:437847/112‑1 (MQ=255)
                                                                           ctcactGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTcgccgc                                                                        >  3:515471/1‑149 (MQ=255)
                                                                                tGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATc                                                                   <  4:421906/149‑1 (MQ=255)
                                                                                tGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATc                                                                   >  2:417296/1‑149 (MQ=255)
                                                                                 gcgcAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCt                                                                  <  3:459314/149‑1 (MQ=255)
                                                                                  cgcAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTg                                                                 >  3:77312/1‑149 (MQ=255)
                                                                                      aCATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATa                                                             <  1:544232/149‑1 (MQ=255)
                                                                                        aTGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACa                                                                                               <  2:320994/113‑1 (MQ=255)
                                                                                        aTGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACa                                                                                               >  1:320994/1‑113 (MQ=255)
                                                                                            ggTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAAc                                                       >  3:434227/1‑149 (MQ=255)
                                                                                                 atatGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGt                                                  >  1:147007/1‑149 (MQ=255)
                                                                                                 atatGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGt                                                  <  1:99074/149‑1 (MQ=255)
                                                                                                  tatGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTg                                                 >  2:144186/1‑149 (MQ=255)
                                                                                                   atGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGa                                                <  4:14053/149‑1 (MQ=255)
                                                                                                       tGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTa                                            <  4:570247/149‑1 (MQ=255)
                                                                                                       tGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTa                                            <  3:606330/149‑1 (MQ=255)
                                                                                                         gTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTATg                                          >  1:105448/1‑149 (MQ=255)
                                                                                                           tACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTATGAt                                        >  3:119458/1‑149 (MQ=255)
                                                                                                            aCTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTATGAtt                                       <  1:363802/149‑1 (MQ=255)
                                                                                                                aaCTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTATGATTGGTc                                   >  1:190151/1‑149 (MQ=255)
                                                                                                                     gCTGTCTTTCGTTTTATTGACGGCAAAACGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTATGATTGGTCGCata                              >  2:268207/1‑149 (MQ=255)
                                                                                                                        gTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGAtt                                              >  1:467578/1‑130 (MQ=255)
                                                                                                                        gTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGAtt                                              <  2:467578/130‑1 (MQ=255)
                                                                                                                        gTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTATGATTGGTCGCATATcg                           >  1:232587/1‑149 (MQ=255)
                                                                                                                           tttCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGaa                                                                                            >  3:467936/1‑81 (MQ=255)
                                                                                                                           tttCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGaa                                                                                            <  4:467936/81‑1 (MQ=255)
                                                                                                                           tttCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTATGATTGGTCGCATATCGCGt                        >  4:159520/1‑149 (MQ=255)
                                                                                                                            ttCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGTCACCGTGCGTGCCTAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACggg                                                    >  3:496724/1‑120 (MQ=255)
                                                                                                                            ttCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACggg                                                    <  4:496724/120‑1 (MQ=255)
                                                                                                                               gTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTATGATTGGTCGCATATCGCGTTTTa                    <  2:105448/149‑1 (MQ=255)
                                                                                                                                  ttATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTATGATTGGTCGCATATCGCGTTTTATGa                 <  2:121048/149‑1 (MQ=255)
                                                                                                                                   tATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACTGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTATGATTGGTCGCATATCGCGTTTTATGAc                <  3:535637/149‑1 (MQ=255)
                                                                                                                                   tATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTATGATTGGTCGCATATCGCGTTTTATGAc                <  4:77312/149‑1 (MQ=255)
                                                                                                                                   tATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTATGATTGGTCGCATATCGCGTTTTATGAc                <  2:145714/149‑1 (MQ=255)
                                                                                                                                    aTTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTATGATTGGTCGCATATCGCGTTTTATGAcg               >  1:36027/1‑149 (MQ=255)
                                                                                                                                     ttGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTATGATTGGTCGCATATCGCGTTTTATGAcgc              <  2:301392/149‑1 (MQ=255)
                                                                                                                                         cGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTATGATTGGTCGCATATCGCGTTTTATGACGCGttt          <  2:36027/149‑1 (MQ=255)
                                                                                                                                          ggCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTATGATTGGTCGCATATCGCGTTTTATGACGCGtttt         >  1:85319/1‑149 (MQ=255)
                                                                                                                                          ggCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTATGATTGGTCGCATATCGCGTTTTATGACGCGtttt         >  1:604524/1‑149 (MQ=255)
                                                                                                                                           gCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACggg                                                    >  4:268062/1‑105 (MQ=255)
                                                                                                                                           gCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACggg                                                    <  3:268062/105‑1 (MQ=255)
                                                                                                                                           gCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTATGAtt                                       <  1:87343/118‑1 (MQ=255)
                                                                                                                                           gCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTATGAtt                                       >  2:497315/1‑118 (MQ=255)
                                                                                                                                           gCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTATGAtt                                       <  1:497315/118‑1 (MQ=255)
                                                                                                                                           gCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGGGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTATGAtt                                       >  2:87343/1‑118 (MQ=255)
                                                                                                                                            cAAAATGTGGCTCAACGAAAATGCCGACGCGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTATGATTGGTCGCATATCGCGTTTTGTGACGCGtgttgt       <  4:515785/149‑1 (MQ=255)
                                                                                                                                              aaaTGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTATGATTGGTCGCATATCGCGTTTTATGACGCGTTTTGTCa     >  3:427548/1‑149 (MQ=255)
                                                                                                                                                 tgtgGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTATGATTGGTCGCATATCGCGTTTTATGACGCGTTTTGTCAGcc  >  1:98691/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CGCAAAGTGATCATCGCCATGGAACATTGCGCCAAAGATGGTTCAGCAAAAATTTTGCGCCGCTGCACCATGCCACTCACTGCGCAACATGCGGTGCATATGCTGGTTACTGAACTGGCTGTCTTTCGTTTTATTGACGGCAAAATGTGGCTCACCGAAATTGCCGACGGGTGTGATTTAGCCACCGTGCGTGCCAAAACAGAAGCTCGGTTTGAAGTCGCCGCCGATCTGAATACGCAACGGGGTGATTTATGATTGGTCGCATATCGCGTTTTATGACGCGTTTTGTCAGCC  >  NC_000913/2324505‑2324798

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