Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,407,898 A→C 5.6% Q113P (CAG→CCG) ‡ alaA → glutamate‑pyruvate aminotransferase; glutamic‑pyruvic transaminase (GPT); alanine transaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,407,8980AC 5.6% 325.0 / 10.2 89Q113P (CAG→CCG) ‡alaAglutamate‑pyruvate aminotransferase; glutamic‑pyruvic transaminase (GPT); alanine transaminase
Reads supporting (aligned to +/- strand):  ref base A (41/43);  new base C (3/2);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

CCGCTCAAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGG  >  NC_000913/2407751‑2408022
                                                                                                                                                   |                                                                                                                             
ccGCTCAAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCa                                                                                                                               <  3:135812/148‑1 (MQ=255)
      aaGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCAtt                                                                                                                         <  2:782828/148‑1 (MQ=255)
       aGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTg                                                                                                                        >  2:918739/1‑148 (MQ=255)
        gggTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGc                                                                                                                      <  1:925025/148‑1 (MQ=255)
         ggTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCt                                                                                                                     <  4:948699/148‑1 (MQ=255)
            aTTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGaa                                                                                                                  >  1:285260/1‑148 (MQ=255)
              tGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACa                                                                                                                >  2:634633/1‑148 (MQ=255)
                 gATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCg                                                                                                             <  1:634633/148‑1 (MQ=255)
                    tCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTTTCGTTCAGGCAATGCAGGCATTGCTGAACAGCgggg                                                                                                          <  1:918739/148‑1 (MQ=255)
                      cAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAgg                                                                                                                             >  2:161291/1‑128 (MQ=255)
                      cAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAgg                                                                                                                             <  1:161291/128‑1 (MQ=255)
                      cAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGAc                                                                                                        >  3:916462/1‑148 (MQ=255)
                       aaaGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCgggg                                                                                                          >  1:660783/1‑145 (MQ=255)
                       aaaGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCgggg                                                                                                          <  2:660783/145‑1 (MQ=255)
                            tCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCgg                                                                                                            <  4:817019/138‑1 (MQ=255)
                            tCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCgg                                                                                                            >  3:817019/1‑138 (MQ=255)
                                    ccGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGtt                                                                                            <  1:798724/146‑1 (MQ=255)
                                    ccGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGtt                                                                                            >  2:798724/1‑146 (MQ=255)
                                    ccGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTcc                                                                                          <  2:285260/148‑1 (MQ=255)
                                     cgcgcgTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGa                                                                                                                   >  2:903919/1‑122 (MQ=255)
                                     cgcgcgTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGa                                                                                                                   <  1:903919/122‑1 (MQ=255)
                                       cgcgTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCAt                                                                                                                          >  4:368034/1‑114 (MQ=255)
                                       cgcgTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCAt                                                                                                                          <  3:368034/114‑1 (MQ=255)
                                       cgcgTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGc                                                                                       >  4:980344/1‑148 (MQ=255)
                                                ccATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCgggg                                                                                                          >  4:484758/1‑120 (MQ=255)
                                                ccATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCgggg                                                                                                          <  3:484758/120‑1 (MQ=255)
                                                  atcatGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTGCCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGa                                                                                                         >  4:662536/1‑119 (MQ=255)
                                                  atcatGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGa                                                                                                         <  3:662536/119‑1 (MQ=255)
                                                   tcatGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCAcc                                                                                    <  2:457023/139‑1 (MQ=255)
                                                   tcatGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCAcc                                                                                    >  1:457023/1‑139 (MQ=255)
                                                    catGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGt                                                                                                 <  4:692290/125‑1 (MQ=255)
                                                    catGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGt                                                                                                 >  3:692290/1‑125 (MQ=255)
                                                    catGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCa                                                                          >  1:528578/1‑148 (MQ=255)
                                                      tGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCAct                                                                        <  2:528578/148‑1 (MQ=255)
                                                      tGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCAct                                                                        <  4:916462/148‑1 (MQ=255)
                                                          gcaCTACCAGGCTCGTGGCATGCGTGATGTTACTGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCCTGCATTGcctga                                                                                                                  >  3:376599/1‑98 (MQ=255)
                                                               aCCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACcgc                                                               >  4:480518/1‑148 (MQ=255)
                                                                ccAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACcgcg                                                              >  1:833697/1‑148 (MQ=255)
                                                                 cAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGc                                                                                                              >  3:678663/1‑99 (MQ=255)
                                                                 cAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGaa                                                                                                     <  3:833368/108‑1 (MQ=255)
                                                                 cAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGaa                                                                                                     >  4:833368/1‑108 (MQ=255)
                                                                 cAGGCTCGTGGCATGCGTGATGGTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGc                                                                                                              <  4:678663/99‑1 (MQ=255)
                                                                  aGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTa                                                                              >  3:383885/1‑130 (MQ=255)
                                                                  aGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTa                                                                              <  4:383885/130‑1 (MQ=255)
                                                                   ggCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGt                                                                                             >  4:979020/1‑114 (MQ=255)
                                                                   ggCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGt                                                                                             <  3:979020/114‑1 (MQ=255)
                                                                   ggCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCgcggcg                                                           <  3:980344/148‑1 (MQ=255)
                                                                     cTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGt                                                         <  3:480518/148‑1 (MQ=255)
                                                                        gTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTc                                                      <  2:833697/148‑1 (MQ=255)
                                                                         tGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTg                                                                                               >  2:853366/1‑106 (MQ=255)
                                                                         tGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTg                                                                                               <  1:853366/106‑1 (MQ=255)
                                                                              tGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGCAATGCCTGCATTGcct                                                                                                                    >  2:1081214/1‑77 (MQ=38)
                                                                              tGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGCAATGCCTGCATTGcct                                                                                                                    <  1:1081214/79‑3 (MQ=38)
                                                                                cGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTGAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACcgcg                                                              >  2:335732/1‑132 (MQ=255)
                                                                                cGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACcgcg                                                              <  1:335732/132‑1 (MQ=255)
                                                                                   gATGTTACCGTGGCAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGtt                                                                                                <  4:965681/95‑1 (MQ=255)
                                                                                   gATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGATGAACAGCGGGGACGAAATGtt                                                                                                >  3:965681/1‑95 (MQ=255)
                                                                                        tACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAg                                                                                                               <  3:956355/75‑1 (MQ=255)
                                                                                        tACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAg                                                                                                               >  4:956355/1‑75 (MQ=255)
                                                                                              ggAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATg                                                                                                  <  3:444612/82‑1 (MQ=255)
                                                                                              ggAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATg                                                                                                  >  4:444612/1‑82 (MQ=255)
                                                                                                    tatTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTg                                                                                                                    <  4:9187/58‑1 (MQ=255)
                                                                                                    tatTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTg                                                                                                                    >  3:9187/1‑58 (MQ=255)
                                                                                                      tttACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCt                        >  1:1036859/1‑148 (MQ=255)
                                                                                                            tCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAg                                                                                                               <  1:759032/55‑1 (MQ=255)
                                                                                                            tCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAg                                                                                                               >  2:759032/1‑55 (MQ=255)
                                                                                                             cGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTc                                                >  1:1026254/1‑117 (MQ=255)
                                                                                                             cGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTc                                                <  2:1026254/117‑1 (MQ=255)
                                                                                                                cAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTAc                                                                             >  3:932623/1‑85 (MQ=255)
                                                                                                                cAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTAc                                                                             <  4:932623/85‑1 (MQ=255)
                                                                                                                    ggTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGtt                                                                                                <  2:1030963/62‑1 (MQ=255)
                                                                                                                    ggTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGtt                                                                                                >  1:1030963/1‑62 (MQ=255)
                                                                                                                      tgtATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCActct                                                                      <  3:658976/86‑1 (MQ=255)
                                                                                                                      tgtATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCActct                                                                      >  4:658976/1‑86 (MQ=255)
                                                                                                                           cGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTTGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTg                     <  2:997335/130‑1 (MQ=255)
                                                                                                                           cGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTg                     >  1:997335/1‑130 (MQ=255)
                                                                                                                              aGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGtt                                                        <  3:428711/92‑1 (MQ=255)
                                                                                                                              aGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGtt                                                        >  4:428711/1‑92 (MQ=255)
                                                                                                                                 ttATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCg                                                     <  1:683307/92‑1 (MQ=255)
                                                                                                                                 ttATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCg                                                     >  2:683307/1‑92 (MQ=255)
                                                                                                                                    tCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCAc                                                                         >  4:270382/1‑69 (MQ=255)
                                                                                                                                    tCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCAc                                                                         <  3:270382/69‑1 (MQ=255)
                                                                                                                                    ctgTTCAGCAATGCCTGCATTGCCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTg                     >  2:958992/3‑121 (MQ=255)
                                                                                                                                    ctgTTCAGCAATGCCTGCATTGCCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTg                     <  1:958992/119‑1 (MQ=255)
                                                                                                                                     cGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCg                                                     >  4:411468/1‑88 (MQ=255)
                                                                                                                                     cGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCg                                                     <  3:411468/88‑1 (MQ=255)
                                                                                                                                           ggCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTgg  >  1:924774/1‑133 (MQ=255)
                                                                                                                                           ggCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTgg  <  2:924774/133‑1 (MQ=255)
                                                                                                                                             cAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGc                                    <  4:930077/97‑1 (MQ=255)
                                                                                                                                             cAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGc                                    >  3:930077/1‑97 (MQ=255)
                                                                                                                                                   |                                                                                                                             
CCGCTCAAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGG  >  NC_000913/2407751‑2408022

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