Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,868,951 A→G 5.3% I190T (ATT→ACT)  pcm ← L‑isoaspartate protein carboxylmethyltransferase type II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,868,9510AG 5.3% 229.9 / 3.2 75I190T (ATT→ACT) pcmL‑isoaspartate protein carboxylmethyltransferase type II
Reads supporting (aligned to +/- strand):  ref base A (32/39);  new base G (2/2);  total (34/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 = 8.61e-01

GCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAATTTCCGGCGGTGCCGCCGTAACAATGATAGCGTCAAAC  >  NC_000913/2868806‑2869097
                                                                                                                                                 |                                                                                                                                                  
gCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGataat                                                                                                                                                   >  3:19641/1‑147 (MQ=255)
gCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGataat                                                                                                                                                   >  7:48350/1‑147 (MQ=255)
          gAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCTCCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTc                                                                                                                                        >  5:15920/1‑148 (MQ=255)
          gAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCATGGCCTCCACGGTATCGATAATAAATTCGCCTc                                                                                                                                        >  8:134620/1‑148 (MQ=255)
                tGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCcgac                                                                                                                                  >  7:297100/1‑148 (MQ=255)
                 ggAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAACCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACCGTATCGATAATAAATTCGCCTCCCcgac                                                                                                                                  >  8:403270/1‑147 (MQ=255)
                  gAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCcgac                                                                                                                                  <  3:58313/146‑1 (MQ=255)
                   aaaaGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAATGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCgacgac                                                                                                                               >  2:176998/1‑148 (MQ=255)
                   aaaaGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTCTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCgacgac                                                                                                                               >  4:199046/1‑148 (MQ=255)
                     aaGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGc                                                                                                                             <  2:75716/148‑1 (MQ=255)
                         tGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCTccc                                                                                                                         <  5:145200/148‑1 (MQ=255)
                           ttCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGt                                                                                                                       <  7:203121/148‑1 (MQ=255)
                            tCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACCTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGt                                                                                                                       >  6:75111/1‑147 (MQ=255)
                              aCGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGttt                                                                                                                     >  7:205900/1‑147 (MQ=255)
                                 aTACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGc                                                                                                                             >  1:187748/1‑136 (MQ=255)
                                 aTACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGc                                                                                                                             <  2:187748/136‑1 (MQ=255)
                                  tACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCaaa                                                                                                                 >  7:79816/1‑147 (MQ=255)
                                  tACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAAt                                                                                                                <  1:302073/148‑1 (MQ=255)
                                      cTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTg                                                                                                            >  1:157917/1‑148 (MQ=255)
                                           cTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGct                                                                                                       <  2:302279/148‑1 (MQ=255)
                                              cGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGctcc                                                                                                     <  4:146929/147‑1 (MQ=255)
                                               ggCAAAATAACCAGGAAAAATCCCGCTATTTCCTCCCGTTTTAAGCCAGCTCCCCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCCCCCGTTTCAAATACTGGTGctcctc                                                                                                   >  2:46771/1‑148 (MQ=255)
                                                gCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGCCGACGCACCCGTTTCAAATACTGGTGCtcctcc                                                                                                  >  1:185461/1‑148 (MQ=255)
                                                gCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCtcctcc                                                                                                  <  8:205900/148‑1 (MQ=255)
                                                gCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTCTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCtcctcc                                                                                                  <  3:199046/148‑1 (MQ=255)
                                                    aaTAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCcgacga                                                                                                                                >  5:274121/1‑114 (MQ=255)
                                                    aaTAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCcgacga                                                                                                                                <  6:274121/114‑1 (MQ=255)
                                                           aGGAAAACTACAGGTATTTCCTCACGTTTTAAGCCAGCTCACGCTTCCCAAAAGGCACAAAGCGCACGGCCTCCAGGGTATAGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACAGGTGCTCCTCCCCTACGGGt                                                                                         >  1:347250/1‑146 (MQ=255)
                                                              aaaaaTCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGCGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGAc                                                                                    >  6:284733/1‑148 (MQ=255)
                                                               aaaaTCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACg                                                                                   <  5:78191/148‑1 (MQ=255)
                                                                  aTCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACgaga                                                                                <  5:284733/148‑1 (MQ=255)
                                                                   tCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGaa                                                                               <  2:157917/148‑1 (MQ=255)
                                                                        gTATTTCCTCACGTTTTAAGCCAGCTCTCCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCg                                                                          <  6:15920/148‑1 (MQ=255)
                                                                         tATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGc                                                                         <  4:50375/148‑1 (MQ=255)
                                                                                 cACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAAt                                                                                                                                               >  8:252113/1‑70 (MQ=255)
                                                                                 cACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAAt                                                                                                                                               <  7:252113/70‑1 (MQ=255)
                                                                                  aCGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTcc                                                                 <  1:81780/147‑1 (MQ=255)
                                                                                     ttttAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAGTAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTgg                                                                                                           >  1:295657/1‑102 (MQ=255)
                                                                                     ttttAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAGTAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTgg                                                                                                           >  1:290291/1‑102 (MQ=255)
                                                                                     ttttAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAGTAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTgg                                                                                                           <  2:295657/102‑1 (MQ=255)
                                                                                     ttttAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAGTAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTgg                                                                                                           <  2:290291/102‑1 (MQ=255)
                                                                                          aGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTc                                                        <  8:392671/148‑1 (MQ=255)
                                                                                           gCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCa                                                       <  5:75111/148‑1 (MQ=255)
                                                                                                 tCACCCTGCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAgcg                                                 >  7:198589/1‑148 (MQ=255)
                                                                                                        tCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTg                                           >  1:138042/1‑147 (MQ=255)
                                                                                                        tCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTg                                           >  8:351199/1‑147 (MQ=255)
                                                                                                          aCTAAAGGGACAAAGCGCACTGCCTCCACGGTATCGATAATAAATTCGCCTCCCAGACGACGCCCCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGATAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGa                                         <  3:91896/147‑1 (MQ=255)
                                                                                                          aCTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGa                                         <  1:46771/147‑1 (MQ=255)
                                                                                                            tAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCcgacg                                                                                                                                 <  4:15019/57‑1 (MQ=255)
                                                                                                            tAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCcgacg                                                                                                                                 >  3:15019/1‑57 (MQ=255)
                                                                                                             aaaGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAAttt                                     >  5:57911/1‑148 (MQ=255)
                                                                                                              aaGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAATTTc                                    >  7:79136/1‑148 (MQ=255)
                                                                                                                 ggACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAATTTCcgg                                 >  2:380965/1‑148 (MQ=255)
                                                                                                                 ggACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAATTTCcgg                                 <  6:57911/148‑1 (MQ=255)
                                                                                                                  gACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTc                                                        >  4:159839/1‑124 (MQ=255)
                                                                                                                  gACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTc                                                        <  3:159839/124‑1 (MQ=255)
                                                                                                                  gACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAATTTCcggc                                <  2:138042/148‑1 (MQ=255)
                                                                                                                     aaaGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAATTTCCGGCGGt                             <  7:23465/148‑1 (MQ=255)
                                                                                                                      aaGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAATTTCCGGCGGTg                            >  3:343378/1‑148 (MQ=255)
                                                                                                                         cgcACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTATGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAATTTCCGGCGGTgccg                         >  5:224740/1‑148 (MQ=255)
                                                                                                                           cACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAATTTCCGGCGGTgccgcc                       <  8:297100/148‑1 (MQ=255)
                                                                                                                                   ccACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCg                                                                    <  8:156035/95‑1 (MQ=255)
                                                                                                                                   ccACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCg                                                                    >  7:156035/1‑95 (MQ=255)
                                                                                                                                   ccACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAATTTCCGGCGGTGCCGCCGTAACAAt               >  4:131292/1‑148 (MQ=255)
                                                                                                                                   ccACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAATTTCCGGCGGTGCCGCCGTAACAAt               <  1:176998/148‑1 (MQ=255)
                                                                                                                                    gaCGGTATGGATTATTAATTGGGGTCCCCGAGGACGCAACGGTTTGAAATACTGGTGCTCCTCCCCTACGGGTAAGACCAGAATCCCGCCTTGGTCCAGCTGCGTCATTTGGGCAGTTGGAATTTCCGGCGGTGCCGCCGTAACAAtc              <  2:347250/147‑2 (MQ=255)
                                                                                                                                    cACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAATTTCCGGCGGTGCCGCCGTAACaa                <  4:343378/146‑1 (MQ=255)
                                                                                                                                    cACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAATTTCCGGCGGTGCCGCCGTAACAAt               <  2:337978/147‑1 (MQ=255)
                                                                                                                                     aCGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAATTTCCGGCGGTGCCGCCGTAACAATGa             <  8:48350/148‑1 (MQ=255)
                                                                                                                                     aCGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAATTTCcggcg                               >  4:248273/1‑129 (MQ=255)
                                                                                                                                     aCGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAATTTCcggcg                               <  3:248273/129‑1 (MQ=255)
                                                                                                                                      cGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAATTTCCGGCGGTGCCGCCGTAACAATGAt            <  3:131292/148‑1 (MQ=255)
                                                                                                                                      cGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAATTTCCGGCGGTGCCGCCGTAACAATGAt            <  1:380965/148‑1 (MQ=255)
                                                                                                                                       ggTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCGCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAATTTCCGGCGGTGCCGCCGTAACAATGATa           <  4:199624/148‑1 (MQ=255)
                                                                                                                                        gTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAATTTCCGGCGGTGCCGCCGTAACAATGATAg          >  3:306426/1‑148 (MQ=255)
                                                                                                                                         tATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGa                                         <  7:187296/116‑1 (MQ=255)
                                                                                                                                         tATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGa                                         >  8:187296/1‑116 (MQ=255)
                                                                                                                                                aataaATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAATTTCcggcg                               <  2:141607/119‑1 (MQ=255)
                                                                                                                                                aataaATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAATTTCcggcg                               >  1:141607/1‑119 (MQ=255)
                                                                                                                                                aataaATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAATTTCcggc                                <  4:251520/118‑1 (MQ=255)
                                                                                                                                                aataaATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAATTTCcggc                                >  3:251520/1‑118 (MQ=255)
                                                                                                                                                aataaATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAATTTCCGGCGGTGCCGCCGTAACAATGATAGCGTCAAAc  >  7:145440/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GCACCCTACTGAACACTGGAAAAGATGTTCACGATACGCTGACCTGCGGCAAAATAACCAGGAAAAATCCAGGTATTTCCTCACGTTTTAAGCCAGCTCACCCTTCACTAAAGGGACAAAGCGCACGGCCTCCACGGTATCGATAATAAATTCGCCTCCCCGACGACGCACCCGTTTCAAATACTGGTGCTCCTCCCCTACGGGTAAGACGAGAATCCCGCCTTCGTCCAGCTGCGTCATTAGCGCAGTTGGAATTTCCGGCGGTGCCGCCGTAACAATGATAGCGTCAAAC  >  NC_000913/2868806‑2869097

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