Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,407,890 Δ1 bp 6.4% coding (330/1218 nt) 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,8900G. 6.4% 330.5 / 16.5 94coding (330/1218 nt)alaAglutamate‑pyruvate aminotransferase; glutamic‑pyruvic transaminase (GPT); alanine transaminase
Reads supporting (aligned to +/- strand):  ref base G (43/45);  new base . (3/3);  total (46/48)
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.47e-01

TACCGCTCAAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGG  >  NC_000913/2407749‑2408022
                                                                                                                                             |                                                                                                                                     
tACCGCTCAAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGAAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATg                                                                                                                                 >  2:925025/1‑148 (MQ=255)
 aCCGCTCGAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGc                                                                                                                                <  3:23460/148‑1 (MQ=255)
  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                                                                                                                             <  1:161291/128‑1 (MQ=255)
                        cAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAgg                                                                                                                             >  2:161291/1‑128 (MQ=255)
                        cAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGAc                                                                                                        >  3:916462/1‑148 (MQ=255)
                         aaaGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCgggg                                                                                                          >  1:660783/1‑145 (MQ=255)
                         aaaGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCgggg                                                                                                          <  2:660783/145‑1 (MQ=255)
                         aaaGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGCAATGcc                                                                                                                               <  4:923227/124‑2 (MQ=255)
                         aaaGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGCAATGcc                                                                                                                               >  3:923227/1‑123 (MQ=255)
                              tCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGcc                                                                                                                               <  2:321552/120‑2 (MQ=255)
                              tCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCgg                                                                                                            <  4:817019/138‑1 (MQ=255)
                              tCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCgg                                                                                                            >  3:817019/1‑138 (MQ=255)
                              tCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAAGGcc                                                                                                                               >  1:321552/1‑119 (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                                                                                                          <  3:484758/120‑1 (MQ=255)
                                                  ccATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCgggg                                                                                                          >  4:484758/1‑120 (MQ=255)
                                                   catcatGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCa                                                                                                                                    <  3:821238/94‑1 (MQ=255)
                                                   catcatGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCa                                                                                                                                    >  4:821238/1‑94 (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                                                                                             <  3:979020/114‑1 (MQ=255)
                                                                     ggCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGt                                                                                             >  4:979020/1‑114 (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                                                                                                                    <  1:1081214/79‑3 (MQ=38)
                                                                                tGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGCAATGCCTGCATTGcct                                                                                                                    >  2:1081214/1‑77 (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                                                                                                               >  4:956355/1‑75 (MQ=255)
                                                                                          tACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAg                                                                                                               <  3:956355/75‑1 (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                                                <  2:1026254/117‑1 (MQ=255)
                                                                                                               cGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTc                                                >  1:1026254/1‑117 (MQ=255)
                                                                                                                  cAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTAc                                                                             >  3:932623/1‑85 (MQ=255)
                                                                                                                  cAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTAc                                                                             <  4:932623/85‑1 (MQ=255)
                                                                                                                      ggTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGtt                                                                                                >  1:1030963/1‑62 (MQ=255)
                                                                                                                      ggTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGtt                                                                                                <  2:1030963/62‑1 (MQ=255)
                                                                                                                        tgtATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCActct                                                                      >  4:658976/1‑86 (MQ=255)
                                                                                                                        tgtATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCActct                                                                      <  3:658976/86‑1 (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                                                     >  2:683307/1‑92 (MQ=255)
                                                                                                                                   ttATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCg                                                     <  1:683307/92‑1 (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                                                     <  3:411468/88‑1 (MQ=255)
                                                                                                                                       cGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCg                                                     >  4:411468/1‑88 (MQ=255)
                                                                                                                                             ggCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTgg  <  2:924774/133‑1 (MQ=255)
                                                                                                                                             ggCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTgg  >  1:924774/1‑133 (MQ=255)
                                                                                                                                             |                                                                                                                                     
TACCGCTCAAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTGCTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGG  >  NC_000913/2407749‑2408022

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