Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,863,313 T→C A181A (GCT→GCC gapA → glyceraldehyde‑3‑phosphate dehydrogenase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,863,3130TC98.2% 399.0 / ‑4.0 108A181A (GCT→GCCgapAglyceraldehyde‑3‑phosphate dehydrogenase A
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (62/44);  total (64/44)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

AGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGTAGGTAAAGTACTGCCAGAACTGAATGGCAAACTGACTGGTATGGCG  >  NC_000913/1863166‑1863460
                                                                                                                                                   |                                                                                                                                                   
aGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACt                                                                                                                                                  >  2:626106/1‑151 (MQ=255)
  gcgcTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTc                                                                                                                                                 >  2:1731885/1‑150 (MQ=255)
          ttCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCg                                                                                                                                        >  1:42530/1‑151 (MQ=255)
          ttCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCg                                                                                                                                        >  2:549798/1‑151 (MQ=255)
          ttCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCg                                                                                                                                        >  2:561394/1‑151 (MQ=255)
          ttCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCg                                                                                                                                        >  2:1467058/1‑151 (MQ=255)
             gACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGtt                                                                                                                                      >  2:1119192/1‑150 (MQ=255)
              aaaaaTATGCGGGCCAGGACTTCGTTTCCAACGCTCCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAATTTATCAACGATAACTCCGGCATCATCAAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAAcc                                                                                                                                         <  1:1640450/144‑1 (MQ=255)
              aCAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCg                                                                                                                                        <  2:526138/147‑1 (MQ=255)
               cAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGAt                                                                                                                                   >  2:685644/1‑151 (MQ=255)
                aaaTATGCTGGCCAGGACATCGTTTCCAACGCTTCCTCCCCCACCAACTGCCTGGCCCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCCCGCTACTACCGCCACTCAGAAAACCGtt                                                                                                                                      <  2:1438970/147‑1 (MQ=255)
                 aaTATGCTAGCCAGGACATCGTTTCCAACGCTTACTTCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATgg                                                                                                                                 >  1:954624/1‑151 (MQ=255)
                  atatGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGc                                                                                                                                >  1:456389/1‑151 (MQ=255)
                  atatGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGc                                                                                                                                >  3:308474/1‑151 (MQ=255)
                      gCTGGCCAGGACATCGTTCCCAACGCTTCCTGCACCACCAACTGCCTGGCCCCGCTGTCTAAAGTTATCACCGATAACTTCGGCATCATCGAAGGTTTAAGGACCCCCGTTCACTCTACTACCGCCACTCAGAAAACCGTTGATGGCCCg                                                                                                                             <  1:1060213/150‑1 (MQ=255)
                           ccAGGACATCGTTTCCACCGCTTCCTCCCCCCCCACCTGCTTGGCTCCGCTGGCTAAAGTTATCAACAATACCTTCGGCATCATCGAAGTTCTGATGCCCACCGTTCACGCTCCTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTca                                                                                                                        <  2:58218/150‑1 (MQ=255)
                             aGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTcaca                                                                                                                      >  2:741465/1‑150 (MQ=255)
                             aGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTcaca                                                                                                                      >  2:753644/1‑150 (MQ=255)
                             aGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTcaca                                                                                                                      >  2:999802/1‑150 (MQ=255)
                              ggACATCTTTTCCAACGCTTCCTGCACCACCACCTGCCTGGCTCCGCTGGCAAAATTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCCCCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACaaa                                                                                                                    <  4:483762/151‑1 (MQ=255)
                              ggACATCGTTTCAAACGCTCCCTGCCCCCCCAACCGCCTGGCTCCTCTGCTTAAAGTTATAAACAATAATTTCGGAATCAAGGAAGGTCTGACGACCACGTTCCCCTCTCTTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACaa                                                                                                                     <  1:1074696/150‑1 (MQ=255)
                               gACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACaa                                                                                                                     >  2:1656548/1‑149 (MQ=255)
                                aCATTGTTTCCACCGCTTCCTTCCCCCCCATCTGCCGGGCCCCGCTGGCTAAATTTATCAACGATACCTTCGGCATCATGAAAGGTCTGATCCCCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGa                                                                                                                  <  1:246901/151‑1 (MQ=255)
                                aCATCGTTTCCAACGCTTCCTGCACCCCCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCCCCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGa                                                                                                                  <  2:1047424/151‑1 (MQ=255)
                                     gTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGgc                                                                                                             >  1:445931/1‑151 (MQ=255)
                                      tttCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGgcg                                                                                                            >  2:38012/1‑151 (MQ=255)
                                         caaaCGCTCCCGGCGCCACCAACGGCCGTGTCCCGCTGCCTAAATTAATCAACTATAGTTTCGTCGTAATCGAGGGTAGGGTGCCCACCTTGCACTCAACTTCCGCCACTCAGAAAACCGTGGATGGCCCGTCTCACAAAGACTGgcg                                                                                                            <  1:1293944/146‑1 (MQ=255)
                                           aaCGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCgcggcg                                                                                                       >  1:1247399/1‑151 (MQ=255)
                                             cGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTCCTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGcggcgg                                                                                                      <  3:107637/150‑1 (MQ=255)
                                              gCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGcc                                                                                                    >  1:1321253/1‑151 (MQ=255)
                                              gCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGcc                                                                                                    >  1:1295059/1‑151 (MQ=255)
                                               cTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCcg                                                                                                   >  2:175587/1‑151 (MQ=255)
                                                  cctccccccccaACTGCCTGCCTCCGCTGGCTAAATTAATCAACGATAACTTGGGCACCATCGAAGGTCTGATGCCCCCCTTCCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCcgcg                                                                                                 <  2:1650831/141‑1 (MQ=255)
                                                  ccTGCCCCCCCAACTGCGTGGCCCCGGGGGTTAAGTTTATCAACGATATCTTCGGCGTCATCAAAGGTGTGATGACCACCGTTCCCGCTACTCCCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCcgc                                                                                                  <  2:1561577/149‑1 (MQ=255)
                                                    tGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCgcggcg                                                                                              >  4:561395/1‑151 (MQ=255)
                                                        acaccaACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGACGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCt                                                                                            >  1:379263/2‑149 (MQ=255)
                                                         caccaACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTcc                                                                                         >  1:1154226/1‑151 (MQ=255)
                                                           caacctgcagggCCCCGCTGGCAAAAGTTATCAACAATAACTTCGGAATCATCGAAGGTCTGATGCCCCCCGTCCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCa                                                                                       <  2:1120332/141‑1 (MQ=255)
                                                             aaCTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTCCTACCGCCACTCAGAAAACCGTTGATGGCCCg                                                                                                                             <  4:52186/111‑1 (MQ=255)
                                                             aaCTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGTCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTacc                                                                                        >  2:202402/1‑148 (MQ=255)
                                                             aaCTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCg                                                                                                                             >  3:52187/1‑111 (MQ=255)
                                                             aaCTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTccc                                                                                        >  1:779130/1‑148 (MQ=255)
                                                                   cTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAAcatc                                                                                >  2:937838/1‑150 (MQ=255)
                                                                     ggCTCCGTTGGCTAAAGTTATCAACGATAACTTGGGCATCATCGAAGTTCTGATGCCCCCCTTTCACGCTCCTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAAcatcat                                                                              <  2:698262/150‑1 (MQ=255)
                                                                      gCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATcc                                                                            >  2:378323/1‑151 (MQ=255)
                                                                      gCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATcc                                                                            >  1:612956/1‑151 (MQ=255)
                                                                        tCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCg                                                                          >  2:240691/1‑151 (MQ=255)
                                                                         ccGCTGGCTAAATTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATccc                                                                           <  3:375263/149‑1 (MQ=255)
                                                                         ccGCTGGCTAAAGTTATCAGCATTATCTGCGTAATCAACGACGGTTTGATCACCCCCGCGCTCGCTCCAACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGt                                                                         <  1:652331/151‑1 (MQ=255)
                                                                          cGCTGGCAAAAGTTATCAACAATATCTTCGCCATCATCAAAGCTCTTATCACCCCCTTTCACCCTCCTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATcc                                                                            <  4:112314/147‑1 (MQ=255)
                                                                             tGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCctc                                                                     >  2:1532158/1‑151 (MQ=255)
                                                                             tGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCctc                                                                     >  1:338166/1‑151 (MQ=255)
                                                                             tGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCct                                                                      >  3:50151/1‑150 (MQ=255)
                                                                              ggCTAAATTTATAAACGATACCTTCGGCATCATCGAAGGTCTGATACCCCCCTTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCATCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCctct                                                                    <  2:1078398/151‑1 (MQ=255)
                                                                                cTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTa                                                                   >  1:770626/1‑150 (MQ=255)
                                                                                 tAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTAcc                                                                 >  1:219197/1‑151 (MQ=255)
                                                                                  aaaGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCCCCGTTCACGCTCCTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCg                                                                <  2:102558/151‑1 (MQ=255)
                                                                                  aaaGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGgcgc                                                                                                           >  4:105639/1‑108 (MQ=255)
                                                                                  aaaGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGgcgc                                                                                                           <  3:105643/108‑1 (MQ=255)
                                                                                  aaaGTTACAAACAATAACTTCGGCATCACCGAAGGTCTGATGACCCCCTTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTAc                                                                  <  2:578755/149‑1 (MQ=255)
                                                                                    aGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGt                                                              >  2:1309083/1‑151 (MQ=255)
                                                                                       taccAACGTTAACTTCGGCATCATGAAAGGTCTGATGCCCTCCGTTCACCCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGt                                                              <  1:1567149/145‑1 (MQ=255)
                                                                                        aTCAACGATAACTTCGGCATCATCGAAGGTCTGATGCCCCCCTTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGtgct                                                           <  4:356234/150‑1 (MQ=255)
                                                                                        aTCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGtgct                                                           >  1:1474572/1‑150 (MQ=255)
                                                                                         tCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGtgctgc                                                         >  2:766882/1‑151 (MQ=255)
                                                                                           aaCGATACTTTCGGCTTCATCGAAGGTCTGATGCCCCCCGTCCCCGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGtgc                                                            <  1:977591/146‑1 (MQ=255)
                                                                                           aaCGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTgctgct                                                        >  2:204960/1‑150 (MQ=255)
                                                                                           aaCGATAACTTCGCCATCATGGAAGGTCTGAGGACCCCCGTTCCCGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTa                                                       <  3:427899/151‑1 (MQ=255)
                                                                                             cGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCCCGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTaaa                                                     <  2:315761/151‑1 (MQ=255)
                                                                                              gATAACTTCGGCATCATCGAAGGGCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATCGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGCGCCGCTaa                                                      >  1:639649/1‑149 (MQ=255)
                                                                                               aTAAGTTCGGAAACATAGAAGGTTTGAGGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGc                                                   <  2:1169052/151‑1 (MQ=255)
                                                                                               aTAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGc                                                   >  4:500431/1‑151 (MQ=255)
                                                                                                tAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCt                                                  >  2:237469/1‑151 (MQ=255)
                                                                                                 aaCTCCGGCATCATCGAAGGTCTGATGCCCCCCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTaaa                                                     <  1:1668121/147‑1 (MQ=255)
                                                                                                    ttCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGt                                                >  2:1521243/1‑149 (MQ=255)
                                                                                                    ttCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGt                                                >  1:804487/1‑149 (MQ=255)
                                                                                                      cGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGTAgg                                             >  2:977326/1‑150 (MQ=255)
                                                                                                      cGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGTAGGt                                            >  1:224793/1‑151 (MQ=255)
                                                                                                       ggcacccgcGAAGTTCTGAAGACCCCCGTCCACGCTACTCCCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGTAGGTa                                           <  1:198397/143‑1 (MQ=255)
                                                                                                         catcatCGAAGTTCTAATGCCCCCCGTTCCCGCTCCTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGTAGGTaa                                          <  1:438845/150‑1 (MQ=255)
                                                                                                          atcatcGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGTAGGTaa                                          >  1:1474769/1‑149 (MQ=255)
                                                                                                          aacatcGAAGTTCTGAGCCCCACCGTTCGCTGTTCAACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTGCCGGTGCTGCTAAACCTGTAGGTAAAg                                        <  1:564895/149‑1 (MQ=255)
                                                                                                           tcatcGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGTAGGTAAAGt                                       >  1:1562877/1‑151 (MQ=255)
                                                                                                             atcGAAGGTCTGATGACCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGTAGGTAAAGTAc                                     >  2:1119946/1‑151 (MQ=255)
                                                                                                                 aaGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGTAGGTAAAGTACTGc                                  >  3:189623/1‑150 (MQ=255)
                                                                                                                   ggTTTGATGCCCCCCTTCCACGCTCCTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGTAGGTAAAGTACTg                                   <  1:173469/147‑1 (MQ=255)
                                                                                                                       tgatgaCCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGTAGGTAAAGTACTGCCAGAACt                           >  1:596066/1‑151 (MQ=255)
                                                                                                                        gatgaCCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGTAGGTAAAGTACTGCCAGAACTg                          <  3:209295/151‑1 (MQ=255)
                                                                                                                        gatgaCCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGTAGGTAAAGTACTGCCAGAACTg                          >  2:1376796/1‑151 (MQ=255)
                                                                                                                           gaCCACCGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGTAGGTAAAGTACTGCCAGAACTGAAt                       >  1:251977/1‑151 (MQ=255)
                                                                                                                             ccaccGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGTAGGTAAAGTACTGCCAGAACTGAATgg                     >  2:257167/1‑151 (MQ=255)
                                                                                                                              caccGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGTAGGTAAAGTACTGCCAGAACTGAATGGc                    >  4:251263/1‑151 (MQ=255)
                                                                                                                                 cGTTCACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGTAGGTAAAGTACTGCCAGAACTGAATGGCaa                  >  3:514853/1‑150 (MQ=255)
                                                                                                                                    tCACGCTACTACAGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGTAGGTAAAGTACTGCCAGAACTGACTGGCAAact               >  3:327032/1‑150 (MQ=255)
                                                                                                                                     cACGCTACTACCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGTAGGTAAAGTACTGCCAGAACTGAATGGCAAactg              <  4:91050/150‑1 (MQ=255)
                                                                                                                                          tactacCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGTAGGTAAAGTACTGCCAGAACTGAATGGCAAACTGACTGGt        >  2:911230/1‑151 (MQ=255)
                                                                                                                                          tactacCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGTAGGTAAAGTACTGCCAGAACTGAATGGCAAACTGACTGGt        >  2:1622337/1‑151 (MQ=255)
                                                                                                                                             tacCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGTAGGTAAAGTACTGCCAGAACTGAATGGCAAACTGACTGGTATg     >  1:677447/1‑151 (MQ=255)
                                                                                                                                              acCGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGTAGGTAAAGTACTGCCAGAACTGAATGGCAAACTGACTGGTATg     <  2:456464/150‑1 (MQ=255)
                                                                                                                                                cGCCACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGTAGGTAAAGTACTGCCAGAACTGAATGGCAAACTGACTGGTATGgcg  >  3:376971/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGGCGCTAACTTCGACAAATATGCTGGCCAGGACATCGTTTCCAACGCTTCCTGCACCACCAACTGCCTGGCTCCGCTGGCTAAAGTTATCAACGATAACTTCGGCATCATCGAAGGTCTGATGACCACCGTTCACGCTACTACCGCTACTCAGAAAACCGTTGATGGCCCGTCTCACAAAGACTGGCGCGGCGGCCGCGGCGCTTCCCAGAACATCATCCCGTCCTCTACCGGTGCTGCTAAAGCTGTAGGTAAAGTACTGCCAGAACTGAATGGCAAACTGACTGGTATGGCG  >  NC_000913/1863166‑1863460

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